MICROCOPY  RESOLUTION  TEST  CHART 

RATIONAL  BUREAU  Of  STANDARDS  1963-A 


I 

i 

I 


> 

' 

» 


Unclassified 


A'-x. '  •  X<x 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

Unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

SSS-R-87-8454 


lb.  RESTRICTIVE  MARKINGS 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release; 
distribution  is  unlimited. 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

AFObK  TK-  a  7  “  13yi 


6a.  NAME  OF  PERFORMING  ORGANIZATION  [6b.  OFFICE  SYMBOL  I  7a.  NAME  OF  MONITORING  ORGANIZATION 
s-rnRRD  (  I  Cf  applicable)  1 


S-CUBED  y 


6c  ADDRESS  {City,  Stitt,  and  ZIP  Cod*) 

P  0  Box  1620^1 
La  Jolla,  CA  92038 


8a.  NAME  OF  FUNDING  /SPONSORING 
ORGANIZATION 


AFOSR/NA 


7b.  ADDRESS  {City,  Stite,  and  ZIP  Cod*) 

Building  410,  Bolling  AFB  DC 
20332-6448 


AFOSR/NA 


1 1 .  TITLE  (Include  Security  Cassil 

(U)  Development  oi 


10.  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  PROJECT  TASK 

ELEMENT  NO.  NO.  NO, 

61102F  2302  C2 


anced/Gorostitutive  Models  for  Plain  and  Reinforced 


Concrete  V/ 


12  GRA>NHegeTmier,  H  E  Read,  K  C  valanis,  H  Murakami 


13a.  TYPE  OF  REPORT 

Annual 


mmm 


IA.  DATE  OF  REPORT  {Year,  Month.  Day) 

'1987  AugubL  Mr  m^rch 


COSAT1  CODES 


GROUP  SUB-GROUP 


^ 'WffSfe  £  necessary  and  identify  by  Nock  number) 

Continuous  Damage  Model;  Dowel  Action; 
-Evelrralrrtm  -Equation-;-  Endochronic  Theory*  M±wte*«ee 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  report  summarizes  research  conducted  by  S-Cubed  to  develop  advance 
constitutive  models  of  plain  and  reinforced  concrete  for  ultimate  use  in  the 
cost-effective  design  and  hardness  assessmnet  of  concrete  protective  structun 
The  specific  goals  of  the  research  are:  (1)  development  of  a  mixture  theory 
which  can  accurately  account  for  steel-concrete  interaction  and  (2)  formulatii 
of  an  advanced  constitutive  theory  for  plain  concrete  which  can  accurately 
portray  its  nonlinear,  inelastic  behavior  including  damage  and  microcracking 
for  arbitrary  loading  histories.  The  importance  of  the  steel-concrete 
interaction  and  the  nonlinear  inelastic  behavior  of  the  palin  concrete, 
including  cracking,  is  emphasized,  • 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 
G3  UNCLASSIFIED/UNLIMITED  3  SAME  AS  RPT 


22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL 
Spencer  T  Wu 


DO  Form  1473,  JUN  86 


121.  ABSTRACT  SECURITY  CLASSIFICATION 

IJn class  if  ied 


^OFFICE, 


Pre wool  r-jiticnt  are  obsolete. 


SECURITY  CLASSIFICATION  OF  THIS  PA' 


Unclassified 


AD- A 187  337 


.*  rvM  ' 


i-.V 


iVV*  «  (?%V 


>4  Division  of  Maxwell  Laboratories,  Inc. 


afoak- ik-  a  7  -  l  3  a l 

SSS-R-87-8454 


DEVELOPMENT  OF  ADVANCED 


CONSTITUTIVE  MODELS 


FOR  PLAIN  AND 
REINFQliiiifcNCRETE 


*.* ' *V « ■.  *Vv ••  V *.V. ♦ ‘.V • Vv •  Vv •  •!:  • 

*.*•  ■•*,*  •  ••  ,»>•,*  •*••»*•«•••  •*•*■•**••*•  *  *  *  ••/*•*,*••* 


^  OCT  2  3  1987  ; 


«l|ipPiolling 

"^jpashington,  DC  2dlii?r:' 


March,  1987 


'  '  boa  appi-r 


;rf*  4  ' 


,•  ,  •  ua 

r  t  . j  -  ■'  ■ 


P.  O.  Box  1620,  Lo  Jolla ,  California  92038-1620 
(619)  453-0060 


VV#* 


*f  V 


■'*>♦  .|!‘V!y 'l»iv  ♦  i  *  i  ■*<*», 


La :  Li  ■ 1  j 


A  Division  of  MmxwoU  Laboratories,  Inc. 


SSS-R-87-8454 

DEVELOPMENT  OF  ADVANCED 
CONSTITUTIVE  MODELS 
FOR  PL^IN  AND 

REINForfBiONCRETE 


. .  \ V- Vv-’^V-V*  v*  *-*  \\** 


******  *’V*.  "*■*’**’ 


^^lashington,  DC  203§i 


March,  1987 


P.  O.  Box  1620,  La  Jolla ,  California  92038-1620 
(619)  453-0060 


ABSTRACT 


This  report  summarizes  research  conducted  by  S-CUBED  to  develop  advanced 
constitutive  models  of  plain  and  reinforced  concrete  for  ultimate  use  in  the  cost- 
effective  design  and  hardness  assessment  of  concrete  protective  structures.  The 
specific  goals  of  the  research  are  two-fold;  namely.  (1)  development  of  a  mixture 
theory  which  can  accurately  account  for  steel-concrete  interaction  and  (2)  formulation 
of  an  advanced  constitutive  theory  for  plain  concrete  which  can  accurately  portray  its 
nonlinear,  inelastic  behavior,  including  developing  damage  and  macrocracking,  for 
arbitrary  loading  histories.  The  importance  of  the  steel-concrete  interaction  and  the 
nonlinear  inelastic  behavior  of  the  plain  concrete,  including  cracking,  is  emphasized. 
The  progress  made  toward  achieving  these  goals  is  described. 
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Section  1 


INTRODUCTION 


1.1  BACKGROUND. 

The  Air  Force  has.  for  some  time  now,  had  a  major  interest  in  reinforced 
concrete  protective  structures,  which  ranges  from  missile  silos  --  to  buried  facilities  for 
personnel  and  equipment  —  to  above  ground  shelters  for  tactical  aircraft.  For 
example,  during  the  past  several  years,  the  Air  Force  has  been  deeply  involved  in  an 
extensive  effort  to  develop  and  assess  the  feasibility  and  relative  effectiveness  of 
various  candidate  modes  for  basing  the  MX  strategic  weapons  system.  In  most  of 
these  basing  modes,  the  key  elements  are  large  reinforced  concrete  structures 
(protective  structures)  which  are  designed  to  protect  a  missile  from  the  shock  loads 
prescribed  by  the  design  attack  scenarios.  The  enormous  costs  involved  in 
constructing  the  large  number  of  such  structures  required  by  the  system  dictates  that 
their  design  be  not  only  safe  but  cost-effective  as  well. 

The  current  methods  for  designing  protective  structures  are  largely  empirical,  with 
heavy  reliance  upon  engineering  experience  gained  from  the  design  of  structures  for 
prior  defense  systems,  and  insight  obtained  from  the  results  of  pertinent  Field  tests 
conducted  on  small-scale  structures.  Full-scale  testing  of  such  structures  under 
simulated  explosive  environments  is  extremely  rare  due  to  the  enormous  costs 
involved.  Analytical  simulations  with  finite  element  codes  are  often  performed,  but 
because  of  the  crudeness  of  the  material  models  for  the  plain  and  reinforced  concrete 
used,  little  confidence  is  placed  in  the  results  so  obtained,  especially  when  nonlinear 
effects  dominate  the  behavior.  At  best,  relatively  simple  elasto-plastic  models  of  plain 
concrete,  normally  developed  from  quasi-static,  low  pressure  laboratory  data,  are 
combined  with  simple  overlay  approaches  to  account  for  the  steel  reinforcement.  The 
resulting  models  of  reinforced  concrete  are  unable  to  properly  account  for  the 
interaction  that  occurs  between  plain  concrete  and  rebar  during  deformation.  Inasmuch 
as  such  models  are  generally  developed  on  the  basis  of  quasi-static,  low  pressure 
data,  they  are  not  capable  of  describing  the  response  of  concrete  to  the  high  pressure, 
impulsive  loading  environments  of  interest.  It  is  therefore  not  surprising  that,  at 
present,  the  collapse  (failure)  load  of  a  missile  silo  cannot  be  analytically  predicted  to 
within  a  factor  of  2.  Clearly,  there  is  a  crucial  need  to  understand  how  complex 
reinforced  concrete  structures  fail  under  complex  loads  applied  at  high  loading  rates. 
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The  inability  of  the  current  constitutive  models  to  adequately  describe  the 
nonlinear  behavior  of  plain  and  reinforced  concrete  is  the  main  reason  for  the  lack  of 
confidence  in  numerical  simulations  of  structural  response.  Without  a  valid  numerical 
simulation  capability,  there  is  little  hope  of  arriving  at  a  cost-effective  design.  The 
inaccuracies  associated  with  the  current  constitutive  models  of  reinforced  concrete  stem 
from  two  main  sources,  namely.  (1)  inadequate  models  of  plain  concrete  behavior, 
especially  for  high  pressures  and  high  rates  of  loading  and  (2)  overly  simplified 
methods  used  in  mixing  plain  concrete  and  steel  rebar  to  obtain  a  model  of  reinforced 
concrete.  In  both  cases,  the  nature  of  the  problem  requires  research  at  a  fundamental 
level,  since  the  basic  technology  is  missing.  Although  extensive  research  has  been 
conducted  over  the  past  fifty  years  on  reinforced  concrete,  the  emphasis  has  not  been 
on  the  high  pressure,  dynamic  response  regime.  Therefore,  if  methods  are  to  be 
developed  which  can  lead  to  cost-effective,  reliable  designs  for  concrete  defense 
structures,  improved  modeling  of  the  response  of  plain  and  reinforced  concrete  to  high 
pressure  loading  environments  is  needed,  as  well  as  corresponding  experiments  to 
validate  the  modeling.  This  requires  a  significant  advancement  in  current  technology. 

1.2  IMPORTANCE  OF  STEEL-CONCRETE  INTERACTIONS. 

Over  the  response  regime  of  interest  to  the  defense  community,  the  behavior  of 
reinforced  concrete  is  dominated  by  complex  nonlinear  interactions  between  steel  and 
concrete.  In  particular,  it  is  well  known  that  such  interactions  can  have  a  major 
effect  on  global  structural  characteristics  such  as  stiffness,  strength,  damping  and 
ductility.  In  addition,  steel-concrete  interactions  can  have  a  primary  influence  on 
damage  accumulation  and  failure  modes. 

To  illustrate  the  above  statement,  consider  the  example  of  reinforced  concrete  in 
direct  tension.  Practical  steel  percentages  (in  the  direction  of  loading)  range  from  3 
percent  ("dense"  layout)  to  0.5  percent  ("sparse”  layout).  If  one  examines  the 
unload-reload  tangent  modulus  associated  with  the  global  stress  and  strain  in  the 
direction  of  loading,  then  one  finds  that  the  initial  (uncracked  concrete)  modulus  drops 
by  approximately  65  percent  for  3  percent  steel,  as  damage  in  the  form  of  progressive 
cracking  and  steel-concrete  bond  degradation,  accumulates  with  increasing  global 
deformation  (extension).  For  0.5  percent  steel,  this  number  becomes  approximately 
95  percent!  For  the  global  deformation  domain  of  practical  relevance,  the  initial 
modulus  is  too  high  for  computational  purposes  and  the  asymptotic  (steel  only) 


modulus  is  too  low.  The  required  modulus  lies  in  the  transitional  region  between  the 
two  extremes.  However,  behavior  in  this  transitional  domain  can  be  predicted  only  if 
steel-concrete  interactions  are  considered.  Fc-  the  example  cited,  the  problem  of 
predicting  the  appropriate  stiffness  is  a  classical  problem  known  as  "tension 
stiffening" . 

Thus,  it  is  evident  that  a  model  of  reinforced  concrete  must,  in  some  sense, 
reflect  steel-concrete  interactions.  In  the  past,  practitioners  have  attempted  to  account 
for  some  of  these  interactions  in  an  effort  to  fit  test  data  by  "adjusting”  the 
constitutive  relations  of  either  the  steel  or  the  concrete  (usually  the  latter).  However, 
such  an  empirical  approach  requires  a  myriad  of  costly  experiments  —  costly  since  the 
specimen  size  must  be  "large"  when  compared  to  the  typical  steel  spacing  (material 
" microstructural"  dimension),  if  a  continuum  description  of  reinforced  concrete  is  to  be 
used.  In  the  final  analysis,  the  empirical  method  does  not  furnish  the  simulation 
accuracy  and  capability  that  are  needed  for  the  proper  treatment  of  current  defense 
problems  concerning  the  design  and  analysis  of  reinforced  concrete  protective 
structures. 

1.3  OBJECTIVES. 

The  primary  objective  of  the  research  described  here  is  to  construct  an  advanced, 
nonlinear,  multiaxia!  constitutive  model  of  reinforced  concrete  that  will  provide 
simulation  accuracy  in  the  nonlinear  response  regime  that  is  superior  to  existing 
models.  A  secondary  objective  is  to  develop  a  rational  method,  based  on  physics,  for 
accurately  accounting  for  the  effect  of  steel-concrete  interactions  on  the  overall 
response  of  reinforced  concrete.  A  final  objective  is  to  develop  an  accurate 
constitutive  model  of  plain  concrete  which  is  capable  of  describing  its  behavior  over  a 
wide  range,  from  the  onset  of  microcracking  to  the  development  of  complete  failure 
through  macrocracking.  These  objectives  are  consistent  with  the 
technological/operational  needs  exemplified  under  subsection  1.1  above. 

1.4  GENERAL  METHODOLOGY. 

An  innovative  methodology  has  been  selected  to  achieve  the  primary  objective  of 
the  research  program.  This  methodology  consists  of  synthesizing  the  global  properties 
of  reinforced  concrete  from  knowledge  of  the  plain  concrete  and  steel  properties,  the 
concrete-steel  interface  properties,  and  the  geometry  of  the  steel  reinforcement.  Such 


a  non-phenomenological  procedure  offers  the  potential  of  maximum  model  simulation 
accuracy/capability  via  inclusion  of  the  important  physical  phenomena,  and  minimizes 
the  number  and  type  of  expensive  tests  necessary  to  define  the  model  parameters. 

A  non-phenomenological  reinforced  concrete  theory  requires,  as  input,  constitutive 
models  for  each  of  the  basic  constituents,  namely,  steel  rebar  and  plain  concrete. 
The  accuracy  of  the  resulting  reinforced  concrete  theory  depends  heavily  on  the 
accuracy  with  which  one  can  model  these  two  constituents.  The  constitutive 
properties  of  steel  are  well  known  and  can  be  adequately  represented  by  reasonably 
simple  elasto-plastic  models.  Plain  concrete,  on  the  other  hand,  is  perhaps  the  most 
complex  structural  material  in  current  use  and.  despite  numerous  efforts  in  recent 
years  to  mathematically  model  its  properties,  there  is  today  still  no  model  which  can 
adequately  describe  its  nonlinear  constitutive  properties  over  a  wide  range  of  behavior, 
including  damage  accumulation  and  post-cracking  response.  In  addition,  a  non- 
phenomenological  reinforced  concrete  theory  requires  constitutive  descriptions  of  the 
steel-concrete  interfaces,  and  a  procedure  for  analytically  mixing  or  homogenizing  steel, 
concrete,  and  steel-concrete  interfaces  to  furnish  a  global  model  of  the  composite. 

As  a  result,  the  major  challenges  that  one  faces  in  attempting  to  develop  an 
accurate  constitutive  theory  for  reinforced  concrete  are  twofold,  namely 
(1)  development  of  an  analytic  mixing  procedure  (mixture  theory)  which  accurately 
accounts  for  the  interaction  between  the  rebar  and  plain  concrete,  and  (2)  formulation 
of  a  constitutive  theory  for  plain  concrete  which  adequately  describes  its  nonlinear, 
inelastic  properties,  including  damage  and  cracking,  for  arbitrary,  multi-axial  load  paths. 
Accordingly,  these  two  major  areas  of  research  have  been  pursued  in  parallel, 
concurrent  efforts  under  the  present  program. 

1.5  PRESENT  REPORT. 

This  final  report  describes  the  progress  made  during  the  contract  period  toward 
achieving  the  goals  described  above.  The  report  has  been  prepared  with  the  intention 
of  providing  the  reader  with  a  general  overview  of  the  research  accomplishments, 
without  becoming  cumbersome  with  extensive  detail.  Those  interested  in  pursuing  the 
various  topics  in  greater  depth  are  referred  to  the  publications  that  have  resulted  from 
this  study  (see  Section  2)  and  to  the  appendices  of  this  report  where  a  number  of 
topics  that  have  not  appeared  in  publication  are  discussed  in  detail. 
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Section  3 


REINFORCED  CONCRETE 

This  section  presents  the  results  of  research  conducted  on  the  theoretical 
description  of  reinforced  concrete.  The  discussion  covers  objectives,  general 
methodology,  mathematical  forms,  and  validation  problems  which  are  intended  to 
demonstrate  the  capability  and  accuracy  of  the  modeling  procedure. 

3.1  BASIC  OBJECTIVES. 

The  nonlinear  global  response  of  reinforced  concrete  protective  structures  to 
intense  dynamic  loads  in  dominated  by  two  fundamental  phenomena:  (1)  the  behavior 
of  the  plain  concrete  under  high  strain  rates  and  complex  multiaxial  stress  states;  and 
(2)  interactions  between  the  reinforcing  steel  and  the  concrete.  This  section  is 
concerned  with  steel-concrete  interaction. 

It  is  generally  recognized  that  steel-concrete  interaction  (SCI)  can  have  a 
dramatic  influence  on  the  nonlinear  response  of  a  reinforced  concrete  (R/C)  structure. 
In  particular,  important  response  characteristics  such  as  global  stiffness,  ductility, 
damage  type  and  evolution,  and  failure  modes/conditions  are  strong  functions  of  SCI. 

In  view  of  the  importance  of  SCI  and  the  current  lack  of  adequate  SCI 
descriptions  in  operational  R/C  computer  codes,  the  present  research  program  was 
undertaken  to  address  this  need.  One  of  the  basic  objectives  of  this  project  was  to 
construct  an  advanced  model  of  R/C  that  properly  reflects  SCI. 

3.2  ACCOMPLISHMENTS. 

A  fundamentally  new  mathematical  framework  for  describing  the  nonlinear 
response  of  reinforced  concrete  has  been  constructed  and  validated  against  a  suite  of 
experimental  data.  This  framework  is  nonphenomenological  on  the  macroscale  in  the 
sense  that  the  global  properties  of  reinforced  concrete  are  synthesized  from  the 
properties  and  geometries  of  the  components  (steel  and  concrete)  and  component 
interfaces.  Such  a  model,  while  far  more  difficult  to  construct  than  a 
phenomenological  model,  offers  the  possibility  of  increased  simulation  capability  and 
accuracy  with  a  minimum  number  of  experiments  necessary  to  evaluate  material 
parameters. 

The  foregoing  modeling  effort  has  proceeded  in  a  sequence  of  stages  that 
represents  increasing  levels  of  generality.  While  this  sequence  is  not  yet  complete,  the 


success  of  the  validations  performed  to-date  point  to  a  significant  breakthrough  in  the 
theoretical  description  of  SCI. 

3.3  METHODOLOGY. 

As  was  noted  previously,  the  nonphenomenological  approach  implies  that  a  model 
is  to  be  constructed  by  synthesizing  the  global  properties  of  R/C  from  the  properties 
and  geometries  of  the  basic  components  of  the  "composite",  i.e..  steel,  concrete, 
interfaces,  and  cracks.  Mathematically,  this  problem  falls  into  the  general  category  of 
"homogenization."  The  term,  as  used  here,  implies  construction  of  a  "macro"  model 
via  "micromechanical"  considerations  together  with  appropriate  smoothing  or  averaging 
operations. 

As  use  of  the  prefixes  "micro"  and  "macro”  above  imply,  homogenization 
typically  enters  the  picture  when  one  wishes  to  determine  the  response  of  a 
continuum  for  which  two  widely  different  length  scales  can  be  identified.  The  large  or 
"macroscale"  is  determined  by  the  specimen  geometry  and/or  the  loading  condition. 
The  small  or  "microscale"  is  determined  by  material  heterogeneity.  In  the  case  of 
reinforced  concrete,  the  typical  steel  spacing  constitutes  the  appropriate  microscale. 

The  homogenization  problem  is  parameterized  by  the  small  ratio  of  the  above 
two  length  scales,  e.  The  fundamental  problem  is  to  determine  the  "proper” 
macroscopic  response  equations  as  6  *  0.  Once  obtained,  it  is  natural  to  introduce 
asymptotic  notions  into  the  analysis  to  determine  a  physically  meaningful  sequence  of 
equations  ordered  in  powers  of  e.  The  higher  order  equations  are  intended  to  provide 
additional  simulation  capability  on  the  macroscale. 

Using  homogenization  concepts,  nonlinear  R/C  model  frameworks  have  been 
constructed  for  both  "dense"  and  "sparse"  unidirectional  steel  layouts.  Typical  such 
layouts  are  depicted  in  Figures  3.1  (a. b).  The  construction  technique  is  based  on  the 
use  of  multivariable  (two-space)  asymptotic  expansions,  a  variational  principle,  and 
certain  smoothing  operations.  The  resulting  models  have  been  cast  in  the  form  of  a 
binary  mixture  which  resembles  an  overlay  of  two  continua:  steel  and  concrete.  The 
continua  interact  via  body  forces  which  are  functionals  of  the  relative  global 
displacements  of  the  continua.  Mathematical  descriptions  of  the  model  and  the 
construction  procedure  can  be  found  in  Murakami  and  Hegemier  (1986).  Hegemier. 
Murakami,  and  Hageman  (1985),  Hegemier  and  Murakami  (1985).  Comprehensive 


treatment  of  the  model  construction  problem  is  furnished  collectively  by  Appendices 
C  -  I. 

3.4  REMARKS  ON  MIXTURE  FRAMEWORK. 

It  was  noted  previously  that  the  new  R/C  model  has  been  cast  in  the  form  of  a 
"mixture."  Some  clarification  concerning  mixture  descriptions  are  presented  below. 

Both  phenomenological  and  nonphenomenological  continuum  models  of  R/C  have 
features  that  resemble  a  binary  mixture.  This  mixture  assumes  the  form  of  a  two- 
phase  medium  with  various  types  of  constraints. 

Currently  operational  computer  programs  for  the  analysis  of  R/C  structures  are 
typically  based  on  phenomenological  models  which  require,  as  a  constraint,  the 
deformation  gradient  or  strain  tensor  of  the  steel  and  concrete  to  be  the  same  at 
spatial  points  where  they  are  overlayed  or  “smeared."  In  what  follows,  such  a 
description  will  be  denoted  as  an  "elementary"  mixture. 

At  the  other  end  of  the  spectrum  are  mixtures  which  allow  relative  deformation 
of  the  components,  i.e.,  they  exhibit  distinct  strain  fields.  Further,  these  models 
incorporate  steel-concrete  interface  physics.  Such  descriptions  will  be  denoted 
"advanced"  mixtures. 


For  the  purpose  of  demonstrating  the  difference  between  advanced  and 
elementary  mixture  descriptions,  consider  the  case  of  a  unidirectionally  reinforced 
specimen  subject  to  global  uniaxial  tension.  Figure  3.2.  Within  the  context  of  a 
quasi-one-dimensional  discussion,  an  advanced  mixture  would  have  the  form: 

(a)  Equilibrium  equations. 


(1)  (1)  +  p  _  o 

n  °11,1  rl  ”  u 


(2)  (2)  .  p  _  0 
*11,1  ri  “  u 


(3-1) 


(b)  Component  constitutive  relations. 


(c)  Interface  constitutive  relation 


pi  ■  F([uil)  •  tux] 3  ui(2)  -  ui(1) 


(3-3) 


Here  a ^  denote  averaged  (over  a  typical  "cell"  containing  an  element  of  steel  and 

a  representative  concrete  cover)  axial  stress,  displacement:  E.  C.  are  tangent  moduli, 
n^  and  n^  are  volume  fractions;  superscripts  1  and  2  refer  to  steel  and  concrete, 
respectively;  and  P  represents  an  interaction  term  which  is  related  to  an  effective 
steel-concrete  interface  shear  stress  r*  =  a*n  according  to: 
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Pj  =  2r*/r 


(3-4) 


where  r^  is  the  steel  radius  and  the  superscript  denotes  "interface."  The 

notation  F  in  (3-3)  implies  a  functional.  Equations  (3-1)  -  (3-3)  resemble  an  overlay 
of  two  materials  (steel  and  concrete)  which  interact  via  a  "body  force"  P  .  These 
equations  involve  four  field  variables:  a. and  u|a^.  a  =  1.2. 


a  =  1.2. 


A  conventional  description  .  or  elementary  mixture,  can  be  obtained  from  the 
foregoing  relations  by  imposing  the  constraint:  llj1^  =  uj2^  =  U  Thus,  the 
displacements  and  hence  the  strains  in  each  material  are  identical.  The  governing 
equations  then  reduce  to: 


(a)  Equilibrium 


a  -  0  =  n^V1)  ♦  n(2)ff(2) 

U  '  a  "  °  ffll  n  ffll 


(b)  Constitutive 

au  =  [n(1)E(1)  ♦  n(2)E(2)]01  1  ,  E  5  E  ♦  C  . 


(3-5) 


The  field  now  involves  only  two  variables: 
be  recovered  from: 

5(1)  =  e(Du  ;<2>  =  i  <2>o 

an  c  ui,i  *  ffn  c  1,1 


and  Uj. 


(3-6) 


The  component  stresses  can 
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Elementary  and  advanced  mixtures  generally  exhibit  drastically  different  simulation 
features  and  capabilities.  For  example,  if  the  concrete  is  modeled  as  elastic-brittle 
fracture,  and  the  steel  as  elastic,  then  the  elementary  model  in  extension  predicts  a 
constant  tangent  modulus.  E(.  Figure  3-3(a).  or  at  best,  a  piecewise  constant  tangent 
modulus  (drop  occurs  when  tensile  failure  strain  of  concrete  is  reached).  Figure  3-3(b). 
Here  E^  =  a/e  where  e  =  Uj  . 

In  contrast,  the  advanced  mixture  model  is  capable  of  exhibiting  realistic  behavior 

as  depicted  in  Figure  3-3(c).  Corresponding  differences  in  the  global  stress-strain 

response  are  shown  schematically  in  Figure  3-4  for  elastic-perfectly  plastic  rebar.  Note 

that,  in  contrast  to  the  elementary  model,  the  advanced  mixture  simulates  both 

stiffness  degradation  and  hysteresis.  The  interaction  term  P  and  its  functional 

dependence  on  [UJ  is  responsible  for  both  phenomena.  [In  Figure  3-4(c).  e  is  now 

defined  as  A £./£.  where  £.  =  initial  steel  specimen  length  and  E  =  (n^a.j1^  + 
(2)  (2)w  ,  o  o  K  6  t  *  11 

In  addition  to  the  obvious  differences  in  global  response  features  between 

elementary  and  advanced  mixtures,  there  is  a  difference  of  a  fundamental  nature  which 

concerns  the  prediction  of  failure.  As  was  noted  previously,  the  elementary  mixture 

assumes  the  same  strains  in  the  steel  and  concrete  at  a  spatial  point.  The  advanced 

mixture,  on  the  other  hand,  allows  distinct  strain  fields.  In  the  former  case,  failure  is 

predicted  when  the  global  extensional  strain  e  reaches  the  rebar  failure  strain,  ej1^.  In 

the  latter  case  (i.e..  with  use  of  the  advanced  mixture)  the  global  and  local  strains 

can  differ  significantly.  Figure  3-5.  Here  one  finds  that  failure  depends  not  only  on 

the  value  of  the  global  strain,  but  also  on  the  rebar  geometry,  the  steel-concrete  bond 

characteristics,  and  the  rebar  hardening  properties.  An  example  of  this  situation  is 

given  by  Figure  3-6.  [In  Figure  3-6.  D  =  rebar  diameter.  H  =  crack  spacing, 

ejs^  5  e]1)  =  rebar  failure  strain,  erf5)  =  Ultimate  rebar  stress,  =  initial  rebar 

t  f  f  .  .y 

yield  stress.]  Generally,  the  elementary  mixture  will  predict  extensional  failure  strains 

considerably  in  excess  of  actual  values;  hence  ductility  predictions  will  be 

nonconservative.  On  the  other  hand,  an  advanced  mixture  can  simulate  failure 

conditions  with  reasonable  accuracy. 

3.5  THEORETICAL  STUDIES. 

The  primary  theoretical  development  involves  a  number  of  key  steps;  these 
include:  (1)  a  judicious  scaling  of  the  problem  variables:  (2)  introduction  of  multiple 
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scales  (coordinates)  which  include  one  set  of  “macrocoordinates”  and  another  set  of 
“microcoordinates."  and  representation  of  all  field  variables  by  a  two-space  asymptotic 
method:  (3)  introduction  of  a  local  periodicity  condition:  (4)  derivation  and  application 
of  a  variational  principle:  and  (5)  conduct  of  an  asymptotic  analysis  to  ascertain 
appropriate  test  (trial)  functions  for  the  variational  principle.  Complete  details 
concerning  the  theoretical  construction  process  are  furnished  in  Appendices  C  -  I. 

To-date,  several  mixture  formulations  have  been  derived  and  explored.  The 
various  models  differ  in  their  form  and  completeness.  The  form  differences  are  a 
consequence  of  the  particular  variational  principle  used  and  the  associated  trial 
functions  employed,  as  well  as  the  description  (i.e..  Eulerian.  Lagrangian.  up-dated 
Lagrangian.  etc.).  A  brief  description  of  these  efforts  follows. 


3.5.1  Description  of  Main  Theoretical  Efforts. 


(1)  Direct  Variational  Method  --  A  very  general  set  of  mixture  equations  has  been 
generated  using  a  "direct”  variational  principle  in  which  only  displacements  are  varied. 
This  work  is  described  in  detail  in  Appendix  C.  The  resulting  model  includes  large 
deformations  and  nonlinear  material  and  interface  behavior.  The  latter  are  assumed  in 
an  incrementally  linear  form  which  includes  a  large  class  of  elastoplastic  and 
elastoplastic-brittle  fracture  relations.  The  theoretical  construction  process  has  been 
performed  for  both  dense  and  sparse  steel  layouts.  These  layouts  are  assumed  to  be 
unidirectional.  They  may,  however,  be  curvilinear. 

(2)  Mixed  Variational  Method  —  The  use  of  a  direct  variational  scheme  in  which 
only  the  displacement  field  is  varied  is  theoretically  convenient  from  the  point  of  a 
model  construction  process,  but  has  several  inherent  drawbacks.  For  example,  the 
displacement  and  stress  microfields  do  not  exhibit  the  same  order  of  accuracy.  Also., 
the  governing  relations  for  the  interaction  terms  do  not  evolve  naturally  from  the 
development.  In  addition,  the  partial  differential  equations  that  result  from  a  direct 
method  tend  to  be  too  “stiff'  in  some  applications.  In  an  effort  to  improve  this 
situation,  the  use  of  a  “mixed"  variational  principle  was  explored  with  considerable 
success.  This  work  is  described  in  detail  in  Appendices  D  -  G.  In  Appendix  D.  the 
basic  mixed  variational  principle  is  established.  In  Appendix  E.  a  mixture  model  is 
established  for  a  dense,  unidirectional  steel  layout  for  3-D  deformations.  In  Appendix 
F.  a  mixture  model  is  developed  for  R/C  plates  with  unidirectional  steel  layouts.  This 
work  is  further  specialized  to  a  R/C  rod  in  Appendix  H:  the  latter  model  is  necessary 


to  properly  interpolate  certain  laboratory  test  data.  The  foregoing  models  include 
nonlinear  material  and  interface  behavior;  however  the  treatments  are  restricted  to 
small  global  or  macro  deformations. 

(3)  Oblique  Steel  Layouts  --  The  models  developed  under  (1)  and  (2)  above  are 
applicable  to  general  global  (macro)  coordinate  systems.  However,  considerable  care 
must  be  exercised  in  this  connection  since  (a)  there  are  two  coordinate  systems 
associated  with  the  analysis;  macro  and  micro,  and  (b)  the  microcoordinates  must  be 
aligned  with  the  steel  layout;  hence  the  microcoordinate  system  is  a  material  system. 
Specific  equations  applicable  to  general  global  coordinate  systems  are  developed  in 
Appendix  G. 

(4)  Large  Microrotations  --  Some  problems  require  the  inclusion  of  large  rotations  of 
the  steel  on  the  microcoordinate  level  even  though  the  global  deformations  may  be 
small.  An  example  of  such  a  problem  concerns  the  transfer  of  global  shear  stress 
across  a  crack  when  rebar  penetrates  the  crack  surface.  For  the  purpose  of  exploring 
problems  of  this  type,  a  special  nonlinear  theory  was  developed  for  which  small 
deformations  are  imposed  on  the  macroscale,  but  moderate  rotations  of  steel  elements 
are  allowed.  This  work  is  described  in  detail  in  Appendix  I. 

3.5.2  Theoretical  Approach. 

An  extremely  brief  sketch  of  the  theoretical  approach  is  provided  next.  The 
reader  is  referred  to  Appendix  C  for  complete  details. 

Consider  a  uniaxial  steel  layout  with  a  (initial)  local  periodic  array  occupying  a 
volumeVo  with  boundary  9Vq  in  the  initial  configuration  as  illustrated  in  Figure  3-7. 
In  the  current  configuration  the  corresponding  volume  and  surface  will  be  denoted  as 
V  and  0V.  respectively.  Let  the  conservation  and  constitutive  relations  for  each 
material  (steel  and  concrete)  be  denoted  collectively  by 

'l  “  (l00^))  =  0  (a  =  1,2)  ,  (3-8) 
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where  a  =  1.2  refer  to  steel,  concrete  respectively  (no  sum  on  a):  x.  t  denote  the 
position  vector,  time  respectively;  is  the  collection  of  all  field  variables;  and  '/'  is 
a  nonlinear  operator.  The  form  (3-8)  represents  a  system  of  equations. 


It  is  essential  to.  at  the  outset,  nondimensionalize  (scale)  the  basic  equations  (3- 
8).  For  this  purpose,  two  fundamental  scales  are  introduced;  A  and  A.  Here  A  is  a 
typical  macrosignal  wavelength,  and  A  represents  the  typical  steel  spacing.  With  the 
aid  of  these  quantities,  the  field  equations  are  scaled,  and  new  independent  variables 


are  introduced  according  to  x  =  j/A.  t  =  t/t^  where  tj  ,  is  a  typical  macrosignal 
travel  time.  In  terms  of  the  nondimensiona!  variables,  the'  field  equations  assume  the 


travel  time.  In  terms  of  the  nondimensiona!  variables,  the  field  equations  assume  the 
form 


*  0 


(3-9) 


The  quantity  e  =  A/A  in  (3-9)  is  the  ratio  of  micro-to-macro  dimensions  of  the 
problem. 


Next,  the  concepts  of  microcoordinates  and  two-space  asymptotic  representations 
of  all  field  variables  is  introduced.  This  proceeds  as  follows:  It  is  expected  that 
stress  and  deformation  will  vary  significantly  with  respect  to  two  basic  length  scales: 

(1)  a  "global"  or  "macro"  length  typical  of  the  body  size  or  loading  condition  and 

(2)  a  "micro"  length  typical  of  a  "cell"  planar  dimension  as  depicted  in  Figure  3.8. 
These  macro-  and  micro-dimensions  will  be  associated  with  the  variables  A.  A 
respectively.  Further,  it  is  expected  that  these  scales  will  differ  by  at  least  one  order 
of  magnitude  in  most  cases.  Thus,  e  <<  1.  This  suggests  the  use  of  multivariable 
asymptotic  techniques.  Hegemier.  et  at.  (1979).  This  approach  commences  by 
introducing  new  independent  variables  according  to 


2*  =  f  00a  ,  *0)  4  0  as  e  +  0  , 


(3-10) 


2  =  *0)a  »  *00  4  1  as  f  4  0. 


For  the  present  problem,  it  will  suffice  to  set 


*0)  =  e  .  *(e)  =  1  . 


(3-11) 
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Thus,  all  field  variables  are  now  assumed  to  be  functions  of  the  "microcoordinates" 
g*  as  well  as  the  "macrocoordinates"  $  =j.  i.e.. 


=  £(a,t;2*;e)  .  (3-12) 

Spatial  derivatives  of  a  function  f(j.t)  then  take  the  form 

n  -  ee  *  «*Ye  (3-w) 


where  7,  7*  denote  gradient  operators  with  respect  to  x  and  coordinates, 
respectively. 

Under  (3-13),  the  field  equations  (3-9)  take  the  form 
^a)  [E W -o  (3-14) 

Next,  local  periodicity  of  all  field  variables  is  assumed  with  respect  to  the 
microcoordinates  j*.  This  condition,  which  is  an  extrapolation  from  linear  analysis  yet 
is  expected  to  provide  good  correlation  with  exact  results,  allows  one  to  consider  a 
typical  cell.  Figure  3.8  for  the  purpose  of  determining  the  variation  of  with 
respect  to  5*.  This  cell  contains  a  representative  area  of  steel,  A^.  and  concrete. 

The  field  equations  (3-14)  are  now  satisfied  in  an  average  sense.  Weak 
solutions  of  (3-14)  are  obtained  via  a  weighted  residual  procedure  which  takes  the 
form: 


II  ''<<0 (f (<0 (a, t ;/;*)]  •  KCa) [is;a‘;^dA'  dv  =  o  (3-is) 

•  v"  1 


where  denote  weighting  functions,  and  •  denotes  an  appropriate  inner  product. 

The  residual  procedure  (3-15)  is  converted  to  a  variational  principle  by  selecting 
a  subset  of  the  weighting  functions  as  variations  of  velocity  (Direct  Method)  or  a 
combination  of  variations  in  velocity  and  stresses  (Mixed  Method).  The  Euler 
equations  of  the  variational  process  are  the  desired  mixture  relations.  The  x* 
dependence  of  the  weighting  functions  is  suggested  by  a  separate  asymptotic  analysis. 


3.5.3  Typical  Mixture  Form. 


In  the  case  of  the  Direct  Method,  the  velocity  field  of  each  constituent  is  given 
by 

v t;x* ;e]  =  (x,t)gK^  )  ♦  0(e2)  (3-16) 

The  dependence  on  $  reflects  variations  of  the  velocity  field  over  the  cell.  i.e..  the 
microstructure.  The  functions  gK^  (sum  over  K  from  1  to  n)  are  given.  The 
variables  U  and  §  are  displacement-type  dependent  variables. 

When  (3-16)  is  used  as  trial  functions,  the  basic  mixture  relations  take  the  form 
£  •  (j(a)  .  (-I)"*1  £  =  0  V,  («.  =  1,2)  (3-17) 

H(“)  ='4xy(a>  .  !jKta)  •  ?W)  •  («  =  1,2)  (3-18) 

E  =  *([6<1)  -  8(2)]  .  (a))  ■  («  -  1,2)  (3-19) 

In  the  above,  is  a  stress  average  over  A1®*.  £  is  an  interaction  term,  and  ^ 
are  functionals.  Equations  (3-17)  -  (3-19).  which  represent  a  static  problem,  have  the 
form  of  an  overlay  of  two  continua  which  interact  via  P.  Equation  (3-17)  is  the 
equilibrium  condition,  equation  (3-18)  is  the  constitutive  relation,  and  equation  (3-19) 
is  the  constitutive  relation  for  the  interaction  term. 

The  mixture  relations  (3-17)  -  (3-19)  account  for  the  behavior  of  concrete,  steel, 
and  steel-concrete  interaction.  The  latter  coupling  occurs  via  the  quantity  £.  Thus, 
for  example,  bond  slip  and  degradation  is  reflected  by  (3-19).  The  operators  Jfand 
.J(  are  defined  once  the  steel,  concrete,  and  steel-concrete  interface  behavior  is 
analytically  described. 

It  is  noted  that  cracks  in  the  concrete,  according  to  the  present  modeling 
scheme,  are  treated  explicitly.  This  furnishes  maximum  simulation  capability.  Another 
level  of  smoothing  (smoothing  of  steel  and  intact  concrete  results  in  the  mixture 


relations  (3-17)  -  (3-19))  is  possible  wherein  the  crack  field  is  smoothed.  The 
problem  with  such  an  operation,  however,  is  that  damage  progresses  from  a  sparse 
crack  field  to  a  dense  field  and  the  entire  range  of  behavior  may  be  important. 
Further,  many  failures  are  localized,  not  distributed. 

3.6  VALIDATION  STUDIES. 

The  purpose  of  this  section  is  to  explore  the  simulation  accuracy  and  capability 
of  the  developed  mixture  theory  framework,  and  to  illucidate  certain  important  features 
of  the  model.  To  aid  in  this  task,  several  basic  problem  types  have  been  selected  for 
detailed  studies.  These  include 

•  Steel-concrete  slip  and  pullout 

•  Behavior  under  direct  tension 

•  Behavior  under  biaxial  stress  states 

•  Behavior  in  shear  (the  dowel  problem) 

•  Behavior  in  direct  compression. 

3.6.1  Simulation  of  Steel- Concrete  Slip  and  Pullout. 

This  problem  type  plays  a  central  role  in  the  evaluation  of  certain  interaction 
terms  in  the  mixture  formulation.  Its  relevance  is  as  follows:  The  concept  of  steel- 
concrete  "bond”  and  "slip"  is  manifested  theoretically  via  Equation  (3-19)  or  Equation 
(14)  of  Appendix  C.  What  is  required  in  these  equations,  and  in  the  mixture  model, 
is  a  local  traction  versus  slip  relation,  e.g..  interface  shear  stress  versus  local  interface 
displacement  difference  between  the  rebar  and  concrete.  However,  such  a  relation 
cannot  be  obtained  directly  via  experimental  techniques.  Although  some  measures  of 
local  slip  can  be  obtained  using,  e.g..  slip-wire  techniques,  the  local  shear  stress 
cannot  be  measured.  One  must,  in  fact,  back  the  local  interface  behavior  out  of  a 
boundary  value  problem.  And  solution  of  a  boundary  value  problem  requires  a 
theoretical  model.  Thus,  a  model  must  first  be  postulated,  then  a  boundary  value 
problem  formulated  which  adequately  represents  the  bond-slip  test,  and  finally  the 
desired  result  backed  out  of  a  suite  of  simulations  of  the  test. 

To  illustrate  the  above  statements,  consider  the  test  schematic  of  Figure  3.9 
which  shows  a  typical  pullout  test.  Suppose  that  the  mixture  theory  developed  herein 
is  adopted  as  the  governing  model  of  the  experiment.  Then,  for  small  deformations 
and  elastic  component  response,  the  theory  developed  in  Appendix  G  reduces  to 
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Figure  3.9.  Typical  pullout  test  specimen. 
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In  the  above,  a  =  1.2  refer  to  steel,  concrete,  respectively:  u|a^.  denote  axial 

displacement  and  axial  stress  (average),  respectively;  K  is  a  tangent  modulus:  n^ 
denotes  the  volume  (area)  fraction  of  material  a:  an<JP()  =  9()/9xlf(*)  =  9()/9t 
where  t  represents  time. 

The  relations  (3-20)  are  the  equilibrium  equations  for  each  material:  the  latter 
interact  via  the  interaction  body  force  which  reflects  interface  shear  transfer 
between  the  steel  and  concrete.  If  a*  denotes  the  interface  shear  stress,  then 

rx 


pi  -  ^  ■ 


(3-23) 


Equations  (3-21)  are  the  constitutive  equations  for  each  material.  For  elastic 
component  response,  the  constants  are  given  by  Young  s  modulus. 

Prior  to  interfacial  slip,  the  mixture  theory  developed  herein  (e.g..  see  Appendix 
H)  predicts  that 


fK  ^T1  £?  1  [2  ♦  ♦  izin  n(1^/n(2)]l  _  K-1 

[  =  8  U>  '  n<2V2>  )  • 


(3-24) 


Thus,  prior  to  slip,  the  behavior  of  the  interaction  term  P^  is  linear  in  [Uj.  Recall 
that  Ufl  represents  the  average  axial  displacement  of  material  a:'  Thus,  even  for  the 
case  of  no  local  slip  at  the  steel-concrete  interface,  the  quantity  [UJ  does  not  vanish, 
i.e.,  llj1)  t  Uj2>.  The  relation  between  the  local  slip  and  [UJ  can  be  deduced  from 
microstructural  relations. 

When  steel-concrete  slip  occurs.  Equation  (3-24)  must  be  generalized.  For 
monotonic  global  deformation,  a  simple  elastic-perfectly  plastic  model.  Figure  3.10(b) 
will  often  suffice.  When  the  elastic  deformations  of  the  concrete  and  steel  are 
subtracted,  this  corresponds  to  a  local  rigid-plastic  description  of  the  interface  bond 
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stress  as  depicted  in  Figure  3.10(a).  Additional  accuracy  can  be  obtained  by  more 
complex  interface  descriptions  such  as  depicted  in  Figure  3.11. 

Let  us  turn  now  to  experimental  versus  theoretical  comparisons.  Using  the 
above  mixture  model,  simulations  of  a  number  of  pullout  tests  from  different 
laboratories  were  performed.  For  this  purpose,  a  boundary  value  problem  was  defined 
wherein  the  actual  specimen  geometry  and  properties  were  accounted  for.  Simulations 
were  conducted  by  numerical  solution  of  (3-20)  -  (3-22).  Typical  results  are 
illustrated  in  Figure  3.12.  The  elementary  elastic-perfectly  plastic  interaction  term 
(rigid-perfectly  plastic  interface  shear  stress  versus  local  slip)  is  observed  to  furnish 
adequate  accuracy  whereas  the  more  complex  multilinear  bond-slip  relation  provides 
excellent  agreement.  Another  comparison  between  simulation  and  experiment  is  given 
in  Figure  3.13.  Here  again,  excellent  accuracy  is  observed. 

Based  upon  a  collection  of  simulation  results,  the  piece-wise  linear  description  of 
local  bond  stress  versus  local  slip  furnished  in  Figure  3.11  was  found  to  provide  an 
adequate  description  of  interface  behavior  for  a  variety  of  concrete  strengths  and 
concrete  covers.  Using  this  description,  simulation  accuracy  was  found  to  vary  from 
the  excellent  results  of  Figure  3.13  to  the  adequate  results  of  Figure  3.14. 

For  hysteretic  deformation,  the  above  description  of  the  interaction  term  must  be 
generalized  considerably.  Several  bond-slip  models  that  have  been  explored  for  this 
purpose  are  depicted  in  Figure  3.15. 

3.6.2  Simulation  of  R/C  Under  Uniaxial  Tension. 

One  important  form  of  damage  in  R/C  occurs  because  the  maximum  global 
(composite)  stress  is  positive,  i.e..  tension.  Such  stress  states  can  result  in 
progressive  crack  accumulation  and  associated  steel-concrete  interface  slip.  One 
manifestation  of  such  damage  is  stiffness  degradation.  This  phenomenon  plays  a 
dominant  role  in  the  bending  and  nonlinear  stretching  of  R/C  beams,  plates,  and 
shells. 

To  understand  how  well  the  mixture  formulation  simulates  the  nonlinear  response 
of  R/C  when  the  maximum  principal  stress  is  positive,  the  classical  tension-stiffening 
problem  was  selected  for  simulation.  Details  of  the  resulting  studies  can  be  found  in 


sal  bond  slip  relation  for  specimens  axially  reinforced  with  deformed  bars  as  a  function 
concrete  strength  and  cover. 


Simulation  versus  experiment  for  pullout  tests. 
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Figure  3.13.  Simulation  of  pullout  test  (average  bond  stress,  r,  is  the 
pullout  force  divided  by  the  steel  surface  area). 


Hegemier.  et  a/..  (1985),  Murakami  and  Hegemier  (1986).  and  Hageman.  et  at.. 
(1986).  In  what  follows  the  basic  features  of  the  problem  and  the  results  are 
sketched. 

Consider  the  problem  of  simulating  the  global  response  of  a  unidirectionally 
reinforced  concrete  specimen  subject  to  a  state  of  global  uniaxial  stress  applied  in  the 
steel  direction  (principal  stresses  that  are  oblique  to  the  steel  will  be  considered  later) 
as  depicted  in  Figure  3.2.  If.  as  in  the  previous  section,  small  deformations  are 
assumed,  then  the  theory  developed  in  Appendix  E  reduces  to 


|(*P)  +  p  =0 
11,1  rl  u  ' 

N(2P)  _  p  _  0 
"ll,l  1  U  ' 

(3-25) 

i(*P)  _  „(«)&(«■) 
11  ~  n  "ll 

(3-26) 

1  ■  K.p[°l]  ■  [°: 

J  5  K(1)  -  5{2)1 

(3-27) 

Equations  (3-25)  -  (3-27)  are  similar  in  form  to  (3-20)  -  (3-22).  They  differ, 
however,  in  the  definition  of  the  tangent  moduli.  Prior  to  bond  slip,  and  for  the 
elastic  response  regime,  the  tangent  moduli  in  (3-26)  and  (3-27)  are  given  by 

E  ^  =  (X  .  2/i)  -  X(a)2/(X  .  ,  (3-28) 


fj_  [2  «  J2>  .  (2£n  nH)/nW)ll 

'  n<2V2)  1 


(3-29) 


where  ft' a ^  are  the  Lame  constants  of  material  a. 

When  steel-concrete  slip  occurs.  Equation  (3-29)  must  be  generalized.  For 
monotonic  global  deformation,  a  simple  elastic-perfectly  plastic  model  is  sufficient.  For 
hysteretic  deformation  detailed  studies  of  steel-concrete  pull-out  specimens  have  guided 
the  construction  of  an  elementary  bond-slip  model  which  is  depicted  in  terms  of  the 
interaction  term  in  Figure  3.16.  For  a  range  of  concretes,  the  following  selection 
of  parameters  is  appropriate: 
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Pcr  =  2^)(JrJcr/£  ,  PF  =  -  P0  =  0.8  P, 


(3-30) 


0.2  P 


cr 


0.1  P 

cr 


The  model  depicted  in  Figure  3.16  is  easily  placed  in  analytical  incremental  form. 
The  description  involves  two  parameters.  K  as  defined  by  (3-29)  and  (ofx)CT  which 
represents  a  critical  steel-concrete  interface  shear  stress. 

Outside  of  the  elastic  response  regime,  the  components  will  be  modeled  as 
follows:  (1)  concrete:  elastic-brittle  fracture:  (2)  steehelastic-plastic  with  strain 
hardening.  This  behavior  is  depicted  in  terms  of  tangent  moduli  in  Figure  3.17. 

(1)  Monotonic  Extension. 

The  response  of  the  R/C  type  noted  above  involves  three  primary  nonlinear 
phenomena:  (1)  steel-concrete  bond  slip.  (2)  plastic  deformation  of  the  rebar,  and 
(3)  progressive  cracking  of  the  concrete.  Items  (1)  and  (2)  are  simulated  by  the 
basic  mixture  relations  once  the  appropriate  tangent  moduli  are  specified.  However, 
item  (3)  requires  comment.  In  the  current  mixture  formulation,  cracks  are  treated 
explicitly,  i.e..  they  are  not  smoothed.  This  furnishes  maximum  simulation  capability 
and  allows  treatment  of  both  sparse  and  dense  crack  fields.  However,  this  approach 
requires  one  to  follow  the  entire  crack  evolution  sequence.  For  the  present  uniaxial 
stress  problem,  only  primary  (through)  cracks  are  considered.  The  development  of  such 
cracks  are  based  on  a  simple  stress  criteria.  The  stress  measure  used  in  this  case  is 
the  average  axial  concrete  tensile  stress.  Njj2a^.  Between  two  cracks,  the  variation  of 
this  stress  is  as  depicted  in  Figure  3.18.  Thus,  new  cracks  always  form  halfway 
between  old  cracks  (neglecting  statistical  variations  in  tensile  strength).  For 
simulation  purposes,  the  boundary  value  problem  depicted  in  Figure  3.19  was  assumed. 
This  Figure  also  shows  the  resulting  crack  sequence. 

Consider  now  monotonic  extension.  Figure  3.20  shows  simulation  versus 
experimental  (Mirza  and  Houde  (1979))  data  for  two  sizes  of  tension  specimens.  The 
specimen  geometry  is  depicted  in  Figure  3.21.  Agreement  is  considered  to  be 
satisfactory  when  one  considers  the  data  scatter  that  is  typical  of  such  tests.  It  is 
emphasized  that  the  bond  slip  law  defined  by  Figure  3.11  was  used  in  the 
simulations.  It  is  noted  that  no  effort  was  made  to  "tune"  the  model  to  the  data. 
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Figure  3.18.  Stresses  on  steel  and  concrete 
and  at  interface. 


Another  comparison  that  was  made  involves  monotonic  extension  of  scaled  R/C 
specimens  with  a  dense  unidirectional  steel  layout.  Figure  3.22.  Experimental 
(Somoyaji  (1979))  versus  theoretical  comparisons  for  load  controlled  tests  are  shown 
in  Figures  3.23  —  3.25.  Figure  3.23  illustrates  the  behavior  of  one  particular  specimen 
whereas  Figures  3.24  -  3.25  depict  comparisons  for  a  suite  of  tests.  Because  of  the 
load-control,  the  tests  in  this  series  were  not  "clean."  In  particular,  expected  sharp 
drops  in  global  stress  due  to  progressive  fracture  of  the  concrete  are  "smeared"  due 
to  the  specimen  loading  procedure.  In  addition,  the  method  of  testing  led  to 
significant  data  scatter.  Nevertheless,  the  experimental-theoretical  comparisons  shown 
demonstrate  that  the  mixture  formulation  is  capable  of  good  simulations  for  a  wide 
range  of  steel  layouts  in  all  major  response  ranges.  Figure  3.26.  The  latter  consists 
of:  Elastic  response  with  no  cracking,  debonding,  or  slip  (Stage  I),  progressive 
cracking  with  debonding  and  slip  (Stage  II),  and  slip  only  (Stage  III).  In  addition  to 
global  stress-strain  response,  reasonable  agreement  has  been  observed  between 
predicted  and  observed  crack  patterns. 


In  addition  to  monotonic  extension,  comparisons  were  also  made  between 
simulations  and  experimental  data  for  hysteretic  extension.  The  case  considered 
involved  cyclic  loading  (tension-compression)  of  full-scale  plate-type  (64  in.  x  64  in.  x 
8  in.)  specimens  reinforced  with  No.  5  rebar.  In  contrast  to  the  scaled  R/C  specimen 
tests,  the  experiments  in  this  series  were  carefully  conducted  under  displacement 
control  (Hegemier.  Nunn  and  Arya  (1978)).  Figures  3.27  and  3.28  indicate  typical 
experimental  versus  theoretical  comparisons  of  global  (composite)  stress-strain 
response.  Figure  3.29  compares  simulated  and  measured  specimen  stiffness 
degradation  as  measured  by  average  unload-reload  slopes.  The  agreement  of  test 
results  and  theoretical  predictions  is  observed  to  be  excellent.  It  is  noted  that  both 
bond  slip  and  progressive  cracking  of  the  concrete  were  active  in  these  tests  and 
simulations. 

In  addition  to  the  above  global  stress-strain  envelop  and  stiffness  comparisons, 
detailed  comparisons  were  made  between  theory  and  experiment  for  each  unload-reload 
cycle.  A  sample  of  these  comparisons  is  furnished  in  Figure  3.30  It  is  evident  that 
the  mixture  model  formulation  is  capable  of  simulating  complex  details  concerning  the 
hysteretic  response  with  remarkable  accuracy. 
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Figure  3.25.  Comparisons  of  experiment  and  theory. 
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Figure  3.29. 


Comparison  of  test  and  simulation:  Tangent  stiffness 
degradation  due  to  cracking  and  debonding. 
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3.6.3  Simulation  of  R/C  Under  Shear  (Dowel  Problem 


The  dowel  problem  is  another  classical  problem  associated  with  the  response  of 
R/C  subsequent  to  cracking.  This  case,  in  its  most  elementary  form,  is  depicted  in 
Figure  3.31  for  a  dense  array,  and  in  Figure  3.32  for  a  "cell".  Here  a  unidirectionally 
reinforced  specimen  with  a  preexisting  crack  perpendicular  to  the  rebar  is  subjected  to 
global  shear.  Of  interest  is  the  prediction  of  the  global  shear  force  or  stress  across 
the  crack  plane  for  a  prescribed  relative  offset  displacement  of  the  concrete.  In 
general,  the  response  involves  nonlinear  geometric  deformation  (moderate  rotations)  of 
the  steel.  Figure  3.33..  transverse  deformation  of  the  steel  due  to  "dragging"  through 
the  concrete,  steel-concrete  bond  slip,  and  aggregate  interlock  between  opposing  crack 
faces. 

For  the  geometry  depicted  in  Figure  3  32.  the  nonlinear  mixture  relations  of 
Appendix  C  can  be  employed  to  deduce  the  appropriate  governing  equations. 
However,  the  description  of  Appendix  C  is  Eulerian  whereas  a  Lagrangian  description 
is  more  convenient  for  this  problem.  Rather  than  presenting  a  general  discussion  of 
the  nonlinear  mixture  relations  in  a  Lagrangian  description  and  then  specializing  to  the 
present  problem,  a  special  derivation  for  the  dowel  case  is  furnished  in  Appendix  I. 
The  basic  equations  are: 

(a)  Equilibrium 
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(b)  Stress  Constitutive: 
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Figure  3.32.  A  dowel  action  test  simulation. 
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where  for  f  <  0  or  f  =  0  and  ?  <  0. 
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where  E./j  denote  Young's  modulus  and  the  shear  modulus,  respectively,  and  for 
f  =  0  and  ?  =  0. 
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C00  =  n(1)/*(1>A  ,  C, 


(3-37) 


(c)  Interaction  Constitutive 


pi  f°r  ipii  >  (pi)»  • 


(i  =  1,2;  no  sum  on  i) 


(3-38a) 


(3-38b) 


where 


f  *  f  ;  1  -  F^h:  • 
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In  the  above.  N.:  p'  is  the  partial  stress  associated  with  the  steel.  A  represents 
the  cell"  area,  P.  is  the  interaction  body  force  in  the  \th  direction.  Mjj  p^  is  a 
weighted  stress  average.  E^  and  denote  Young's  modulus  and  the  shear  modulus 
of  the  steel.  H  and  H.  are  hardening  parameters,  and  Sj  is  a  displacement-type 
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variable.  An  assumption  of  moderate  rotations  of  the  steel  is  implicit  in  these 
relations. 


(d)  Boundary  Conditions: 

For  the  case  of  a  smooth,  lubricated  crack  interface,  the  appropriate  boundary 
conditions  are 

=  0,  N^P)  =  0,  sOatXj.a.  (3-39) 

(e)  Forcing  Functions 

The  "forcing  functions"  of  the  problem  are 

py  py  ■  wd)l-«2a)  <3-*o> 


where  the  offset  A  is  prescribed. 


The  relations  (3-39)  above  exclude  interface  shear  transfer  (1ST).  An  elementary 
procedure  for  inclusion  of  1ST  is  as  follows:  First.  ( 3- 3 1 ) - ( 3-40)  are  solved 
incrementally  for  u!1^  and  T|la^  at  =  0  where  t|13^  =  N j*a).  T^Ia'  =  + 

U2W  Njj1'.  Second,  an  effective  normal  stress  p  =  po  +  (V|la^/A)  is  defined  on 
the  crack  face  at  x^  =  0  where  pQ  is  a  confining  pressure  (see  Figure  3.33)  which  is 
part  of  the  test  set-up  for  combined  dowel  action  (DA)  and  1ST  experiments.  Third, 
the  effective  normal  stress  p  is  applied  to  the  concrete  interface  at  x}  =  0  and  an 
1ST  law  is  postulated  in  the  form  f(T^2a^.  p.  A)  =  0  and  this  relation  is  solved  for 
T^2a^  for  a  specified  p  and  A.  This  latter  quantity  reflects  1ST.  According  to  this 
procedure.  1ST  influences  the  total  effective  shear  stress,  r  across  Xj  =  0  but  does 
not  alter  the  DA  contributions.  T^la^.  Here 


T  = 


|(la)  +  j(2a) 
12  +  2 


(3-41) 


(1)  Monotonic  Deformation. 

Let  us  now  examine  the  simulation  capability  of  the  foregoing  model  for 
monotonic  deformation.  The  first  case  to  be  considered  is  a  simulation  of  the 
lubricated  monotonic  dowel  tests  by  Paulay,  Park  and  Phillips  (1974).  The  test  set¬ 
up  and  specimen  are  illustrated  in  Figures  3.34  and  3.35.  Wax  was  used  to  lubricate 
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the  joint  surface.  Comparisons  of  theoretical  and  experimental  results  are  shown  in 
Figure  3.36  for  three  different  steel  volume  fractions.  Agreement  is  observed  to  be 
good.  The  "shear  stress"  in  Figure  3.36  is  based  on  (3-41)  with  lj2a^  =  0. 
Simulation  of  the  test,  which  approximates  pure  DA.  was  performed  using  (3-31  )-(3- 
40)  and  an  appropriate  incremental  numerical  procedure.  The  model  parameters 
employed  in  the  simulation  are  given  in  Table  3-1.  The  forms  of  the  interaction 
terms  used  are  depicted  in  Figures  3.37  and  3.38. 


The  second  case  to  be  considered  includes  1ST:  the  tests  were  performed  by 
Paulay.  Park  and  Phillips  (1974)  using  the  test  set-up  and  the  test  specimens  (of 
different  joint  surface  preparation)  depicted  in  Figures  3.34  and  3.35.  Figure  3.39 
shows  comparisons  of  experimental  data  and  simulations  for  several  different  interface 
conditions.  For  simulation  purposes,  an  1ST  relation  of  the  form 
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was  adopted.  The  model  parameters  used  in  the  simulation  are  given  in  Table  3-2. 
Agreement  between  simulations  and  test  data  is  observed  to  be  excellent. 

Another  case  of  interest  concerns  the  data  produced  by  Karagozian  and  Case 

(1973).  In  these  experiments  the  maximum  slip  A  was  much  larger  than  that 

associated  with  the  aforementioned  tests.  The  test  set-up  and  the  specimens  used 

are  depicted  in  Figures  3.40  and  3.41.  In  order  to  include  the  effect  of  confining 

pressure  p  on  the  1ST  relation,  r  and  r  are  expressed  as: 

max  * 

r  =  979  +  1.53  p  (psi) 

max  'r  ' 

r00  =  76  +  0.84  p  (psi). 

The  remaining  model  parameters  employed  for  simulation  purposes  are  given  in  Table 
3-3.  comparisons  between  simulations  and  test  results  are  illustrated  in  Figure  3  42 
for  a  range  of  confining  pressures.  Agreement  between  predictions  and  experimental 
data  is  observed  to  be  reasonable. 
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Figure  3.36..  Comparison  of  Paulay  load-slip  data  with  simulations 


800 


m 


m 


■A#* 


Scobbled 


Keyed 


600 

a. 


CO 

co 

Id 

CC 

w  400 


retarder 
^no  bond 


Jrowelled  plus  xxx 


trowelled 


0.02 


0.04  0.06 

SLIP  (in) 


a  SBOI 
°  KBOI 
+  ZBOI 
°  RBOIX 
*  TBOI 

0.08 


Figure  3.39. 


Comparison  of  simulations  and  load-slip  data  for 
different  surface  conditions. 
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Figure  3.40.  Test  set  up  used  by  Karagozian  and  Case  (1973). 
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Figure  3.41. 


Test  specimen  used  by  Karagozian  and  Case  (1973). 
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(3)  Hysteretic  Deformation 


A  critical  test  of  the  model  simulation  capability  is  represented  by  hysteretic 
deformation.  For  validation  purposes  the  experiments  of  Jimenez.  Gergely  and  White 
(1978)  were  selected.  The  test  set-up  and  specimen  are  illustrated  in  Figures  3.43 
and  3.44.  Thin  brass  sheets  were  used  to  lubricate  the  joint  surface.  A  typical 
simulation  versus  experiment  is  shown  in  Figure  3.45.  The  agreement  is  observed  to 
be  good  when  one  considers  the  complexity  of  the  response.  The  model  parameters 
employed  in  the  simulation  are  provided  in  Table  3-4. 


TABLE  3-4 

PARAMETERS  USED  IN  SIMULATION  OF  TESTS 
Jimenez,  Gergely,  and  White 
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Figure  3.43. 


Test  set  up  used  by  Jimenez,  et  a!.,  (1978). 
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Figure  3.44. 


Tests  :  C4 -7A,  C4-9A  and  C4-9B 
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Test  specimen  used  by  Jimenez,  et  a!.,  (1978). 


3.7  STUDIES  IN  PROGRESS. 


A  number  of  studies  remained  in  progress  as  of  the  writing  of  this  report.  Two 
of  the  more  important  efforts  are  described  below. 

3.7.1  Biaxial  Stress  States. 

The  validation  studies  described  previously  were  restricted  with  respect  to  the 
stress  states  (global)  that  were  considered.  For  example,  in  the  discussion  of  R/C 
under  tension,  the  global  stress  state  considered  was  uniaxial.  Most  practical 
problems,  however,  involve  either  global  biaxial  stress  states  (R/C  plate)  or  global 
triaxial  stress  states  (fully  3D  R/C  layout). 

The  models  developed  thus  far  are  capable  of  simulating  the  response  of  R/C 
under  either  biaxial  or  triaxial  conditions.  In  an  effort  to  validate  the  theory  for  such 
cases,  the  problems  depicted  in  Figures  3.46,  3.47  were  considered.  In  all  cases,  the 
rebar  layout  is  oblique  to  the  global  principal  stress  directions.  In  Figure  3.46(a)  the 
objective  was  to  understand  how  the  rebar  angle  l  effects  the  global  stiffness 
degradation  as  the  crack  field  evolves.  Figure  3.46(b)  includes  the  effect  of  a  second 
(compressive)  principal  global  stress  as  well  as  The  problem  depicted  in  Figure 
3.47  is  representative  of  a  R/C  shear  wall  with  two  extreme  cases  of  boundary 
conditions. 

In  each  of  the  foregoing  cases,  the  intent  was  to  compare  simulations  with 
laboratory  data.  For  this  reason,  relevant  test  data  on  both  model  R/C  and  full-scale 
R/C  tests  have  been  collected.  These  data  are  currently  under  study. 

3.7.2  Description  of  Concrete  in  Compression . 

This  subject  matter  concerns  the  behavior  of  reinforced  concrete  in  compression. 
To  model  this  behavior,  it  is  necessary  to  adequately  describe  the  plain  concrete 
component  in  the  mixture  theory.  This  task  was  addressed  in  considerable  detail.  A 
complete  presentation  of  the  work  is  provided  in  Appendix  J. 

It  is  well  known  that  the  ductility  of  concrete  increases  as  the  mean  compressive 
stress  increases.  The  current  approach  followed  by  most  investigators  is  to  model 
this  effect  via  a  conventional  elastic-plastic  theory.  However,  substantial  increases  in 
global  ductility  of  concrete  are  observed  well  below  the  brittle-ductile  transition. 
Thus,  it  is  unlikely  that  plastic  flow  is  responsible  for  this  effect.  A  more  plausible 
explanation  concerns  "sliding"  across  distributed  microcracks. 
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Recent  analytical  micromechanical  studies  appear  to  support  the  following  view  of 
the  compression  problem:  under  low  confining  stress  microcracks  tend  to  combine  to 
form  a  single  macrocrack  or  fault  at  failure:  as  the  confinement  is  increased,  the 
microcracks  tend  to  form  a  distribution  which  resembles  a  slip  system.  Although  final 
failure  may  again  occur  by  faulting,  relative  sliding  along  such  a  system  prior  to 
faulting  can  lead  to  dramatic  increases  in  ductility. 

In  an  effort  to  incorporate  the  slip  system  premise,  a  new  model  of  concrete 
was  constructed.  This  model,  the  details  of  which  are  given  in  Appendix  J.  includes 
an  arbitrary  number  of  slip  systems.  It  is  planned  to  conduct  validation  studies  to 
determine  the  simulation  capability  and  accuracy  of  the  theory  during  a  future  research 
period. 
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Section  4 


PLAIN  CONCRETE 

This  section  describes  the  research  conducted  on  the  theoretical  description  of 
the  behavior  of  plain  concrete.  The  discussion  covers  objectives,  accomplishments, 
general  mathematical  forms  and  validation  problems,  which  are  intended  to 
demonstrate  the  capability  and  accuracy  of  the  models. 

4.1  BASIC  OBJECTIVES 

The  ultimate  success  of  the  mixture  theory  for  reinforced  concrete,  discussed 
above,  hinges  very  strongly  on  the  ability  of  the  constitutive  model  of  plain  concrete, 
used  with  the  theory,  to  accurately  describe  the  behavior  of  plain  concrete  over  the 
wide  range  of  response  that  can  be  expected  in  practice.  Any  deficiencies  in  the  plain 
concrete  model  will  surely  be  reflected  —  and  possibly  enhanced  —  when  the  model  is 
used  in  conjunction  with  the  mixture  theory  to  describe  reinforced  concrete.  Clearly, 
the  accuracy  of  the  mixture  theory  can  never  be  greater  than  the  accuracy  of  the 
plain  concrete  model  used  with  it.  Because  of  this,  every  effort  should  be  made  to 
develop  an  accurate  model  of  plain  concrete. 

During  the  past  ten  years,  considerable  research  has  been  devoted  to 
understanding  and  modeling  the  constitutive  behavior  of  plain  concrete.  From  this  has 
come  a  variety  of  different  constitutive  models  for  plain  concrete  which  provide 
reasonably  accurate  descriptions  of  plain  concrete  behavior  for  stress  paths  which 
generally  do  not  differ  greatly  from  the  standard  paths  followed  in  the  usual 
laboratory  testing  Largely,  due  to  the  lack  of  available  appropriate  data,  little  proof¬ 
testing  of  the  models  has  been  done  to  examine  their  predictive  capabilities  for 
complex  stress  or  strain  paths  that  are  expected  to  occur  in  practice.  Many  of  the 
models  are  limited  to  the  stress  range  below  failure,  since  they  contain  no  provision 
for  treating  cracking.  A  few  of  the  models  do.  however,  attempt  to  treat  cracking 
and  post-cracking  behavior  as  well.  Generally  speaking,  however,  the  probems  of 
developing  a  constitutive  model  of  plain  concrete  capable  of  realistically  describing  the 
full  spectrum  of  behavior  from  pre-cracking  to  cracking  and  finally  post-cracking 
response  still  remains  basically  unresolved.  The  goal  of  the  present  research  is  to 
address  this  need. 

4.2  ACCOMPLISHMENTS. 

A  new  constitutive  model  of  plain  concrete  was  developed  and  validated  with 
remarkable  success  against  a  variety  of  proof  tests.  The  model  is  based  on  the 
endochronic  theory  of  plasticity  and  exhibits  the  important  features  of  nonlinear 
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concrete  behavior  for  stress  states  below  failure.  Details  of  the  model,  together  with 
example  validations,  are  given  in  Section  4.3.  An  approach  for  extending  the  model 
to  encompass  cracking  and  post-cracking  behavior  has  recently  been  made  by  Valanis 
(1987).  This  approach  will  be  the  basis  of  our  future  work  to  develop  an 
endochronic  plastic  fracturing  model  of  plain  concrete. 

In  the  course  of  exploring  various  approaches  for  mathematically  representing 
cracking  in  plain  concrete,  a  new  continuous  damage  model  was  developed  which  has 
certain  advantages,  some  of  which  are  not  found  in  other  models.  Details  and 
features  of  the  model  are  described  in  Section  4.4  and  more  fully  in  Appendix  A.  A 
number  of  examples  are  given  to  illustrate  the  various  features  of  this  new  model. 

In  view  of  the  Air  Force  s  interest  in  the  shock  loading  of  concrete  defense 
structures,  and  the  need  in  design  studies  for  a  rate-dependent  constitutive  model  of 
plain  concrete,  a  comprehensive  review  of  strain  rate  effects  in  plain  concretes  was 
conducted.  The  review  focussed  mainly  on  experimental  issues  related  to  strain  rate 
testing  and  on  the  extent  to  which  the  data  can  be  trusted.  These  issues  are 
discussed  in  Section  4. 

Finally,  certain  fundamental  issues  concerning  the  experimental  basis  for 
constitutive  models  for  brittle  materials  have  been  explored  and  are  discussed  in 
Section  4.6.  Problems  related  to  standard  quasi-static  laboratory  testing  methods  are 
treated.  Included  are  the  topics  of  strain  softening,  strain  hardening,  failure  states 
and  failure  modes. 

In  those  cases  where  the  research  accomplishments  have  been  reported  in  detail 
in  the  published  literature,  only  the  main  results  are  summarized  here  The  interested 
reader  is  referred  to  the  appropriate  papers  for  further  details.  The  references  for  this 
section  are  listed  at  the  end  of  the  section. 

4.3  A  NEW  CONSTITUTIVE  MODEL  FOR  PLAIN  CONCRETE 

A  new  constitutive  model  of  plain  concrete  was  developed  by  Valanis  and  Read 
(1985.1986)  which  appears  to  have  remarkable  capabilities  for  predicting  the  nonlinear 
inelastic  behavior  of  concrete  for  stress  states  below  failure.  The  model  is  formulated 
on  the  basis  of  the  endochronic  theory  of  plasticity  and.  as  such,  does  not  require  a 
yield  surface  nor  the  specification  of  loading  or  unloading  criteria,  as  in  classical 
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plasticity.  It  predicts  that  plastic  flow  will  occur  from  the  onset  of  loading,  a  feature 
which  makes  the  model  attractive  and  appropriate  for  describing  the  behavior  of  a 
material  such  as  concrete  which  does  not  exhibit  a  well-defined  yield  point.  Basically, 
the  model  is  isotropic,  rate-independent  and  satisfies  the  Clausius-Duhem  dissipation 
inequality.  It  realistically  portrays  the  major  features  of  nonlinear  inelastic  behavior 
exhibited  by  plain  concrete,  including  shear-volumetric  coupling,  effect  of  hydrostatic 
compression  on  shear  response,  hardening,  hysteresis  and  stress-path  dependence.  For 
the  case  of  isothermal  deformation,  the  governing  equations  of  the  new  model  are  as 
follows: 


a 


'0 


dfiP  =  dS  -  48/2/*  , 


de  -  da/K  , 


(4-1) 

(4-2) 

(4-3) 

(4-4) 

(4-5) 


dz#  =  dz/Fs  ,  dzH  =  dz/kF^  .  (4-6a,b) 

Here  j  denotes  the  deviatoric  stress  tensor,  a  is  the  hydrostatic  stress  (pressure).  gp 
represents  the  plastic  component  of  the  deviatoric  strain  tensor  g.  while  ep  is  the 
plastic  component  of  the  volumetric  strain  e.  Moreover,  n  and  K  are.  respectively, 
the  shear  and  bulk  moduli,  while  k  is  a  constant  which  determines  the  magnitude  of 
shear-volumetric  coupling.  The  double  bars  surrounding  a  symbol  denote  its  norm, 
while  single  bars  denote  absolute  value.  Furthermore.  Fs  and  FH  are.  respectively, 
shear  and  hydrostatic  hardening  functions.  In  addition,  z  denotes  the  intrinsic  time 
scale,  while  z?  and  zH  are.  respectively,  the  intrinsic  times  for  shear  and  hydrostatic 
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behavior.  Finally.  p(z)  and  ^(z)  are  weakly  singular  kernel  functions  satisfying  the 
condition  p(0)  =  ^(0)  =  **.  but  integrable  in  the  domain  0  <  <  *. 
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The  weakly  singular  nature  of  the  kernel  functions  is  a  crucial  feature  of  the 
model  for  two  reasons: 

(1)  It  provides  for  closure  of  hysteresis  loops  in  the  uniaxial  or  shear 
stress-strain  space,  however  small  they  may  be.  and 

(2)  it  ensures  that,  at  points  of  unloading  or  reloading,  the  response  is 
always  elastic. 

The  specific  forms  of  the  kernel  functions  used  in  the  present  model  are: 


-a  z 
r  s 


(4-7) 


'  H  # 

where  the  constants  A  .  B  a  ,  and  P  are  all  positive  and  finite.  With  such  forms,  it 
can  be  shown  (Valanis  and  bead.  1^84)  that  each  of  the  hereditary  integrals  in  (1) 
and  (2)  is  reducible  to  a  system  of  coupled  linear  ordinary  differential  equations:  this 
greatly  simplifies  the  computational  strategy  for  dealing  with  the  model. 

The  function  FH  reflects  the  effect  of  compaction  on  hydrostatic  behavior  and 
therefore  is  taken  to  depend  on  the  plastic  volumetric  strain.  ep.  The  function  F  .  on 
the  other  hand,  accounts  for  the  effect  of  hydrostatic  pressure  on  shear  response  and 
failure,  as  well  as  the  form  of  the  trace  of  the  failure  surface  in  the  ir-plane. 
Generally  speaking.  F$  will  depend  on  both  a  and  the  third  invariant  of  the  deviatoric 
stress  tensor.  J'.  In  the  present  study,  however,  where  the  behavior  of  concrete  for 
stress  histories  that  do  not  produce  significant  macrocracking  is  considered,  the  model 
predictions  are  not  likely  to  be  sensitive  to  the  fine  details  of  the  trace  of  the  failure 
surface  in  the  f-plane:  this  follows  from  the  fact  that  the  response  regime  of  interest 
is  not  in  the  close  proximity  of  the  failure  surface.  Therefore,  to  simplify  the  model. 
F^  is  taken  to  depend  only  on  a. 


An  extension  of  the  model  to  include  the  effect  of  J,  on  failure  is  discussed  in 
Appendix  B.  M 
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It  should  be  noted  that  the  endochronic  model  described  above  is  very  different 
from  previous  endochronic  models  of  concrete  developed  by  Bazant  and  Bhat  (1976) 
and  Bazant  and  Shieh  (1978)  on  the  basis  of  earlier  versions  of  endochronic  theory. 

In  the  model  discussed  here,  proper  closure  of  hysteresis  loops  is  guaranteed,  and  the 
instantaneous  response  at  points  where  loading,  unloading  and  reloading  begin  is 
elastic.  The  basic  plastic  flow  properties  of  the  model  have  been  explored  recently  by 
Trangenstein  and  Read  (1982)  and  Murakami  and  Read  (1986:1987). 

i 

The  foregoing  endochronic  model  was  applied  to  an  extensive  set  of  laboratory 
data  recently  generated  by  the  University  of  Colorado,  using  a  true  triaxial  device 
(Scavuzzo.  et  a/..  1983).  The  test  program  consisted  of  six  different  series  of  stress- 
controlled  tests,  each  of  which  was  designed  to  explore  a  different  facet  of  concrete 
behavior.  Altogether,  the  response  of  concrete  to  over  45  complex  stress  paths  was 
investigated.  In  all  cases,  the  stress  paths  were  such  that  no  significant  cracking 
took  place  during  the  tests.  j 

The  model  was  fit  to  a  very  small  subset  of  the  data,  after  which  it  was  proof- 
tested  by  exercising  it  around  over  20  complex  stress  paths  and  predicting  the 
corresponding  deformation  histories.  None  of  the  data  from  the  complex  stress  path 
tests  selected  for  proof-testing  were  used  in  fitting  the  model  parameters  and  no 
optimization-techniques  were  employed. 

I 

The  proof-tests  were  remarkably  successful  and  revealed  the  powerful  predictive  I 

capability  of  the  model.  Figures  4.1  to  4.4.  for  example,  illustrate  the  results  of  the  ! 

proof-tests  for  four  of  the  many  complex  stress  paths  considered.  The  imposed  stress  J 

paths  are  depicted  at  the  top  of  each  figure.  It  is  seen  that  the  comparisons 
between  the  model  predictions  and  the  measurements  are  excellent,  considering  the 
relatively  large  data  scatter  normally  expected  with  concrete.  The  ability  of  the  model 
to  accurately  describe  the  behavior  of  plain  concrete  around  the  complex  stress  paths 
not  used  in  fitting  the  model  reveal  its  remarkable  predictive  capabilities.  Various 
other  examples  of  the  proof-tests  can  be  found  in  the  papers  by  Valanis  and  Read 
(1985.1986). 

The  above  endochronic  model  has  several  limitations  which  restrict  its  application 
to  stress  states  below  failure.  First,  the  model  does  not  have  the  capability  to  I 

describe  macrocracking  and  its  effect  on  material  behavior.  Secondly,  since  the  model  J 
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is  isotropic,  it  cannot  account  for  the  anisotropy  that  can  be  expected  to  develop 
when  significant  macrocracking  with  preferred  orientation  occurs.  Finally,  when  there 
is  shearing  in  the  presence  of  a  fixed  hydrostatic  pressure,  the  model  exhibits  only 
compaction.  Therefore,  dilatant  behavior  usually  observed  when  cracked  material  is 
sheared  cannot  be  described  by  this  model. 

Recently,  however.  Valanis  (1987)  has  developed  an  approach  for  generalizing  the 
endochronic  concrete  model  which  overcomes  all  of  the  above  limitations.  By  using  a 
mapping  between  undamaged  and  damaged  material,  he  found  the  form  taken  by  the 
endochronic  concrete  model  in  the  presence  of  cracking.  This  new  model  has  been 
termed  the  endochronic  plastic  fracturing  model  (Valanis.  1987).  The  effect  of 
cracking  is  reflected  through  a  second  order  integrity  tensor  g.  which  is  governed  by 
an  evolution  equation  that  relates  the  cracking  to  the  tensile  strain  history.  The 
plasticity  of  the  non-cracked  material  is  described  by  the  endochronic  concrete  model. 
The  resulting  model  appears  to  have  the  desired  features  for  describing  the  full  gamut 
of  concrete  behavior,  including  plastic  flow  and  cracking.  It  exhibits  anisotropy  if 
cracking  develops  in  preferred  directions,  as  well  as  dilatancy  during  cracking.  The 
anisotropy  is  accounted  for  in  a  clever  way  through  the  use  of  the  second-order 
tensor  g,  rather  than  by  a  fourth-order  tensor,  as  is  the  usual  approach.  The  use  of 
a  second-order  tensor  for  this  purpose  greatly  reduces  the  analytic  complexity  of  the 
model  and  simplifies  the  experimental  determination  of  the  material  parameters.  Crack 
initiation  and  growth  are  governed  by  the  evolution  equation  for  g.  In  its  present 
form,  the  model  is  restricted  to  distributed  damage,  since  it  does  not  contain  a 
characteristic  length.  The  endochronic  plastic  fracturing  model  appears  to  offer 
considerable  promise  for  describing  the  general  behavior  of  plain  concrete  and.  as  such, 
will  be  the  focus  of  our  future  efforts  in  this  direction. 

Many  of  the  important  questions  faced  by  the  Air  Force  today  regarding  the 
hardness  of  concrete  defense  structures  requires  a  knowledge  and  understanding  of  the 
likelihood  that  a  given  structure  will  fail  under  a  given  loading  environment.  To 
address  this  issue  in  the  most  cost-effective  manner,  analytic  methods  must  be 
available  to  allow  system  designers  to  systematically  vary  design  parameters  in  an 
expedient  and  efficient  manner  so  that  an  optimum  structural  design  can  be  achieved. 
Finite  element  methods  are  now  available  for  this  purpose,  but  the  results  obtained 
with  these  methods  can  only  be  as  good  as  the  constitutive  relations  used  in 
conjunction  with  them.  To  examine  the  hardness  of  concrete  structures  requires  a 
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constitutive  model  of  plain  concrete  which  not  only  is  capable  of  describing  the 
inelastic  behavior  of  intact  unfractured  material  but  is  capable  of  describing  developing 
damage  and  cracking  and  their  effect  on  the  overall  constitutive  behavior.  For  the 
endochronic  concrete  model  to  be  useful  in  addressing  these  important  issues,  it  must 
be  able  to  treat  elastic-plastic  and  cracking  behavior.  In  the  following  section,  our 
initial  efforts  to  develop  an  advanced  cracking  model  for  this  purpose  are  described. 

4.4  A  NEW  CONTINUOUS  DAMAGE  MODEL. 

The  approach  that  we  have  taken  to  develop  a  cracking  model  appropriate  for 
use  in  finite  element  framework  utilizes  the  notion  of  an  equivalent  homogeneous 
continuum  element  to  simulate  the  response  of  a  finite  volume  of  cracked  or  damaged 
material,  as  depicted  in  Figure  4.5.  In  effect,  this  means  that  instead  of  attempting 
to  track  and  treat  individual  cracks  in  a  finite  volume  of  material  and  to  account  for 
their  effect  on  the  overall  response  of  the  element,  we  postulate  an  equivalent 
homogeneous  continuum  element  whose  constitutive  properties  are  such  that  they 
provide  the  same  response  to  applied  loads  as  the  actual  cracked  material  element. 
The  problem  then  is  to  find  the  appropriate  constitutive  relation  for  the  equivalent 
continuum  element. 

With  the  above  goal  in  mind,  a  new  and  very  promising  model  for  the 
equivalent  homogeneous  continuum  was  developed  under  this  contract  by  Valanis 
(1985).  Although  the  model  described  in  this  section  applies  to  brittle  elastic  solids, 
the  basic  approach  is  extendable  to  elastic-plastic  fracturing  solids,  as  noted  in  the 
preceding  section.  The  basic  form  of  the  model  for  elastic-fracturing  solids  is 

2  =  H-8) 

where  g  and  £  denote,  respectively,  certain  averages  of  the  stress  and  strain  taken 
over  the  volume  of  the  material,  while  g.  termed  the  integrity  tensor,  describes  the 
orientation  of  damage,  as  well  as  its  extent  and  character.  It  should  be  noted  that 
the  model  is  also  capable  of  describing  the  evolution  of  damage  in  a  material  that  has 
known  pre-existing  damage  in  the  form  of  oriented  cracks. 

A  unique  feature  of  the  model,  which  sets  it  apart  from  other  cracking  models, 
is  that  it  can  describe  the  development  of  anisotropy  due  to  cracking  without  requiring 
the  introduction  of  a  fourth-order  stiffness  tensor  into  the  formulation.  Through  a 


Figure  4.5  Representation  of  a  cracked  multi-phase  material  element  by  an  equivalent 
homogeneous  continuum  element. 


clever  formulation  described  below,  general  anisotropy  is  represented  by  a  constitutive 
equation  that  contains  only  second  order  tensors.  The  resulting  model  can  describe 
the  full  spectrum  of  cracking,  from  the  accumulation  of  microdamage  up  to  complete 
fracture.  In  addition,  it  appears  that  the  model  can  describe  the  standard  modes  of 
cracking,  as  well  as  the  "splitting  mode”  under  uniaxial  compression  (Horii  and 
Nemat- Nasser.  1985).  a  mode  which  most  of  the  existing  fracture  models  cannot  treat. 
Finally,  the  model  predicts  dilatancy  during  cracking,  in  accord  with  experimental 
observation,  as  we  shall  show  in  the  sequel.  In  the  following,  a  brief  outline  of  the 
model  will  be  given.  A  more  extensive  discussion  of  the  model  is  given  in 
Appendix  A. 

The  fracturing  theory  is  formulated  below  within  the  context  of  irreversible 
thermodynamics,  although  Valanis  (1987)  has  recently  shown  that  the  same  model  can 
be  developed  by  considering  the  mapping  from  undamaged  to  damaged  material.  We 
now  introduce  the  free  energy  j  per  unit  volume.  Since  fracture  produces  a 
reduction  in  material  integrity,  an  integrity  tensor  g  is  introduced  which  has  the 
following  properties:  It  is  a  second-order  symmetric  tensor  which  is  equal  to  the  unit 
tensor  £  in  the  undamaged  virgin  material  and  to  the  null  tensor  in  the  fully  fractured 
material,  i.e..  material  which  cannot  support  tensile  stress  in  any  direction.  We  thus 
set 

1  •*[*.£}  (4-9) 


where  isothermal  conditions  and  small  strains  are  assumed. 

In  a  thermodynamic  sense.  £  plays  the  role  of  an  internal  variable,  in  which  case 
we  have  the  relations 


,  .M  n  =  _  9* 

2  9£  '  »  3£ 


(4-9a) 


where  g  is  the  stress  tensor  and  Q  the  internal  force,  dual  to  g.  which  drives  the 
fracturing  process. 


In  view  of  the  properties  of  the  following  relations  must  hold 

0  ,  *(0,4}  =  0  , 


(4-9b) 


*  This  nomenclature  is  preferred  to  free  energy  density,  in  view  of  the  concept  of 
"irreducible"  material  volume  which  must  be  introduced  when  dealing  with 
heterogeneous,  multi-phase  materials,  cich  as  plain  concrete. 


the  first  meaning  that  a  fully  fractured  material  cannot  contain  free  energy  and  the 
second  that  an  unstrained  material  must  have  zero  free  energy  (both  relative  to  the 
reference  state).  The  relations  below  must  also  hold: 


2(2,9)  =  9  .  2(9,*)  =  9 
8(2.9)  =  9 


(4-9c) 


The  first  two  equations  stipulate  that  the  stress  must  vanish  in  the  fully  failed 
material  and  that,  since  the  material  remains  elastic  in  the  damaged  state,  the  stress 
will  vanish  at  zero  strain.  The  last  equation  stipulates  that  the  internal  force  Q 
causing  fracture  must  vanish  in  the  fully  failed  material. 

The  free  energy  f  is  expanded  in  a  Taylor  series  in  g  and  terms  no  higher  than 
quadratic  are  retained.  In  addition,  it  is  assumed  that  f  depends  on  £  in  a  quadratic 
manner.  On  this  basis,  it  follows  that 

*  *  Cijk(L£ij£k<L  ■  M-*0 

where  £(^)  is  purely  quadratic  in  If  the  material  is  isotropic  in  the  virgin  state,  it 
may  be  represented  in  terms  of  the  outer  products  of  the  unit  tensor  and  £.  Since  £ 
is  purely  quadratic  in  it  must  have  the  form 

Cijk2  =  X  *i/k«.  *  2/i  *ik*j2  '  (4"9®) 


in  view  of  the  symmetries  imposed  upon  it  by  the  symmetry  of  g  via  Eq.  (4-9d). 
The  constants  in  Eq.  (4-9e)  must  indeed  be  the  Lam6  constants  of  the  virgin  material 
since 


C. 


i  jk£. 


=  X 


6. -6.0 

i  j  kH 


sikV 


(4-10) 


when  *  =  jj.  Therefore,  from  Eqs.  (l-9d)  and  (4-9e)  we  can  write 
*  *  2  M*ij£ij)2  *  <4  *ik*j2.£ikEj2.  - 


I 

(4-11) 


*  This  assumption  has  been  validated  through  a  completely  independent  (non¬ 
thermodynamic)  approach  based  upon  a  mapping  from  undamaged  to  damaged 
material  (Valanis,  1987).  ^ 


which  can  be  combined  with  the  first  of  Eqs.  (4-9a)  to  give 

ffij  =  X  *  2/4  (4"12^ 

To  complete  the  specification  of  the  model,  it  is  necessary  to  establish  an  association 
between  the  mathematical  structure  of  £  and  observed  physical  damage.  This  is 
discussed  below. 

In  the  initial  version  of  the  model,  attention  was  focussed  on  the  case  in  which 
developing  damage  is  totally  directed,  that  is.  an  advancing  crack  produces  damage 
only  in  a  direction  normal  to  the  crack.  Attention  to  this  case  was  originally 
motivated  by  the  observation  that  in  cementitious  materials,  such  as  plain  concrete, 
when  subjected  to  simple  stress  histories,  the  crack  patterns,  while  locally  quite 
tortuous,  nevertheless  are  nominally  normal  to  the  direction  of  maximum  tensile  strain, 
as  illustrated  in  Figure  4-6.  With  these  notions  in  mind,  an  evolution  equation  for  £ 
was  introduced  on  the  basis  of  the  following  considerations. 

When  the  material  is  not  strain  rate  sensitive,  an  increment  of  tensile  strain  will 
produce  an  increment  of  damage.  Specifically,  given  an  increment  in  strain.  Ag.  let 
Aea  be  its  eigenvalue  and  n®  its  eigenvectors  (a  =  1.2.3).  If  Ag  and  g  are  coaxial, 
then  Ac®  will  constitute  an  increase  in  tensile  strain  if  Ae®  >  0  and  e®  l  0.  where 
e®  denote  the  eigenvalues  of  g.  In  the  more  general  case  where  Ag  and  |  are  not 
coaxial.  Ae®  is  said  to  be  an  increase  in  tensile  strain  if  Aea  >  0  and  e..n®n®  s  e£  l 
0  (a  not  summed).  Now  let  Ag  be  the  change  in  the  integrity  tensor  g  due  to  the 
increment  dg.  and  let  be  the  eigenvector  of  Ag.  With  the  above  in  mind,  we  now 
introduce  the  following  two  propositions: 

(i)  The  tensors  Ag  and  Ag  are  coaxial.  The  eigenvectors  of  Ag  are 
then  equal  to  the  eigenvectors  of  Ag.  i.e.. 

0®  =  cj  («-13) 

(ii)  The  eigenvalues  of  Ag  are  assumed  to  be  related  to  the  eigenvalues 
of  Ag  according  to  the  following  expression  which  has  its  roots  in 
the  physics  of  annihilation: 


Figure  4.6  Observed  fracture  patterns  in  plain  concrete. 
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kAe“.  Ae  >  0.  e  “  £  0 
n 


10  otherwise 


and 
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In  the  above  equations,  k  and  m  are  positive  material  constants.  In  view  of  the 
above  developments,  we  can  write 


a 


ananfl 

ninj 


(4-17) 


By  adjoining  this  evolution  equation  for  £  to  Eq.  (4-12).  the  constitutive 
description  for  a  brittle  elastic  solid  which  develops  directed  damage  is  complete,  once 
the  material  parameters  X./j.k  and  m  have  been  specified.  Applications  of  this  theory 
to  several  cases  involving  homogeneous  and  non-homogeneous  triaxiai  stress  fields  are 
described  in  Appendix  A.  As  shown  there,  the  theory  leads  to  results  which  are  in 
basic  agreement  with  experimental  observation  and  with  theoretical  solutions  from 
linear  elastic  fracture  mechanics. 


Recently,  in  an  effort  to  adapt  the  model  specifically  to  plain  concrete,  a  more 
general  form  of  the  evolution  equation  for  g  was  developed  by  Valanis  (1986)  which 
accounts  for  both  directed  and  non-directed  damage.  This  representation  of  g  was 
motivated  by  the  results  from  the  careful  study  done  by  Diamond  and  Bentur  (1984) 
which  clearly  showed  that  directed  cracking  in  cementitious  materials  is  usually 
accompanied  by  the  formation  of  small,  randomly  oriented  cracks  that  develop  near 
the  crack  tip  (see  Figure  4.6d).  From  a  statistical  standpoint,  the  damage  produced 
by  these  small  cracks  can  be  considered  isotropic,  and  we  denote  such  damage  by  the 
scalar  Q.  Thus  Q  reduces  the  stress  carrying  capacity  of  a  material  by  the  same 
degree  in  all  directions.  As  a  result,  the  integrity  tensor,  which  we  now  denote  by 
gjj.  must  be  reduced  by  the  presence  of  R  and  hence  must  be  a  state  function  of  Q. 
Accordingly,  to  account  for  both  directed  and  isotropic  damage,  we  write 


(4-18) 
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*0  -  <Qk<“)  . 


subject  to  the  following  conditions  which  must  always  apply: 


<(0,0)  .  <  ,  <(i,{*)  •  2 


(4-19) 


Here,  g  denotes  the  unit  tensor  and  is  defined  by  Eq.  (4-15).  To  satisfy 
conditions  (4-19).  we  take  g  to  be  of  the  form: 


<0  *  «(0)<(f“)  . 


(4-20) 


where  g  is  the  integrity  tensor  for  directed  damage  discussed  above  and  w(Q)  is  a 
monotonically  decreasing  function  of  Q  such  that  w(0)  =  1  and  w(l)  =  0  A  suitable 
function  which  satisfies  the  above  constraints  is 


*(0)  *  (1  -  Q)m 


(4-21) 


where  a  l  0.  Thus,  the  integrity  tensor  g  considered  in  the  previous  section  is  now 
modified  by  the  presence  of  isotropic  damage,  as  specified  by  Eqs.  (4-20)  and  (4-21) 
The  constitutive  equation  thus  becomes 


»ij  £kl  *  2*ikVkt) 


(4-22) 


where  A  and  ft  have  the  same  meaning  as  before.  Inasmuch  as  an  expression  for  the 
evolution  of  g  with  damage  was  given  earlier,  an  evolution  equation  for  u  is  needed  in 
order  to  complete  the  constitutive  representation:  this  is  considered  below. 

Physical  considerations  suggest  that  isotropic  damage  in  cementitious  materials  is 
caused  primarily  by  shearing  which  gives  rise  to  local  decohesion  between  the  mortar 
and  aggregate  and/or  the  paste  and  the  sand  particles.  Therefore,  following  similar 
reasoning  that  gave  rise  to  Eq.  (4-14).  we  take  dQ  to  be  proportional  to  a  positive 
monotonic  function  f  of  the  isotropic  integrity  1  -  Q.  i.e.. 


dQ  =  f(l  -  Q)d£ 


where 


d£  *  c| |dSl 


(4-23) 


(4-24) 


Here.  C  >  0,  g  is  the  deviatoric  strain  tensor,  and  the  double  bars  around  a 
symbol  denote  its  Euclidean  norm.  Bearing  in  mind  that  the  ultimate  value  of  Q  is 
unity,  a  suitable  form  of  f  is 

f  =  (1  -  Q)b  (4-25) 

where  b  >  0.  Upon  substitution  of  Eq.  (4-25)  into  Eq.(4-23).  and  solving  the 
resulting  differential  equation,  we  obtain  the  result: 

1 

1  -  0  -  [1  -  (1  -  b)f]1_b  ,  (4-26) 

which  is  similar  in  form  to  Eq.  (A-26)  given  earlier  in  connection  with  directed 
damage.  For  the  special  case  of  b  =  1.  the  above  expression  reduces  to 

1  -  0  =  e"b*  (4-27) 

Note  that  C  in  Eq.  (4-24)  plays  the  role  of  an  isotropic  fracture  toughness  parameter 
which  influences  the  rate  of  isotropic  degradation.  The  specification  of  the 
constitutive  model  that  accounts  for  both  directed  and  isotropic  damage  is  thus 
complete. 

To  illustrate  the  capability  of  the  above  model  to  describe  the  observed  behavior 
of  a  real  cementitious  material,  let  us  consider  the  case  of  axial  compression  of  a 
plain  concrete  prism,  which  has  been  studied  experimentally  by  Van  Mier  (1984).  We 
denote  by  a  and  the  stress  and  strain  in  the  axial  direction  and  by  e  (=e  )  the 

strain  in  the  transverse  direction.  For  this  case,  it  can  be  shown  (see  Appendix  A) 
that  the  model  discussed  above  leads  to  the  following  set  of  coupled  equations: 

2a 

=  E  cja  -  (a  -  b)  C*e2)1-b  (4-28) 

1 

-r/Cl  =  c2{l  -  (1  -  m)k  ej1""  (4-29) 

Equations  (4-28)  and  (4-29)  were  applied  to  the  data  reported  by  Van  Mier  (1984) 
for  the  case  of  axial  compression  of  a  prism  having  a  height  of  50  mm  and  lateral 
dimensions  of  100  mm  From  these  data,  the  following  values  of  the  parameters  in 
the  above  equations  were  determined. 
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~  =  1.5  (l-b)C*  =  40 

Figures  4.7  to  4.9  illustrate  the  manner  in  which  the  model,  with  the  values  of  the 
parameters  listed  above,  describes  the  data  from  Van  Mier  (1984).  Figure  4.7.  for 
instance,  shows  the  dependence  of  the  axial  strain  on  the  lateral  strain 
revealing  the  excellent  agreement  between  the  model  and  the  data.  Figure  4.8  depicts 
the  relation  between  the  axial  strain  and  the  volumetric  strain  ek((,  showing  that 
the  material  initially  experiences  some  small  compaction  before  it  becomes  dilatant 
with  increasing  compressive  strain.  Finally.  Figure  4.9  shows  the  dependence  of  the 
ratio  a.  I  (a.)  on  the  axial  strain  e..  Here,  lo.)  denotes  the  maximum  value  of 
a  measured  in  the  tests.  This  figure  reveals  that  the  model  predicts  softening 
beyond  the  peak  stress  which  is  in  reasonable  agreement  with  the  data.  Overall,  the 
agreement  between  the  model  and  the  data  is  considered  quite  good. 


Remark* 

The  fracture  model  presented  in  this  section  and  further  elaborated  on  in 
Appendix  A  is  currently  in  the  development  stage.  Based  on  the  evidence  presented 
in  the  preceding  sections,  the  model  clearly  has  a  number  of  attractive  and  unique 
features  which  make  it  worthy  of  continued  development. 

•  It  can  treat  material  that,  in  the  initial  state,  is  either  uncracked 
and  isotropic  or  has  preexisting  cracks. 

•  It  model  has  the  ability  to  develop  anisotropy  as  damage  occurs  in 
preferred  directions.  Initially,  the  model  is  isotropic. 

•  The  development  of  anisotropy  is  accounted  for  in  the  model  in  a 
very  convenient  manner  involving  only  second  order  tensors.  By 
comparison,  the  usual  anisotropic  model  requires  a  fourth-order 
tensor,  which  greatly  increases  the  complexity  of  the  model,  as  well 
as  the  evaluation  of  the  material  parameters. 
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Figure  4.7  Relation  between  axial  and  tranaverse  strains  under  conditions  of 
simple  compression. 


Figure  4.8  Dilatancy  during  simple  compression. 


Axial  stress-strain  relation  showing  softening  in  compression. 
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•  The  model  describes  the  full  spectrum  of  developing  damage, 
beginning  with  the  gradual  accumulation  of  distributed  damage  and 
continuing  up  to  complete  fracture. 

•  In  addition  to  the  standard  modes  of  fracture,  the  model  can 
describe  the  splitting  modes  under  uniaxial  compression,  which  most 
other  fracture  models  are  unable  to  treat. 

*  During  the  development  of  increasing  damage,  the  model  initially 
exhibits  compaction  and  then  dilatancy.  in  agreement  with 
experimental  observations  on  plain  concrete. 

*  The  model  appears  to  be  attractive  for  use  in  finite-element  (FE) 
analysis.  In  the  present  study,  we  developed  a  computer  program 
for  the  model  and  incorporated  it  into  a  FE  code,  which  was  then 
used  to  analyze  several  fracture  problems  described  in  Appendix  A. 

In  its  present  form,  the  above  fracture  model  applies  specifically  to  elastic- 
fracturing  solids.  The  approach  has.  however,  been  recently  extended  by  Valanis 
(1987)  to  plastic  fracturing  solids,  as  noted  earlier.  In  our  future  effort,  we  will 
pursue  this  new  approach  with  the  goal  of  developing  an  endochronic  plastic  fracturing 
model  for  plain  concrete.  The  damage  will  be  accounted  for  in  the  model  by  the  £ 
tensor,  while  the  endochronic  concrete  model  will  describe  the  concrete  plasticity.  The 
damage  portion  of  the  model  will  be  calibrated  from  experimental  data  obtained  for 
plain  concret  at  low  pressures,  where  the  deformation  is  primarily  due  to  cracking. 
Experiments  conducted  at  high  pressures,  above  the  brittle-ductile  transition  where  the 
deformation  is  mainly  due  to  plastic  flow,  will  be  used  to  calibrate  the  plastic  portion 
of  the  model.  Once  the  model  has  been  developed  in  this  manner,  it  will  then  be 
proof-tested  against  data  obtained  at  intermediate  pressures. 

It  should  also  be  noted  that,  as  presently  formulated,  the  model  is  not  sensitive 
to  specimen  size.  As  a  result,  its  application  is  essentially  restricted  to  the  case  in 
which  damage  develops  in  a  uniformly  distributed  manner  over  a  specimen.  In  this 
case,  the  behavior  does  not  depend  on  the  size  of  the  specimen.  On  the  other  hand, 
if  the  damage  is  not  uniformly  distributed  but  localizes,  a  size  effect  emerges  and 
then  must  be  reflected  in  the  model.  The  notion  of  a  size  effect,  and  its 
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incorporation  into  fracture  models  for  plain  concrete,  is  currently  a  topic  of 
considerable  research  and  debate  within  the  mechanics  community.  We  plan  to 
pursue  this  issue  in  our  future  efforts. 

Finally,  when  the  applications  of  interest  involve  rapid  impulsive  loading  of 
concrete,  as  is  the  case  here,  the  model  should  reflect  the  possible  effect  of  strain-rate 
dependence  on  the  behavior.  Toward  this  end.  a  critical  review  of  the  existing  state- 
of-the-art  regarding  strain  rate  effects,  and  their  measurement,  for  plain  concrete  was 
prepared  as  part  of  this  study.  A  discussion  of  the  strain-rate  issues  and  our 
findings  are  presented  in  the  following  section. 

4.5  STRAIN-RATE  EFFECTS  IN  PLAIN  CONCRETE. 

Understanding  the  manner  in  which  concrete  defense  structures  respond  to 
rapidly  applied,  impulsive  loadings  —  such  as  produced,  for  example,  by  an  explosive 
source  --  continues  to  be  a  major  concern  for  the  Air  Force.  Generally  speaking, 
typical  threat  scenarios  can  be  expected  to  induce  peak  strain  rates  in  concrete 
structures  that  range  from  10  to  10  sec  . 

One  of  the  complicating  factors  that  arises  when  attempts  are  made  to  model 
such  behavior  is  that  concrete  responds  differently  under  high  rates  of  loading  than 
when  loaded  quasi-statically.  Unfortunately,  virtually  all  of  the  existing  data  used  to 
develop  and  calibrate  concrete  constitutive  models  have  been  obtained  from  quasi-static 
laboratory  tests  The  existing  evidence,  which  is  inconsistent  and  confusing  as  we 
will  show,  nevertheless  indicates  that,  as  the  rate  of  loading  increases,  concrete 
becomes  more  brittle  (less  ductile)  and  its  strength  increases.  Because  of  this, 
constitutive  models  calibrated  to  quasi-static  laboratory  data  are  not  likely  to  be 
adequate  for  describing  the  behavior  of  plain  concrete  at  high  rates  of  loading. 

It  appears  that  little  is  known  about  the  mechanism  that  causes  rate  dependence 
in  plain  concrete.  Some  of  the  outstanding  questions  which  need  to  be  addressed 
experimentally  include  the  following:  Is  the  rate  effect  due  mainly  to  cracking  and  the 
fact  that  cracks  propagate  at  a  finite  velocity?  Or  is  it  due  to  plastic  flow?  Is  the 
rate-effect  observed  both  below  and  above  the  brittle-ductile  transition?  Answers  to 
these  questions  would  prove  most  valuable  in  future  effort  to  introduce  rate  effects  in 
the  model  in  an  appropriate  manner. 


With  the  goal  of  eventually  introducing  rate-dependence  into  the  endochronic 
concrete  model,  we  conducted,  as  part  of  this  study,  a  critical  review  of  the  existing 
literature  on  strain  rate  effects  in  concrete,  focussing  primarily  on  the  following 
question:  Does  the  existing  data  base  on  strain  rate  effects  in  plain  concrete  have 
credibility  and  accurately  portray  the  true  strain  rate  dependence,  or  are  the  rate 
effects  that  have  been  inferred  in  the  past  from  various  experimental  devices  primarily 
due  to  system  effects,  inhomogeneous  deformation  of  the  test  specimen  or  overly 
simplified  methods  of  data  analysis?  In  view  of  the  importance  to  the  Air  Force  of 
understanding  the  dynamic  response  of  plain  concrete  in  the  range  of  strain  rates  from 
10  1  to  102  sec’1,  the  review  was  primarily  focussed  on  this  range. 

Inasmuch  as  detailed  discussions  of  the  results  from  this  critical  assessment  of 
rate  effects  in  concrete  have  been  given  in  a  paper  (Hegemier  and  Read.  1985)  and  in 
a  topical  report  (Read.  1984).  we  shall  summarize  the  important  findings  below. 
Those  interested  in  further  details  are  referred  to  the  references  cited  above,  where 
more  extensive  discussions  are  given. 

Drop-Hammer  Tests. 

With  two  exceptions  (Birkimer.  1968:  Malvern,  et  a!.,  1985).  all  of  the  data  in 
the  literature  on  the  strain  rate  dependence  of  plain  concretes  within  the  range  from 
10  1  to  102  sec’1  have  been  obtained  from  drop-hammer  tests.  In  such  tests,  which 
are  schematically  depicted  in  Figure  4-10.  a  concrete  specimen  is  impacted  on  its 
upper  flat  surface  by  a  flat  heavy  hammer  which  is  released  from  a  prescribed  height 
and  falls  under  the  action  of  gravity  The  test  specimen  rests  as  shown  on  an  anvil, 
which  is  usually  fabricated  from  a  high  strength  steel  and  remains  elastic  during  the 
impact  process.  By  varying  the  height  of  drop,  the  impact  velocity  —  and  hence  the 
nominal  strain  induced  in  the  specimen  —  can  be  changed.  Although  not  shown  in 
Figure  4-10,  there  is  a  braking  device  which  catches  the  hammer  on  rebound  after 
impact  and  prevents  it  from  making  a  second  contact  with  the  impacted  specimen,  or 
from  crushing  the  remains  of  a  failed  specimen. 

Figure  4-11  summarizes  the  results  from  all  the  known  drop-hammer  studies 
conducted  on  plain  concrete  which  have  been  reported  in  the  literature.  Here,  the 
effect  of  strain  rate  on  the  ratio  of  the  dynamic  strength.  FQ.  to  the  static  strength. 
F  .  is  shown.  All  of  the  data  points  which  fall  in  the  range  from  10’6  to  10’2  sec  1 


Figure  4.10  Typical  features  of  drop-hammer  devices. 


igure  4.11  The  effect  of  strain  rate  on  the  ratio  of  the  dynamic  to  static 
strengths  of  plain  concrete. 


were  obtained  with  conventional  hydraulic  testing  devices,  while  the  data  between  10’1 
and  102  sec  *  are  from  drop-hammer  tests,  except  for  the  single  data  point  from 
Birkimer's  investigation,  which  was  obtained  with  a  specially  designed  hydraulic  testing 
machine. 

An  inspection  of  Figure  4-11  reveals  that  there  are  large  inconsistencies  between 
the  results  from  the  various  investigations.  Watstein's  data,  for  example,  show  the 
dynamic  strength  of  concrete  increasing  sharply  with  strain  rate,  while  Atchley  and 
Furr  find  a  greater  increase  in  strength  at  the  low  strain  rates  and  less  at  the  high 
rates.  Birkimer's  data  reveal  that  the  strength  becomes  independent  of  strain  rate  for 
strain  rates  above  10’1  sec’1,  which  is  in  direct  conflict  with  Watstein's  data  which 
show  a  very  dramatic  effect  of  strain  rate.  Hughes  and  Watson,  on  the  other  hand, 
find  virtually  no  effect  of  strain  rate  up  to  8  sec’1,  while  Hughes  and  Gregory 
reported  that.  “No  correlation  was  apparent  (from  their  data)  between  the  average 
rate  of  straining  and  either  uniaxial  compressive  strength  or  mix  proportions  in  the 
mixes  used." 

The  results  displayed  in  Figure  4-11  clearly  illustrate  the  disturbing 
inconsistencies  in  the  existing  data,  especially  for  the  range  of  strain  rates  from  10’1 
to  102  sec’1.  Although  we  will  not  attempt  to  go  into  the  details  here,  which  can  be 
found  in  the  references  cited  near  the  beginning  of  this  section,  the  data  shown  in 
Figures  4.11  appear  to  be  strongly  influenced  by  the  differences  in  drop-hammer  design 
between  the  various  investigations,  gauge  placement,  specimen  shape,  size  and  end 
capping,  as  well  as  methods  of  data  interpretation.  In  addition,  it  is  also  possible,  in 
view  of  the  conclusions  reached  by  Janach  (1976)  and  Young  and  Powell  (1979).  that 
increased  confinement  of  the  specimen  due  to  radial  inertia  effects  present  at  the 
higher  rates  of  loading  can  give  rise  to  an  "apparent”  strain  rate  effect.  Also,  in  all 
probability,  the  test  specimens  do  not  deform  in  a  homogeneous  manner  during  the 
gauge  read  time,  which  violates  one  of  the  crucial  requirements  for  inferring 
constitutive  properties. 

In  view  of  the  above  considerations,  it  appears  that  there  are  a  number  of 
serious  complications  with  the  drop-hammer  test  that  make  the  extraction  of  credible 
constitutive  property  information  from  it  very  difficult.  Because  of  this,  the  strain 
rate  effect  usually  inferred  from  such  tests  is  probably  more  apparent  than  real,  due 
in  large  part  to  system  effects.  Further  studies,  both  analytically  and  experimentally, 
need  to  be  done  to  address  the  issues  raised  above  and  to  determine  the  usefulness 
of  drop-hammer  tests  for  determining  constitutive  properties. 


Split  Hopkinson  Bar  Tests. 

Another  test  device  which  has  recently  been  used  to  investigate  strain  rate 
effects  in  plain  concrete  is  the  split  Hopkinson  bar  (SHB).  In  a  SHB  test,  a 
cylindrical  specimen  is  positioned  on  axis  between  two  long  circular  high  strength  steel 
bars  of  the  same  diameter,  as  shown  in  Figure  4-12.  A  compressive  pulse  is 
introduced  into  the  input  bar  by  various  means  and  travels  through  the  test  specimen 
into  the  output  bar.  Strain  gauges  mounted  on  the  input  and  output  bars  near  the 
test  specimen  measure  the  strain  histories.  If  it  is  assumed  that  the  deformation  of 
the  specimen  during  the  read  time  of  the  gauges  is  homogeneous  and  due  to  uniaxial 
stress,  then  the  dynamic  constitutive  properties  of  the  specimen  can  be  inferred  from 
the  measured  strain  histories.  The  homogeneity  assumption  is  therefore  the  key  issue. 
Parenthetically,  we  note  that  the  type  of  loading  applied  to  a  specimen  by  the  SHB 
is  the  same  as  that  produced  by  a  drop  hammer  device,  i.e..  nominally  uniaxial  stress. 

The  SHB  has  been  widely  used  in  the  past  to  study  the  strain  rate  dependence 
of  various  metals  and  other  ductile  solids  and.  for  those  materials,  it  is  generally 
considered  to  be  a  credible  approach.  For  brittle  materials,  however,  the  use  of  the 
SHB  has  raised  serious  questions  which  to  date  have  not  been  resolved.  The  careful 
experimental  studies  by  Janach  (1976)  and  Young  and  Powell  (1979)  on  rocks,  for 
example,  showed  that  the  observed  strain  rate  effects  are  most  likely  not  the  true  rate 
dependence  of  the  rocks  but.  instead,  are  due  to  strong  lateral  inertia  effects  arising 
from  dilatancy  (bulking)  during  failure. 

To  better  appreciate  the  complexities  associated  with  the  use  of  a  SHB  in 
conjunction  with  brittle  solids,  consider  the  nature  of  the  loading  and  its  subsequent 
effect  on  the  behavior  of  a  brittle  specimen.  Because  the  impedance  of  test 
specimens  is  usually  smaller  than  that  of  the  high  strength  steel  bars,  the  first  wave 
to  enter  the  specimen  usually  does  not  load  it  to  the  peak  stress.  Only  after  a 
number  of  reflections  of  the  wave  between  the  specimen-bar  interfaces  does  the  axial 
stress  reach  its  peak.  The  loading  time  is  proportional  to  the  specimen  length  and  to 
the  number  of  reflections  to  reach  the  peak  stress. 

If  the  specimen  is  a  brittle  material,  it  will  begin  to  fail  near  the  lateral  surface 
if  the  axial  stress  becomes  sufficiently  high,  even  though  it  has  not  reached  its  peak 
value.  Thus,  unlike  ductile  materials,  failure  in  brittle  materials  can  begin  near  the 
lateral  surface  while  the  axial  loading  process  is  still  in  progress.  Such  behavior  is.  of 
course,  highly  inhomogeneous. 


Cylindrical  Specimen 
of  Plain  Concrete 


Figure  4.12  The  standard  split  Hopkinson  bar  configuration. 
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Figure  4.13  Split  Hopkinson  bar  configuration  used  by  Young  and  Powell  to 
study  radial  inertia  effects  in  rocks. 


The  behavior  of  brittle  specimens  in  a  SHB  test  is  further  complicated  by 
dilatancy  (bulking)  during  failure,  which  does  not  occur  in  ductile  materials  During 
the  process  of  brittle  failure  in  compression,  cracks  develop  and  break  the  material 
into  subregions,  and  these  subregions,  consisting  of  still  intact  material,  must  separate 
in  order  that  the  radial  unloading  process  can  be  completed  The  unloading  process 
occurs  through  separations  of  these  subregions  of  cracked  material  and  is  accompanied 
by  frictional  sliding  and  the  formation  of  internal  voids,  which  produce  an  inelastic 
bulking  of  the  failing  material.  When  this  occurs,  radial  inertia  can  become  very 
important  because  the  reaction  forces  to  radial  acceleration  of  the  fractured  subregions 
produce  a  confining  pressure  on  the  not  yet  unloaded  interior  material.  Inasmuch  as 
the  strength  of  most  brittle  materials  increases  with  confining  pressure,  the  effective 
confining  pressure  produced  by  the  radial  acceleration  will  increase  the  strength  of  the 
interior  material  When  this  is  the  case,  it  becomes  difficult,  but  not  impossible  as 
discussed  below,  to  separate  the  true  rate  effect  from  the  apparent  rate  effect  due  to 
radial  inertia  effects. 

While  the  work  of  Janach  (1976)  was  the  first  to  point  out  the  importance  of 
considering  radial  inertia  effects  in  SHB  tests.  Young  and  Powell  (1979)  were  the  first 
to  actually  measure  them  Using  embedded  wire  loops  and  an  externally  applied 
longitudinal  magnetic  field,  as  depicted  in  Figure  4  13.  radial  velocity  histories  were 
measured  at  interior  points  of  rock  specimens.  From  these  histories,  the  radial 
accelerations  were  determined,  as  well  as  the  confinement  produced  by  the  radial 
stress.  Although  the  values  of  the  peak  axial  strain  rates  achieved  in  the  tests  were 
not  reported  by  Young  and  Powell,  it  is  highly  probable,  in  view  of  the  high  stresses 
generated  in  the  specimens  on  impact,  that  the  axial  strain  rates  may  have  been  as 
high  as  10*  to  105  sec**.  Young  and  Powell  concluded  from  this  study  that  the 
effects  of  radial  inertia  were  sufficiently  large  to  account  for  the  observed  increases  in 
compressive  strength.  They  cautioned  that  it  is  not  justifiable  to  attribute  the 
increased  strengths  of  the  rock  specimens  considered  to  a  strain  rate  effect. 

A  different  conclusion  was  reached  recently  by  Malvern  et  al  (1985).  who 

performed  a  series  of  SHB  tests  on  four  different  types  of  high  strength  plain 

concrete.  The  results  are  depicted  in  Figure  4.14.  where  the  variation  of  the  strength. 

o  with  strain  rate  is  shown  for  four  different  high  strength  concretes.  To 
max 
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accommodate  the  large  tests  specimens  required  for  testing  concrete,  an  unusually 
large  SHB  system  was  constructured  for  this  purpose.  In  an  effort  to  obtain  insight 
into  the  possible  effect  of  radial  inertia,  two  additional  strain  gauges  were  mounted  on 
the  specimen  surface  midway  between  its  ends.  One  gauge  measured  the  axial  strain 
while  the  other  measured  the  hoop  (e^)  strain.  By  numerically  differentiating  a 
smoothed  history  of  e^,  the  history  of  the  strain  rate  le^l  was  constructed  Upon 
comparing  this  with  the  history  of  the  axial  stress,  it  was  concluded  that  radial 
inertia  could  not  be  the  cause  of  the  increased  strengths  observed  in  the  tests. 
However,  there  were  some  features  of  the  data  that  could  not  be  explained,  which  led 
Malvern  et  a/..  to  comment  that.  "There  are  some  very  important  questions  involved 
in  these  dynamic  tests  for  which  we  do  not  yet  have  answers."  It  is  also  somewhat 
surprising  that  they  find  that  tests  conducted  on  longer  specimens  show  higher 
strengths  than  those  on  shorter  specimens  at  the  same  strain  rates.  A  number  of 
considerations  would  lead  one  to  expect  that  the  opposite  might  be  true. 

The  conclusion  reached  by  Janach  (1976)  and  Young  and  Powell  (1979) 
regarding  the  importance  of  radial  inertia  effects  in  SHB  tests  on  rocks  appears  to  be 
in  disagreement  with  the  conclusions  reached  by  Malvern  et  a/..  (1985)  from  SHB 
tests  on  plain  concrete.  These  differing  conclusions  may.  however,  not  be  inconsistent 
if  one  takes  into  account  the  possible  large  differences  in  the  strain  rates  between  the 
studies.  Malvern's  tests  were  done  at  strain  rates  in  the  range  from  10  to  102  sec  \ 
while  Young  and  Powell's  tests  appear  to  have  been  done  at  much  higher  rates, 
possibly  as  high  as  104  to  105  sec  .  It  is  to  be  expected  that  radial  inertia  effects 
will  grow  in  importance  as  the  strain  rate  increases.  Consequently,  it  may  be  that 
Malvern's  tests  were  done  at  sufficiently  low  strain  rates  that  radial  inertia  was  not 
important,  while  Young  and  Powell's  tests,  done  at  considerably  higher  rates,  were  in 
a  range  where  radial  inertia  dominated  the  response.  The  studies  by  Malvern,  et  al. 
(1985)  and  by  Young  and  Powell  (1979)  have  provided  valuable  information,  but  there 
is  still  considerable  uncertainty  and  unresolved  issues  in  the  way  brittle  specimens 
respond  in  a  SHB  test. 

Gas  Gun  Tests. 

Gas  guns  have  been  used  in  a  very  limited  manner  in  the  past  to  study  the 
dynamic  behavior  of  concrete.  In  fact,  only  two  such  studies  have  been  reported,  the 
first  of  which  was  conducted  by  Gregson  (1971)  and  a  more  recent  one  by  Gupta  and 
Seaman  (1979).  Of  these  two,  only  the  study  by  Gupta  and  Seaman  provides  direct 


insight  into  the  effect  of  strain  rate  on  the  behavior  of  plain  concrete  The  advantage 
that  the  gas  gun  has  over  the  drop-hammer  and  SHB  devices  is  that,  because  it 
subjects  specimens  to  uniaxial  strain  deformation  it  is  free  of  radial  inertia  effects. 
Its  disadvantage  at  least  as  far  as  the  present  strain  rate  range  of  interest  (10  1  to 
102  sec  *)  is  concerned,  is  that  it  typically  produces  strain  rates  in  concrete  that  are 
above  this  range.  Nonetheless,  gas  gun  data  on  concrete  at  the  higher  strain  rates 
can  serve  as  a  guide  in  modeling  the  effects  at  the  lower  rates,  since  a  suitable 
model  should  properly  extrapolate  to  the  higher  rates. 

In  the  gas  gun  study  by  Gupta  and  Seaman  (1979).  Lagrangian  gauge  analysis 
techniques  were  used  to  unfold  the  dynamic  stress-strain  relationship  from  the 
compressive  loading  data,  and  the  results  so  obtained  are  depicted  in  Figure  4  15(a). 
Here,  the  measured  dynamic  stress-strain  paths  for  two  compressive  loading  shots  are 
shown  together  with  a  corresponding  quasi-static  path  (also  obtained  under  uniaxial 
strain  conditions)  An  inspection  of  this  figure  reveals  that  the  dynamic  stress-strain 
curve  is  well  above  the  static  stress-strain  curve,  indicating  the  presence  of  a 
substantial  rate  effect.  The  corresponding  particle  velocity  histories  measured  at 
several  depths  in  the  high  stress  shot  (93.000  psi)  are  shown  in  Figure  4.15(b) 
From  this,  one  can  observe  the  rapid  attenuation  of  the  pulse  as  it  propagates 
through  the  specimen,  as  well  as  the  increasing  dispersion  of  the  wave  front. 

Using  the  data  from  the  high  stress  shot,  the  peak  stress  and  the  corresponding 
strain  rate  were  determined  for  each  of  the  four  wave  profiles  shown  in  Figure  4.15- 
(b)  The  results  are  shown  in  Figure  4.11  presented  earlier  as  the  solid  and  open 
triangles.  According  to  Seaman  (1984).  there  is  reasonable  confidence  in  the  data 
denoted  by  the  solid  triangles,  but  less  confidence  in  that  represented  by  the  open 
triangles.  Note  that  we  have  taken  the  liberty  to  plot  in  this  figure  the  data  from 
two  different  modes  of  deformation,  i.e..  drop-hammer  data  (uniaxial  stress)  and  gas 
gun  data  (uniaxial  strain).  The  validity  of  presenting  such  data  on  a  single  plot  has 
not  been  explored. 

4.6  EXPERIMENTAL  METHODS  AND  FAILURE  MODES. 

In  developing  constitutive  models  for  concrete  and  rock  on  the  basis  of  data 
from  laboratory  tests,  theoreticians  often  tacitly  assume  that  certain  fundamental 
premises,  such  as  homogeneous  fields  and  material  continuity,  are  satisfied  In 


many  cases,  appropriate  information  concerning  these  premises  is  not  supplied  by  the 
experimenters.  This  is  not  necessarily  the  fault  of  the  experimentalist,  since  he  may 
not  have  been  aware  that  certain  features  of  a  test  were  important  to  the  theoretician 
and  needed  to  be  measured  and  reported.  Also,  it  often  happens  that  experiments 
are  conducted  before  a  model  has  been  developed  and  the  data  needs  are  defined.  In 
most  all  cases,  the  interaction  between  theoreticians  and  experimentalists  appears  to 
be  insufficient,  and  has  led  to  a  propagation  of  misinformation  for  certain  aspects  of 
response.  One  such  case  involves  the  post-peak  response  of  concrete  and  rock  below 
the  brittle/ductile  transition  (Read  and  Hegemier  1984:  Hegemier  and  Read.  1984). 
Another  concerns  the  dependence  of  concrete  and  rock  properties  on  strain  rate,  as 
discussed  in  the  preceding  section. 

To  explore  these  issues  in  some  depth,  an  extensive  examination  of  experimental 
devices  and  test  data  normally  used  in  the  development  of  constitutive  models  for 
concrete  and  rock  was  conducted  as  part  of  the  present  study.  Various  issues 
pertaining  to  strain  softening,  strain  hardening,  failure  states,  failure  modes  and  strain 
rate  effects  were  considered.  Since  most  of  the  results  from  this  study  have  been 
described  in  detail  elsewhere  (Hegemier  and  Read.  1985).  we  shall  focus  here  primarily 
on  questions  pertaining  to  failure,  which  have  not  been  touched  on  before. 

First  of  all.  it  should  be  noted  that  there  is  no  unanimity  among  researchers  on 
the  definition  of  failure.  In  some  cases,  failure  is  defined  in  terms  of  the  maximum 
attainable  stress  state.  In  others,  failure  is  defined  by  the  onset  of  dilatancy.  In  still 
other  cases,  failure  is  defined  by  a  critical  state  of  strain.  Moreover,  the  term 
"failure"  in  the  literature  is  often  applied  to  brittle  failure  via  cleavage  or  shear 
faulting,  and  the  onset  of  ductile  flow.  The  current  presumption  is  that  differences 
between  definitions  of  failure  do  not  lead  to  significant  differences  in  failure  data. 

Comparison  of  data  from  tests  in  which  the  stress  paths  are  prescribed  reveals 
that  test  devices  which  employ  "unconstrained”  boundary  conditions  (of  a  type  which 
do  not  lead  to  early  failure  via  extrusion  of  lubricating  materials)  furnish  failure  data 
in  stress  space  that  is  reasonably  clustered.  In  contrast,  measured  strains  at  failure 
exhibit  vast  data  scatter.  It  is  suspected  that  part  of  this  phemonenon  may  be  due 
to  inherent  difficulties  in  measuring  strains  in  unconstrained  fluid-cushion  type  test 
systems.  However,  significant  data  scatter  is  also  observed  from  specimen  to  specimen 


within  the  same  test  system.  In  view  of  the  above,  it  would  appear  preferable  to 
describe  failure  in  stress  space  rather  than  strain  space.  This  does  not.  however, 
solve  the  problem,  since  deformation  characterization  is  equally  as  important  as  stress 
characterization  in  most  problems.  In  fact,  in  many  cases,  deformation  plays  the 
major  role. 

The  mode  of  failure  is  also  an  important  piece  of  information  to  the  theoretician. 
Unfortunately,  failure  modes  are  poorly  documented  in  the  literature.  While  much  effort 
has  been  devoted  to  exploring  failure  criteria  in.  for  example,  stress  space,  little 
correlation  has  been  established  between  points  on  a  failure  surface  and  corresponding 
failure  modes. 

The  paucity  of  data  on  failure  modes  is  compounded  by  the  fact  that  failure 
modes  are.  as  one  would  expect,  sensitive  to  the  experimental  test  setup  —  in 
particular,  test  boundary  conditions.  In  the  biaxial  tests  of  mortar  and  concrete 
conducted  by  Andenaes.  Gerstle  and  Ko  (1977).  the  specimens  in  the  fluid  cushion 
boundary  condition  series  failed  by  either  splitting  (tensile  splitting)  of  the  cube 
parallel  to  the  free  surface,  by  corner  or  edge  failure,  or  by  a  combination  of  both.  A 
different  mode  of  failure  was  observed  in  their  tests  conducted  with  steel  platens. 
Here,  the  failure  surfaces  formed  characteristic  shear  cones  at  the  free  surfaces,  with 
an  angle  varying  between  12  degrees  and  18  degrees  with  the  major  principal  stress 
direction. 

Information  concerning  failure  modes  in  triaxial  stress  states  is  furnished  by  Mills 
and  Zimmerman  (1970).  who  conducted  uniaxial,  biaxial  and  triaxial  experiments  on 
cubical  concrete  specimens,  using  lubricated  boundary  conditions.  Two  types  of  basic 
triaxial  tests  were  considered.  Type  I  tests  consisted  of  cases  in  which  the  largest 
principal  stress  is  the  major  principal  stress  and  the  minor  and  intermediate  principal 
stresses  are  smaller  and  equal.  Type  II  tests  were  conducted  in  which  the  major  and 
intermediate  principal  stresses  are  equal  and  larger  than  the  minor  principal  stress. 
The  failed  cubes  from  the  uniaxial  and  Type  I  triaxial  tests  were  found  to  be  cracked 
in  planes  parallel  to  the  load.  In  the  Type  II  tests,  the  cubes  were  found  to  be 
cracked  in  planes  parallel  to  the  plane  of  the  two  larger  loads.  The  cracks  in  the 
tests  with  intermediate  stress  conditions  indicated  splitting  in  the  direction  of  the 
smallest  principal  stress;  however,  the  orientation  of  the  crack  planes  with  respect  to 
the  plane  of  the  largest  and  intermediate  principal  stress  did  not  exhibit  predictable 
trends.  The  above  mentioned  failure  modes  are  illustrated  schematically  in  Figure 
4.16. 
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Summary  of  typical  failure  modes  for  concrete  specimens 


The  recent  experimental  study  by  Van  Mier  (1984),  however,  reveals  that  the 
above  picture  of  failure  modes  in  concrete  may  be  somewhat  oversimplified.  In  this 
study,  a  new  true  triaxial  device  with  servo-controls  on  all  three  axes  was  used  to 
study  the  responses  of  cubes  and  prisms  of  a  medium  strength  plain  concrete  to 
uniaxial  and  triaxial  compression.  The  loads  were  applied  to  the  test  specimens 
through  brush-bearing  platens  which  deform  them  with  reduced  friction. 

Uniaxial  compression  tests  were  performed  on  prisms  of  different  slenderness- 
ratios  (h/d  =  0.5.  1.0  and  2.00).  For  the  prisms  with  h/d  =  0.5  and  1.0.  failure 
occurred  through  fracture  of  the  specimen  into  cone-shaped  pieces  in  alternating 
positions,  where  the  heights  of  the  cones  were  restricted  by  the  height  of  the 
specimen.  In  the  case  of  the  prisms  with  h/d  =  2.0,  eight  tests  were  done  and.  of 
these,  four  failed  by  localized  shear  fracture,  while  the  other  four  failed  by  formation 
of  a  multitude  of  axial  cleavage  cracks.  Thus,  even  when  specimens  of  the  same 
concrete  mix  and  slenderness  ratio  were  tested  in  the  same  device,  the  fracture  mode 
was  not  unique.  The  cracking  patterns  for  several  of  the  uniaxial  compression  tests 
are  shown  in  Figure  4.17. 

Triaxial  tests  were  done  only  on  cubical  specimens  When  the  applied  loading 
was  such  that  extension  occurred  in  only  one  direction  (planar  mode),  the  specimen 
failed  through  the  development  of  localized  shear  bands.  On  the  other  hand,  when 
loading  caused  extension  in  two  directions  (cylindrical  mode),  failure  occurred  by 
distributed  fracture  throughout  the  sample.  Typical  cracking  patterns  for  the  triaxial 
tests  are  also  depicted  in  Figure  4.17. 

The  results  from  van  Mier’s  study  dramatically  underscore  the  complex  nature  of 
failure  in  plain  concrete.  Inasmuch  as  van  Mier  took  special  care  to  reduce  to  the 
extent  possible  friction  between  the  test  specimens  and  the  loading  device  by  using 
brush-bearing  platens,  the  observed  modes  of  failure  are  believed  to  reflect  true 
material  behavior  and  not  system  effects.  From  these  results,  we  see  that  the  mode 
of  failure  depends  on  both  a  specimen's  slenderness-ratio  and  the  relative  magnitudes 
of  the  applied  principal  stresses. 

In  conclusion,  it  is  clear  that  failure  modes  are  strongly  influenced  by  a  number 
of  factors,  including  test  boundary  conditions.  When  unconstrained  boundary 
conditions  are  used,  similar  general  trends  in  failure  modes  are  observed.  Nonetheless. 


quantitatively  there  is  considerable  data  scatter  with  respect  to  failure  modes, 
especially  at  low  confining  pressures.  The  above  problems  are  further  compounded  by 
generally  poor  descriptions  of  the  failure  mode,  failure  sequence,  and  state  of  the 
specimen  during  and  after  failure.  For  further  elaboration  on  the  topics  discussed  in 
this  section,  see  the  paper  by  Hegemier  and  Read  (1985). 
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APPENDIX  A 


NEW  CONTINUOUS  DAMAGE  MODEL  FOR  CONCRETE  1 

I 

i 

A.l  INTRODUCTION-  j 

! 

The  purpose  of  this  appendix  is  to  describe,  in  detail,  a  damage-cracking  model 
that  has  been  developed  under  the  present  program.  The  model  was  formulated  by 
Valanis  (1985)  and  is  presently  under  development.  The  discussion  given  below 
largely  follows  that  given  by  Valanis  (1985)  but  does  not  reflect  various  topics  that 
are  presently  under  study,  which  include  (a)  the  addition  of  plasticity,  (b)  size  effects 
and  (c)  strain  rate  dependence.  | 

A. 2  FORMULATION  OF  MODEL.  ! 

I 

I 

Consider  a  linearly  elastic  material  which  is  isotropic  in  its  virgin  unstrained 
state  and  undergoes  small  deformation  at  isothermal  conditions.  As  the  deformation 
proceeds,  damage  gradually  develops,  which  in  turn  reduces  the  material  integrity.  For  | 

sufficiently  large  deformation,  the  accumulated  damage  can  lead  to  complete  fracture.  1 

in  which  the  material  cannot  support  tensile  stress  in  a  direction  normal  to  the 
damage.  We  shall  formulate  this  damage-cracking  process  within  the  context  of  ' 

irreversible  thermodynamics  and,  for  this  purpose,  we  introduce  the  free  energy  ♦  per 
unit  volume.  Since  fracture  leads  to  a  reduction  in  material  integrity,  we  introduce  an 
integrity  tensor  g  which  will  be  discussed  at  length  later  and  which  is  symmetrical 

such  that  g  =  g  (the  unit  tensor),  when  the  material  is  in  its  virgin  (undamaged) 

state,  and  g  =  Q  when  the  material  has  fully  failed,  i.e..  when  the  material  cannot 
support  stress  in  any  direction.  We  thus  set: 


*  =  *(£,*)  (A-l) 

where  £  denotes  the  strain  tensor.  In  a  thermodynamic  sense,  g  now  plays  the  role 
of  an  internal  variable,  in  which  case  we  can  write  the  following  relations: 


2  = 


9£ 


(A-2) 


*  This  nomenclature  is  preferred  to  free  energy  density,  in  view  of  the  concept  of 
"irreducible"  material  volume  which  must  be  introduced  when  dealing  with 
heterogeneous  multi-phase  materials,  such  as  plain  concrete. 


where  g  is  the  stress  tensor  and  g  denotes  the  internal  force,  dual  to  £  which  is 
driving  the  fracturing  process. 


In  view  of  the  definition  of  g.  the  following  relations  must  hold: 


*(£,2)  =  0  ,  t(2,4)  =  o  ,  (A— 3) 

the  first  meaning  that  a  fully  failed  material  cannot  contain  free  energy  and  the 
second  that  an  unstrained  material  must  have  zero  free  energy;  both  conditions  are 
relative  to  the  reference  state.  Furthermore,  the  following  conditions  must  also  hold: 


Z(£,Q)  =  2  *  Z(2  >i)  =  2  (A~4) 

8U*S)  =  2  (A-5) 


Equations  (A-4)  require  that  the  stress  vanish  in  the  fully  failed  material  and  that, 
since  the  material  remains  elastic  in  the  damaged  state,  the  stress  will  vanish  at  zero 
strain.  Equation  (3.5)  stipulates  that  the  internal  fracture  causing  force  g  must 
vanish  in  the  fully  failed  material. 

If  we  now  expand  *  in  a  Taylor  series  in  g.  retain  terms  no  higher  than  the 
quadratic,  and  observe  relations  (3.3).  we  obtain  the  following  form  of  *: 


i 


♦  =  C 


i  j  k£.C 


ijek£ 


(A-6) 


where  =  £(£).  It  is  further  assumed  that  £  depends  upon  £  in  a  purely  quadratic 
manner.  Since  the  material  is  assumed  to  be  isotropic  in  its  virgin  state.  £  may  be 
represented  in  terms  of  outer  products  of  the  unit  tensor  g  and  Furthermore, 
since  £  is  purely  quadratic  in  g.  it  has  no  other  representation  than 


C. 


i  j  k£. 


=  X 


r'Jr 


k  l 


♦  2/* 


ikrj£' 


(A— 7) 


in  view  of  the  symmetries  imposed  upon  it  by  the  symmetry  of  g  implied  in  Eq.  (A- 
6).  The  constants  X  and  2 n  in  Eq.  (A-7)  must  indeed  by  the  Lame  constants  of  the 
virgin  material,  since 


*  This  assumption  has  recently  been  verified  by  Valanis  (1987).  using  a  non¬ 
thermodynamic  approach  involving  a  mapping  from  undamaged  to  damaged 
material. 


when  £  =  g.  Thus,  we  can  write 


*  =  2  X(*ij£ij)2  +  *  *ik*jteijCk£  ' 
so  that,  in  view  of  Eq.  (A-2a),  g  us  given  by  the  expression: 

aij  =X'i/k£ek£*2'‘'ikW 


(A-9) 


(A-10) 


This  equation  defines  the  constitutive  response  of  an  elastic  fracturing  material,  once 
we  establish  how  the  integrity  tensor  g  is  related  to  the  fracturing  process.  In 
particular,  an  equation  is  needed  which  describes  the  manner  in  which  g  evolves  with 
accumulated  damage,  and  this  is  considered  below. 

A. 3  EVOLUTION  EQUATION  FOR  ^(DIRECTED  DAMAGE)- 

When  the  material  is  not  strain  rate  sensitive,  an  increment  of  tensile  strain  will 

produce  an  increment  of  damage.  Specifically,  given  an  increment  in  strain,  dg.  let 

dea  be  its  eigenvalues  and  n.fl  its  eigenvectors.  If  dg  and  g  are  coaxial,  then  de° 

will  constitute  an  increase  in  tensile  strain  if  de°  >  0  and  eB>  0.  where  e®  denote 

the  eigenvalues  of  g.  In  the  more  general  case  where  dg  and  |  aje  not  co-axial.  dea 

is  said  to  be  an  increase  in  tensile  strain  if  deB  >  0  and  e..n.  n.  =  e  >  0  (a  not 

ij  i  j  a  — . 

summed).  Now  let  djg  be  the  change  in  the  integrity  tensor  g  due  to  the  increment 
dg.  and  let  jjJ  be  the  eigenvectors  of  dg.  We  now  specify  that  d^  and  dg  are  co¬ 
axial.  so  thatr 


a 

D, 


(A— 11) 


The  physical  significance  of  Eq.  ( A- 11)  is  examined  in  the  discussion  following  Eq. 
(A- 14) 


To  completely  characterize  the  evolution  of  the  relation  between  d^a  and  dea 
needs  to  be  specified.  In  this  section,  we  consider  a  type  of  damage  that  is 
completely  "directed",  that  is.  it  develops  in  a  specific  direction  without  being 
accompanied  by  non-directed  (or  isotropic)  damage.  This  type  of  damage  is  the 
variety  dealt  with  in  linear  elastic  fracture  mechanics  (LEFM)  and  appears  to  typify 
that  which  occurs  in  single  phase  brittle  materials.  Later  on  in  Section  A. 7.  we  shall 
consider  the  more  general  case  in  which  the  advancing  damage  consists  of  both 


directed  and  non-directed  (isotropic)  components.  Such  damage  is  characteristic  of 
multi-phase  materials,  such  as  plain  concrete,  as  discussed  by  Diamond  and  Bentur 
(1984). 


Restricting  attention  now  to  completely  "directed”  damage,  we  adopt  the 
following  expression  relating  d^a  and  dea: 

<v*  =  -  [K  ]  V  (*-«) 


where 


kdea  ,  for  dea  >  0  and  ea  2  0 
.  0  ,  otherwise 


(A-13) 


Here,  k  and  m  are  positive  constants,  and  4a  =  4..nfln°  la  not  summed).  The  form 
of  Eq.  (A-12)  has  its  basis  in  the  physics  of  annihilation,  but  may  be  changed  to 
provide  a  more  realistic  description  of  damage  accumulation  for  a  particular  material. 


In  view  of  Eq.  (A-ll).  we  can  therefore  write 


^ij  *  -  £  to)’-?"*  if.  (»->«) 

a 

With  this  evolution  equation  for  the  constitutive  description  of  the  elastic  fracturing 
solid  is  now  complete,  once  the  material  parameters  X.  /i.  k  and  m  have  been 
specified. 

Physical  Interpretation  of  Equations 

The  physical  meaning  of  the  above  equations,  and  their  relation  to  the  fracturing 
process,  is  best  understood  by  reference  to  Eq.  (A-10).  which  relates  the  stress  tensor 
to  the  strain  tensor.  Let  f/1  be  the  eigenvectors  of  ^  and  let  g  and  £  be  referred  to 
a  system  of  coordinates  Xfl.  Then,  the  components  of  g  and  £  on  the  Xfl  coordinate 
system  can  be  written  as: 


In  this  system  of  coordinates,  g  is  diagonal,  i.e..  =  0  for  r  i  s.  and  Eq.  (A- 10) 

takes  the  form: 


a  =  X  4  4  e  *  2u  4  4  c 

rs  rrsr»n  mn  r  rrmrsn  mn 


(A-17) 


Since  £  is  diagonal,  we  can  write: 

_  3 

ff,.  =  X  iL,  Vi  e  ♦ 

11  rll  ^-rrrr  rr  rll  11 
r=l 


a21  =  °12  ~  2*i  ^11^22C12 


(A-18) 


^31  =  ^13  =  2/i  ^11^22C13 

Let  us  now  single  out  the  eigenvector  N1  and  examine  what  happens  to  the  stress 
when  the  eigenvalue  ^(=  vanishes.  Note  from  Eqs.  (A-18)  that,  if  ^  =  0.  it 

follows  that  =  <?i2  =  a  =  0.  that  is.  the  solid  cannot  support  stresses  on  a 
plane  normal  to  N1,  on  which  5  ^  =  0.  The  physical  meaning  of  this  is  that 
the  decrease  of  represents  damage  normal  to  the  plane  of  N1  and.  as  such,  it  is  a 
measure  of  the  plane  microcracks  that  have  developed  normal  to  N1  and/or  the 
increase  in  the  size  of  such  cracks,  so  that  when  =  0.  a  plane  crack  in  the 
accepted  sense  has  formed  across  the  material  element  on  a  plane  with  normal  N1.  so 
that  the  element  cannot  support  tensile  or  shear  stresses  on  that  plane. 

On  the  basis  of  the  above  observations,  the  physical  meaning  of  Eqs.  (A-ll). 
(A-12)  and  (A-13)  now  becomes  clear.  An  increase  in  tensile  strain  dea  on  a  plane 
with  normal  na  causes  planar  damage  d^a  on  that  plane,  where  d^a  is  given  by  Eq. 
(A-12).  in  accordance  with  the  observation  that  "planar  microcracks  form  perpendicular 
to  the  direction  of  the  principal  tensile  strain." 


The  constant  k  in  Eq.  (A-13)  is  a  material  parameter  that  reflects  the  fracture 
resistance  or  fracture  toughness  of  the  material.  Eq.  (A-12)  is  representative  of 
processes  involving  annihilation  of  populations  and  catastrophic  processes  in  systems 
where  the  increment  of  annihilation  is  proportional  to  the  state  of  integrity,  which  is 
given  here  by  the  factor  (^)m  This  form  is  representative  of  such  systems,  even 
though  any  monotonically  decreasing  function  f(#;  |).  such  that  f(0)  =  0.  will  likely  do 
as  well  Finally,  we  note  that  the  damage  tensor  Q  can  be  expressed  in  terms  of  £ 
as  follows: 


The  physical  interpretation  of  the  governing  equations  is  now  complete. 


A. 4  APPLICATION  TO  SOME  SIMPLE  CASES  OF  HOMOGENEOUS 
DEFORMATION. 

In  this  section,  analytic  solutions,  based  upon  the  constitutive  model  described 
above,  are  presented  for  several  simple  cases  of  homogeneous  deformation  including 
simple  tension,  simple  compression  and  simple  shear.  These  solutions  provide  valuable 
insight  into  the  characteristics  of  the  model  and  reveal  its  remarkable  predictive 
capabilities,  despite  the  fact  that  it  involves  only  four  material  parameters. 

A. 4.1  Homogeneous  Triaxial  Strain  Fields. 

As  a  preliminary  development,  consider  the  case  in  which  the  strain  field  is 
homogeneous  and  consists  only  of  principal  strains,  i.e..  e..  =  0  for  i  t  j.  Thus.  d£ 
and  g  are  always  coaxial  and  g  is  coaxial  with  £  by  virtue  of  Eq.  (A  ll).  As  a 
result,  g  is  diagonal.  In  view  of  Eq.  (A-10).  g  is  also  diagonal.  The  principal 
stresses  are  given  below  in  terms  of  the  principal  strains  and  the  principal  values  of 

t 


ol  =  X  ^rer  ♦ 

°2  ~  X  *2  21  *r€r  *  2/i  *2e2  (*-20) 

=  X  ^2-  *rer  *  2^  *3C3 
r 

The  appropriate  evolution  equation  for  g  is  obtained  from  Eq.  (A-14)  and  is  given  in 
terms  of  the  principal  components  of  g  as  follows: 


where 


(k  de.  ,  if  de-  >  0  , 

<Kj  -  ' 

L  0  ,  otherwise 

We  now  consider  two  special  cases 
and  simple  compression. 

A.4.2  Simple  Tension. 

Consider  the  case  of  monotonically  increasing  simple  tension.  In  this  case,  c  — 
a3  =  0.  and  integration  of  Eqs.  (A-21)  with  =  e3  and  the  initial  conditions  ^  = 
^3  =  1.  shows  that  ^  always.  Hence,  we  find  that 

^2C2  =  "  2(X^0  e\h~  "  ^1C1  ^A‘22^ 


£i^° 

of  the  above  equations,  namely,  simple  tension 


and 


a 


1 


where 


c  _  (2ti+3\)u 


(A-23) 


(A-24) 


From  Eq.  (A-21a).  the  evolution  equation  for  ^  is 


=  0 


subject  to  the  initial  condition  ^(0)  =  1.  Thus,  the  solution  is 


*1 


-kel 

e  ,  if  m  =  1 

1 

jl  -  (l-mJkCjj^  m  ,  if  m  i  1 


(A-25) 


(A-26) 


Combining  Eqs.  (A-23)  and  (A-26).  we  obtain  the  expression: 


-  2ke, 


,  if  m  =  1 


ffl  = 


E  e^|l-(l-m)ke^j^"m  ,  if  m  i  1 


(A-27) 


Inspection  of  Eq.  (A-27)  reveals  that,  for  0  <  m  <  1,  a ^  is  not  a  monotonic  function 
of  but  reaches  a  maximum  and  then  goes  to  zero  for  finite  e^.  For  1  m  <  3. 
Oj  again  reaches  a  maximum  and  then  decays  monotonically  to  zero  as  ■»  «*.  For 
m  =  3.  a  monotonically  increases  to  an  asymptotic  constant  value  as  e 3  *  *.  while 
for  3  <  m  <  *.  a  increases  monotonically  with  The  influence  of  m  on  the 

constitutive  response  is  shown  in  Figure  A. la.  Thus,  in  the  range  0  m  <  3.  the 
model  exhibits  softening  in  tension  in  agreement  with  experimental  data. 

2 

The  model  also  unloads  elasticity  with  an  effective  modulus  of  E^j  which  is 
diminished  by  the  accumulated  damage.  Thus,  the  model  has  all  of  the 
characteristics  of  a  fracturing  elastic  solid.  It  follows  from  Eq.  (A-22)  that  <  0 
and  hence  ^  =  1.  i.e.,  no  damage  is  produced  in  the  transverse  direction  due  to  the 

uniaxial  tension  loading. 

A.4.3  Simple  Compression. 


Consider  now  the  case  of  monotonically  increasing  uniaxial  compression.  Here. 
a2  =  a3  =  0.  while  and  are  compressive.  Again,  Eqs.  (A-21)  together  with  e2 

=  e3<  and  the  initial  conditions  f  f  =  1.  lead  to  the  result  that  ^  ^ 

always.  Similar  to  the  case  of  simple  tension,  we  find  that 


= 


(A-28) 


(A-29) 

The  damage  evolution  in  this  case  is,  however,  entirely  different  from  that  for  simple 
tension.  Since  both  de 3  and  6^  are  compressive  at  all  times,  we  find  from  Eq.  (A- 
21)  that  d^  =  0  and  thus  ^  =  1  always.  However,  in  view  of  Eq.  (A-28).  e2  is 
now  tensile  and  hence  (and  03)  will  increase  with  the  consequence  that  damage 
will  develop  on  planes  wnich  are  parallel  to  the  axis  of  compression,  in  accordance 
with  experimental  observation:  this  is  the  so-called  axial  splitting  mode  (Horii  and 
Nemat-Nasser,  1985). 
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(a)  Stress-strain  relation  in  simple  tension 
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Figure  A.l 


(b)  Axial  strain  vs.  lateral  strain  in  simple 
compression 


Predicted  responses  for  simple  tension  and  simple  compression,  showing  the 
effect  of  the  parameters  m  on  the  resulting  behavior. 


».V  W.W7»T7*T^’ 


From  Eq.  (A-21).  we  can  write 

d^2  +  de2  =  0  '  (A-30) 

the  solution  of  which  is 


[  c"ke2  ,  m  =  1 

*2  =  j  _±_  (A-31) 

[  |l  -  (l-m)ke2j1-m  ,  m  t  1 

Thus,  in  view  of  Eq.  (A-29)  and  the  condition  ^  =  1.  the  axial  stress  strain  relation 
is 

=  E  £j  ,  (A-32) 

so  that  the  response  in  the  axial  direction  is  purely  elastic. 

To  find  the  relation  between  the  strains,  we  use  Eqs.  (A-28)  and  (A-31) 
together  with  the  fact  that  ^  =  1.  It  then  follows  that 


-  uex  = 


-ke. 


e2e 


,  m  =  1 


(A-33) 


[e2|l  -  (l-m)ke2j1-,n  ,  m  t  1 
Here,  we  may  distinguish  three  different  cases: 

(a)  0  <  m  <  1 

The  strain  e ^  is  not  a  monotonic  function  of  It  reaches  a  maximum 

absolute  value  and  then  goes  to  zero  for  a  finite  value  of  Specifically. 


l£lLax  =  Vv  (n+l) 


n+1 


(A-34) 


where  n  =  1/ (1-m) . 


Again,  the  strain  is  not  a  monotonic  function  of  It  reaches  a  maximum 
absolute  value  again  given  by  Eq.  (A-34)  for  1  <  m  and  then  goes  to  zero  for 
infinite  e2_  For  m  =  1 


l£lJmax  “  i/ek  ' 

where  e  denotes  the  base  of  the  natural  logarithm. 


(A-35) 


(c)  2  <  m  <  • 

In  this  case.  Cj  is  a  monotonic  function  of  and  approaches  the  following 
limiting  condition  for  increasing  lateral  strain: 


I  im  I  I 
|£l|max 


\l 


l  « 


m  =  2 
m  >  2 


The  above  cases  are  depicted  graphically  in  Figure  A.  1(b). 

A.4.4  Collapse  of  a  Block  Under  Axial  Compression. 


(A-36) 


It  follows  from  the  above  discussion  that  a  block  consisting  of  elastic  fracturing 
material  will  collapse  under  homogeneous  axial  compression  if  0  <  m  <  2.  since  e1 
has  an  upper  bound  given  by  the  expression: 


Umax 


2-m 

-  .1  f  1  I 1-m 

vV.  12-mJ 


(A-37) 


The  limiting  cases  are 

,  m  =  1 

lit  >  ■  * 2 

The  collapse  stress  may  be  calculated  directly  from  Eq.  (A-32),  with  the  result. 


'1 1  max 


Q_ 

i/ke 
1 


(A-38) 


2-m 

i'll*.  *  s  fey1"1  <a-39> 

Thus,  the  theory  predicts  the  collapse  of  a  block  under  axial  compression  due  to 
damage  on  planes  parallel  to  the  axis  of  compression.  Again,  this  is  the  axial 
splitting  mode. 
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A.4.5  Ratio  of  Collapse  Stresses  in  Tension  and  Compression. 

At  this  point,  it  is  of  interest  to  explore  the  rationality  of  the  above  results  by 
comparing  the  ratio  ac/al.  where  crc  and  a1  are  collapse  stresses  in  compression  and 
tension,  respectively,  obtained  analytically,  with  that  observed  experimentally  for 
concrete,  even  though  the  end-grip  conditions  in  the  tests  may  not  be  ideal,  as 
assumed  in  solutions  obtained  here.  Restricting  attention  to  the  case  m  =  1.  it 
follows  from  Eqs.  (A-37)  and  (A-39)  that 


(A-40) 


Since  u  ~  0.2  for  concrete,  we  have 


-  10  , 
T.  ~  9 


(A— 41) 


which  is  close  to  experimental  observation  (Raphael.  1984). 

A.4.6  Simple  Shear. 

Consider  now  the  case  of  monotonically  increasing  shear  strain  under  conditions 
of  simple  shear  and  small  strain.  The  configuration  of  interest  is  depicted  in  Figure 
A. 2.  Inasmuch  as  the  only  non-zero  strain  component  is  we  can  write 


k£.ck£  =  2  ^12e12 


^ik^jCek£=  i  l^j  2  +  ^i2^jl)el2 


Therefore,  in  this  case,  the  general  constitutive  relation 


<v,  =  x  *ij<Wkc  +  ^  *ik*iPe 


(A-42) 


(A-43) 


ik*j£.ck£ 


(A-44) 


reduces  to  the  form 


aii  =  2[X  ^i j^l2  +  +  ^i2^jl) ] e] 


(A-45) 


Therefore,  we  can  write 
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al2  =  2[X  *12  *  ^(*11*22  *  ^12]  ]e: 
*11  =  2(X  *  ^0*ll*12e12 
*22  ~  2[X  *  2^)*22*12e12 


(A-46) 

(A-47) 


(A-48) 


For  a  small  shear  strain,  the  maximum  tensile  strain  will  lie  in  a  direction  that 
makes  a  45  degree  angle  with  each  of  the  coordinate  axes.  Consequently,  cracks  will 
develop  in  the  material  as  shown  in  Figure  A. 2.  The  components  of  the  unit  normal 
Q  to  the  crack  are  nj  =  n2  =  l/JT 

Since  there  is  only  one  cracking  direction,  we  set  a  =  1  in  the  evolutionary 
equation  for  d^..  and  write 


d*u  *  -  k('-»)n,de  "i"j 


(A-49) 


where 


4  =  4.  .n.n. 

rn  rij  1  j 


de  =  de. 


(A-50) 


Since  n^  =  -  1/J2,  it  follows  from  Eq.  (A-49)  that 

d*ll  =  d*22  =  d*12 


(A-51) 


Hence 


-rn  -  r\2  • 

1  J0 


since  4..  -  5..  initially. 
r'J  'J 

*11  =  1  *  * 12  ' 

and  Eq.  (A-50a)  gives 
*n  s  *11  *  *12  ' 


(A-52) 


From  Eq.  (A-53)  we  find 


(A-53) 


(A-54) 
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since  =  ^2  according  to  Eq.  (A-51)  . 


Substitution  of  Eq.  (A-53)  into  Eq.  (A-54)  yields 

K  *  1  *  2  #12 

which,  when  combined  with  Eq.  (A-49) ,  leads  to  the  result 

d*ij  *  -  kP  *  2^12]"  d£12"i"j 

Therefore,  since  n^  =  =  l/f2,  we  can  write 

d^12  =  _  2  t1  *  2  ^12) m  de12 

which  can  be  integrated  to  give,  for  m  =  1: 


(A-55) 


(A-56) 


(A-57) 


,  1  "kc12  1 

^12  "  2  le  " 


(A-58) 


and  for  m  ^  1: 


a  I  1  m-1 

f12  "  2  |[k(m-l)e12«-lj  “  1 

Recalling  Eq.  (A-53)  and  the  fact  that  022  =  Eq.  (A-46)  can  be  placed  in  the 
following  form: 

°12  1  2[<X  *  Mth  *  2^12  *  £12 

Therefore.  Eqs.  (A-58)  to  (A-60)  describe  the  relationship  between  aJ2  and  e 
From  Eqs.  (A-58)  and  (A-59)  note  that 


(A-59) 


(A-60) 


1 '»  aio  =  "  o  (x  +  2/t) e 


(A— 61) 
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which  in  view  of  Eq.  (A-60)  leads  to  the  result: 


1  im  ffio  =  o  (x  ♦  2^) e - , 


(A-62) 


’t  », 


I*  ’ 

k 


iV.. 


*<  *  ^  -.i  ' 


Hence,  when  X  <  2/*.  the  ultimate  slope  given  by  Eq.  (A-62)  will  be  less  than  the 
initial  slope.  2 p.  as  shown  in  Figure  A. 3.  The  case  X  >  2fi  on  the  other  hand,  leads 
to  the  ultimate  slope  being  greater  than  the  initial  slope,  i.e..  a  hardening,  which  is 
difficult  to  visualize  from  a  physical  standpoint.  These  limiting  cases,  however,  fall 
outside  of  the  assumption  of  small  strains  made  here  and  therefore  may  not  be 
physically  meaningful. 

A. 5  APPLICATION  TO  NON-HOMOGENEOUS  STRAIN  FIELDS. 

In  an  effort  to  explore  the  capability  of  the  continuous  damage  model  to  predict 
the  proper  trends  for  fracture  problems  involving  non-homogeneous  strain  fields,  we 
developed  a  computer  program  for  the  model  and  incorporated  it  into  a  finite  element 
code.  We  then  conducted  a  limited  numerical  study  of  a  fracture  problem,  having  a 
non-homogeneous  strain  field,  for  which  a  theoretical  solution  is  available  from  linear 
elastic  fracture  mechanics  (LEFM).  The  problem,  depicted  in  Figure  A. 4(a),  consists 
of  a  long  flat  plate  of  width  w.  which  has  a  crack  of  length  2a  and  is  acted  upon  by 
a  uniform  tensile  stress  a.  The  problem  is  to  find  the  magnitude  of  a  which  causes 
the  crack  to  grow  catastrophically  and  thereby  cause  fracture  of  the  plate. 

The  finite  element  grid  used  in  the  numerical  study  consisted  of  uniform,  initially 
square  zones  of  unit  length.  Because  of  symmetry,  only  the  right  half  of  the  problem 
was  considered.  Vertically,  there  were  21  zones  while  16  zones  were  used  in  the 
horizontal  direction.  Therefore,  the  width  w  had  the  value  32.  The  effect  of  an 
initial  crack  in  the  x-direction  was  introduced  into  the  grid  by  simple  setting  ^  =0 

in  those  zones  through  which  the  crack  passed.  In  the  study.,  three  different  values 
of  initial  crack  length  were  considered,  namely,  a  =  1.  3  and  6. 

The  results  from  the  numerical  study  are  shown  in  Figure  A. 4(b).  where  the 
stress  af  to  fracture  the  plate  is  plotted  against  the  dimensionless  ratio  a/w.  Also 
shown  on  this  figure  is  the  corresponding  theoretical  prediction  from  LEFM  (see 
Knott.  1973).  which  is  given  by  the  expression: 


w  tan 


(")] 


(A. 63) 


where  Kj(..  the  mode  I  fracture  toughness,  was  set  to  the  value  7.07  so  that  the 
theoretical  curve  would  pass  through  the  central  numerical  point  shown  in  the  figure. 
As  an  inspection  of  the  figure  reveals,  the  continuous  damage  model  predicts  a  decay 


(a)  Geometric  configuration  of  flat  plate. 
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Figure  A. 4. 
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(b) Dependence  of  fracture  stress  af  on 
the  ratio  a/w. 


Response  of  a  thin  elastic  plate  with  a  crack  to  an  applied 
uniform  tensile  stress. 
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of  af  with  increasing  a/w.  which  is  qualitatively  in  agreement  with  the  theoretical 
result  from  LEFM. 

Exact  agreement  between  the  model  prediction  and  LEFM  is  not  to  be  expected, 
however,  for  several  reasons.  First,  the  model  produces  some  damage  in  every  zone 
of  the  grid  from  the  onset  of  the  tensile  loading,  while  in  LEFM  theory  it  is  assumed 
that  the  damage  remains  confined  to  the  original  crack  until  there  is  catastrophic 
crack  extension.  Secondly,  for  the  material  region  just  ahead  of  the  crack,  the  two 
approaches  lead  to  significantly  different  stress  fields.  In  LEFM.  the  material  ahead  of 
the  crack  remains  linearly  elastic  and  undamaged  until  overall  failure  occurs,  while,  in 
the  numerical  simulation  with  the  model,  the  zones  ahead  of  the  crack  experience 
increasingly  greater  damage  as  a  grows.  Thus,  the  stress  fields  in  the  material  just 
ahead  of  the  crack  are  very  much  different  in  the  two  cases.  In  view  of  these 
considerations,  the  differences  between  the  two  predictions  shown  in  Figure  A. 4  are 
not  surprising.  The  general  trends,  however,  are  similar  and  indicate  that  the  present 
model,  when  used  in  conjunction  with  a  finite  element  code,  leads  to  realistic  results. 
Closer  agreement  between  the  two  predictions  could  probably  be  obtained  by  revising 
the  evolution  equation  for  £. 

A. 6  POST-FAILURE  BEHAVIOR  AND  CRACK  CLOSURE. 

It  was  shown  in  Section  A. 3  that  if  an  eigenvalue  £A  of  £  vanishes  in  a 
material  element  then  the  element  cannot  support  stress  on  a  plane  normal  to  NA. 
the  latter  being  the  eigenvector  of  £  associated  with  £A.  In  a  structural  sense  the 
element  has  developed  damage  equivalent  to  a  crack  on  the  plane  normal  to  NA.  It 
follows  that  the  element  has  developed  kinematic  freedom  in  directions  normal  and 
tangential  fo  NA  in  the  following  sense.  Let  ^eap  be  the  components  of  the  strain 
increments  A£  defined  relative  to  a  system  of  coordinates  x“  coincident  with  the 
eigenvectors  Na  of  £.  Thus 
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Let  A  =  1.  Perusal  of  Eqs.  (A-18a.b.c)  show  that  if  l  ^  =  0.  then  not  only  ! 

are  =  a y  =  a 31  =  <x2J  =  0  as  discussed  previously  but  the  strains  e  j 

e12  =  e2j.  =  eJ3  become  indeterminate  from  the  constitutive  properties  of  tne  , 

element.  \ 


A. 7  DIRECTED  AND  ISOTROPIC  DAMAGE. 


It  has  recently  been  observed  (Diamond  and  Bentur.  1984)  that,  in  multi-phase 
materials,  such  as  plain  concrete,  the  formation  of  directional  cracks  is  usually 
accompanied  by  the  presence  of  small  cracks  randomly  oriented,  which  from  a 
statistical  viewpoint  amount  to  isotropic  damage,  denoted  here  by  Q.  In  this  sense  Q 
reduces  the  stress  carrying  capacity  of  a  material  by  the  same  degree  in  all  directions. 
Thus  the  integrity  tensor,  now  denoted  by  must  be  reduced  by  the  presence  of  Q 
and  hence  it  must  be  a  state  function  of  Q.  Thus,  in  symbolic  form: 


ia  *  <Q(“.y  (*-«) 

subject  to  the  following  conditions  which  must  always  apply: 

*Q(0,0)  =  fi 

(A-66) 

=2 


where  £  is  the  unit  matrix. 

To  satisfy  conditions  (A-66)  we  adopt  a  relation  of  the  multiplicative  type.  i.e..  we 
set 

<0  =  w(Q)*(£a)  (A-67) 

such  that  w(0)  =  1.  w(l)  =  0.  Furthermore  u  must  be  a  monotonically  decreasing 
function  of  Q.  A  suitable  function  that  satisfies  the  above  constraints  is 

w(G)  =  (1  -  G)a  (A-68) 


where  a  l  0.  Thus  the  integrity  tensor  jg.  considered  in  the  previous  sections,  is  now 
modified  by  the  presence  of  isotropic  damage  as  specified  by  Eqs.  (A-67)  and  (A-68) 


The  constitutive  equation  thus  becomes 


*  =  I  **  *ik^jtL£ijek<J 


(A-69) 


and 


aij  =  u2{x*]jKi£ki  *  2fi  ^ik^jteke) 

A.7.1  The  Rale  Equation  for  Q. 


(A-70) 


Physical  considerations  suggest  that  isotropic  damage  in  cementitious  materials  is 
caused  mainly  by  shearing  which  gives  rise  to  local  decohesion  (between  the  mortar 
and  the  aggregate,  and/or  the  paste  and  the  sand  particles).  Thus  following  similar 
reasoning  that  gave  rise  to  Eq.  (A-12)  and  noting  that  the  ultimate  value  of  Q  is 
unity,  we  set  dQ  to  be  proportional  to  a  positive  monotonic  function  of  the  isotropic 
integrity  (1-Q),  i.e.. 


dQ  =  f(l  -  G)df 


(A-71) 


where  dc  is  defined  in  terms  of  the  increment  in  the  deviatoric  strain  tensor  de..  as 

•1 


=  c2(d.ij<to.j) 


(A-72) 


and  c  >  0  is  a  material  parameter  A  suitable  form  of  the  function  f  is  as  follows: 


f (l  -  A)  =  (l  -  Q)' 


(A-73) 


where  b  >  0.  Solution  of  Eq.  (A-71)  in  the  light  of  Eq.  (A-73)  gives  rise  to  the 
expression 


1  -  Q  =  (1  -  (1  -  b)j)1_D  , 

which  is  of  the  same  form  as  Eq.  (A-26).  Specifically,  if  b  =  1  then 
(1  -  0)  =  e-b^  . 


(A-74) 


(A-75) 


We  note  that  c  in  Eq.  (A-72)  plays  the  role  of  an  isotropic  fracture  toughness 
parameter  which  influences  the  rate  of  isotropic  degradation. 

A- 20 


A. 7.2  Isotropic  Damage  in  Homogeneous  Triaxial  Strain  Fields. 


In  the  presence  of  isotropic  damage  the  counterparts  of  Eqs.  (A-20)  are: 


“2  =  X*l^kck  +  2/i  *1£1 
w 

(A-76) 

~2  ~  X^2^kek  *  2>i  *2e2 

ki 

(A-77) 

~2  =  X*3^k£k  *  *3£3 

(A-78) 

In  that  which  follows,  we  illustrate  the  application  of  the  above  equations  to 
two  special  cases  of  interest,  namely,  simple  tension  and  simple  compression. 

(A)  SIMPLE  TENSION 


Here  =  0  and  the  analysis  follows  very  closely  that  of  Section  A. 4. 2. 

Specifically,  in  view  of  Eqs.  (A-78): 


*2e2  =  '  "  *lel 


(A-79) 
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(A-80) 


However,  as  opposed  to  the  analysis  in  Section  A. 4. 2  two  functions  ^  and  u  are  now 
necessary  to  determine  the  axial  stress  In  the  absence  of  explicit  coupling 

between  directional  and  isotropic  damage,  ^  is  still  given  by  Eq.  (A-26)  in  which 
case 


1 

^  =  {l  -  (1  -  m)ke1}1_m  (A-81) 

However,  we  expect  that  k  in  Eq.  (A-81)  will  depend  on  (1  since  clearly  the  material 
resistance  to  cracking  will  be  greatly  influenced  by  the  presence  of  the  isotropic 
damage.  Specifically  we  expect  that 
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where  k  is  an  increasing  function  of  Q.  In  this  event  Eq.  (A-25)  can  no  longer  be 
integrated  directly  and  Eq.  (A-81)  will  not  hold. 

Evolution  of  Isotropic  Damage. 

In  the  presence  of  axial  loading  =  0) 


*r  = 


which  in  the  light  of  Eq.  (A-79)  becomes 


(A-83) 
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ds  *  cll 

de^  +  v  d 

k  £iJ 

(A-84) 


where  use  was  made  of  Eq.  (A-12).  Thus,  in  general,  we  have  the  following  coupled 
equations  that  govern  simultaneously  the  evolution  of  directional  and  isotropic  damage 
in  simple  tension: 


d^  =  -  k(Q)de1 


(A-85) 


dfi  =  (1  -  Q)b  cjf  (l  ♦  i*x  - 
Approximate  Solution  of  Eqs.  (A-85)  and  (A-86). 


(A-86) 


To  obtain  a  qualitative  understanding  of  the  simultaneous  effects  of  isotropic  and 
directional  damage  we  solve  Eqs.  (A-85)  and  (A-86)  approximately.  Specifically  we 
uncouple  the  equations  by  setting  k(fl)  equal  to  a  constant  k  and  setting  equal  to 
zero  the  terms  on  the  righthand  side  of  Eq.  (A-86)  that  are  multiplied  by  v,  since  in 
the  region  of  softening  0  <  1  and  the  terms  and  k/j*  are  considerably  less  than 

unity.  Then  ^  is  given  by  Eq.  (A-81)  while  fi  is  given  by  Eq.  (A-87).  i.e.. 
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(1  -  b) 


[3  I1' 

Ci 2  el) 


(A-87) 


in  view  of  Eq.  (A-74).  Thus  may  now  be  found  upon  use  of  Eqs.  (A-80).  (A-81) 
and  (A-68)  and  is  given  below: 
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(A-88) 
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=  E(1  -  (1  -  n,)k£1]1-m  ejl  -  (1  -  b)cj|  ej‘ 


In  the  very  specific  case  where  m  =  b  =  1 


a\  ■  E£ie 


*  acllk 


(A-89) 


In  other  words,  directional  and  isotropic  degradation  compound  the  effect  of  softening 
in  simple  tension.  As  we  shall  see  this  is  not  the  case  in  compression.  For 
comparison  when  isotropic  damage  is  absent,  see  Eq.  (A-27b). 


(B)  SIMPLE  COMPRESSION 


In  this  case  Eqs.  (A-79)  and  (A-80)  apply  but  now  ^  =  1.  Thus  in  this 

specific  case 


*9£9  =  *  * 


(A-90) 


(A-91) 


As  opposed  to  tension,  softening  in  compression  is  due  entirely  to  volumetric  damage. 
The  appropriate  rate  equations  in  the  presence  of  monotonic  loading  are: 


t2  =  -  *2  k(0)de2 


(A-92) 


*  “  -  °)b'fl  (‘  *  ;  -  ^  £2)*: 


(A-93) 


where  Eqs.  (A-71).  (A-73).  (A-90),  and  (A-92)  were  used.  Eqs.  (A-92)  and  (A-93) 
are  coupled  and  cannot  be  integrated  directly.  We  obtain,  therefore,  approximate 
solutions  of  these  equations  in  order  to  study  their  qualitative  trends.  To  that  end 
we  let  k(fl)  be  a  constant  k  and  define  c*  by  the  relation 


t1  *  ;  -  ;  * 


(A-94) 


1: 


Equations  (A-92)  and  (A-93)  may  now  be  integrated  to  give: 


1 

+2=  {l  -  (1  -  m)k£2}1-m  (A-95) 

1 

1  -  Q  a  {l  -  (1  -  bjc'ej1"15  (A-96) 

1 

■Wj  =  €2^2  =  -  (!  -  "0k  e2)1’ffl  (A-97) 

2a 

=  ECl{l  -  (1  -  b)c*e2}1‘b  (A-98) 


In  Figure  A. 5  we  show  an  experimental  relation  between  the  axial  and  transverse 
strains  in  the  case  of  a  concrete  block  under  axial  compression,  according  to  Van 
Mier  (1984).  and  the  analytical  relation  according  to  Eq.  (A-97)  when  v  —  0.35, 
k  =  70  and  m  =  2.4.  Also  shown  is  a  compressive  stress-strain  relation  for  the 
same  material  and  specimen.  The  approximate  analytical  axial  stress-axial  strain 
relation  obtained  on  the  baste  of  Eqs.  ^(A-97)  and  (A-98)  is  shown  in  Figure  A. 6 
when  2a/l-b  =  1.5  and  (l-b)c  =  40  (c  assumed  constant).  For  the  same  model,  a 
comparison  between  the  predicted  and  measured  volume  strain  e^k  as  a  function  of 
the  axial  strain  e1  is  shown  in  Figure  A. 7.  From  this  figure,  it  is  evident  that  the 
model  describes  the  observed  dilatancy  quite  well. 

A.7.3  Discussion. 

The  extensive  and  painstakingly  precise  investigations  by  Van  Mier  into  the 
compressive  deformation  and  fracture  of  concrete  make  it  abundantly  clear  that  the 
softening  branch  of  the  axial  stress-strain  curve  represents  the  behavior  of  a  geometric 
structure  rather  than  that  of  a  uniform  material.  The  experiments  in  tension  by 
Gopalaratnam  and  Shah  (1985)  lead  to  similar  conclusions.  There  are  therefore 
significant  statistical  variations  in  the  behavior  of  concrete  in  the  "softening  regime." 

The  continuum  damage  theory  which  we  have  proposed  in  the  previous  sections 
captures,  we  believe,  many  of  the  essential  features  of  the  constitutive  behavior  of 
fracturing  concrete.  Further  aspects  of  this  approach  shall  be  pursued  in  our  future 
research. 
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Relation  between  axial  and  transverse  strains  under  conditions  of 
simple  compression. 
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APPENDIX  B 


RECENT  ADVANCEMENTS  IN  THE  ENDOCHRONIC  CONCRETE  MODEL 

This  Appendix  describes  several  advancements  that  have  recently  been  made  in 
the  endochronic  model  described  earlier  in  Section  4.3  which  extend  its  range  of 
applicability  to  high  pressures  and  provide  a  more  realistic  description  of  the  failure 
surface.  Details  of  these  extensions  of  the  model  are  given  in  the  following  sections 
of  this  appendix. 

B.l  A  THIRD  INVARIANT  FAILURE  SURFACE  FOR  PLAIN  CONCRETE. 

Because  the  fine  details  of  the  failure  surface  were  not  expected  to  have  a 
significant  effect  on  the  behavior  of  the  endochronic  model  for  peak  stress  states 
sufficiently  far  from  failure,  a  simple  failure  surface  was  adopted  in  the  model 
discussed  in  Section  4.3.  with  the  function  Fs  taken  to  depend  linearly  on  a.  When 
this  is  the  case,  it  can  be  shown  that  the  resulting  failure  is  of  the  Drucker-Prager 
type,  with  the  trace  of  this  surface  in  the  )r-plane  being  a  circle.  Thus,  such  a 
surface  does  not  account  for  the  possible  effect  on  failure  of  the  intermediate  principal 
stress.  A  more  realistic  failure  surface  for  plain  concrete  allows  Fs  to  depend  not 
only  on  a  but  also  on  the  third  invariant  of  the  deviatoric  stress  tensor.  J'.  The 
incorporation  of  such  a  failure  surface  into  the  endochronic  model  is  described  in  this 
section  and  illustrated  through  application  to  a  suite  of  plain  concrete  failure  data. 

When  the  shear-kernel  p( zs)  is  of  such  a  form  that  a  failure  surface  exists,  the 
shear  hardening  function  F  defines  the  properties  of  the  failure  surface.  For  isotropic 
media,  the  most  general  form  of  the  failure  surface  can  be  expressed  in  terms  of  the 
three  independent  invariants  of  the  stress  tensor,  namely: 

o  =  tr(Z) 

J2  =  2  2  :  2  (B_1) 

.  iJI  J3 

2  (J2)2/3  ' 


B- 1 


where 


J3  =  3  det  (5) 


(B— 2) 


Here,  the  notation  det($)  denotes  the  determinant  of  the  components  of  s.  The  most 
general  form  of  the  failure  criterion  for  isotropic  media  can  therefore  be  expressed  in 
the  form: 

n~2  =  1  (a,J).  (B-3) 


It  was  shown  earlier  by  Hegemier.  Read.  Valanis  and  Murakami  (1986)  that  the 
hardening  function  Fs  can  be  expressed  in  terms  of  the  function  fJ(a.J)  as  follows: 


=  H  i  fa.j) 

s  M 

w 


where 


(B-5) 


Therefore,  to  define  Fs  for  a  particular  material,  the  function  f(oJ)  must  be  specified 
This  is  done  below  for  plain  concrete 

There  are  a  number  of  advanced  failure  criteria  for  plain  concrete  which  have 
been  proposed  in  the  literature  that  have  the  general  form  of  Eq.  (B-3)  (Mills  and 
Zimmerman.  1970:  Wiliam  and  Warnke,  1974:  Ottosen.  1977:  Lade.  1982: 
Peyton.  1983;  Podgorski.  1985).  Considering  the  large  scatter  in  data  between 
experimenters,  as  well  as  between  different  testing  devices  (see  Hegemier  and 
Read.  1986).  most  of  the  advanced  failure  criteria  noted  above  provide  satisfactory 
descriptions  of  existing  plain  concrete  failure  data.  For  critical  assessments  of  most 
of  the  above  failure  models,  the  reader  is  referred  to  Ottosen  (1977)  and 
Podgorski  (1985). 


Inasmuch  as  all  of  the  failure  models  noted  above  are  reasonably  accurate,  and 
none  appear  to  have  significant  advantages  in  predictive  capability  over  the  others,  we 
select  the  model  proposed  by  Peyton  (1983)  for  specific  consideration  in  this  section 
because  of  our  greater  familiarity  with  it.  A  description  of  this  failure  model  is  given 
below,  where  the  corresponding  form  of  the  function  if(ff.J)  is  defined.  The  model  is 
then  applied  to  failure  data  for  plain  concrete. 


The  failure  criterion  proposed  by  Peyton  (1983)  is  given  by  the  expression: 


(B-6) 


where  r  and  a  are.  in  general,  smooth  monotonic  functions  of  the  hydrostatic  stress 
a.  i.e.. 

r  =  r(a)  . 

(B-7) 

a  =  a(a)  . 

The  role  of  the  function  t(o)  is  to  describe  the  effect  of  pressure  on  the  meridians  of 
the  failure  surface,  while  the  function  a(<r)  determines  the  manner  in  which  the  trace 
of  the  failure  surface  in  the  x-plane  changes  with  a.  In  order  for  the  trace  to  be 
convex,  it  is  necessary  that  0<a$l.  If  a  is  a  continuous  decreasing  function  of  a.  the 
trace  of  the  failure  surface  in  the  jr-plane  changes  smoothly  from  a  triangle  (a=l)  to 
a  circle  (a=o)  with  increasing  a. 


To  express  Eq.  (B-6)  in  the  form  of  Eq.  (B-3).  we  note  that  by  defining  an 
angle  6  in  the  x-plane  as  shown  in  Figure  B-l  it  follows  that: 

J  =  -sin  3 8  (B-8) 

Solving  Eq.  (B-lc)  for  J3  and  combining  the  result  with  Eq.  (B-6)  gives 


3/2 


3f3 


=  1  , 


(B-9) 


which  is  of  the  form: 


if  we  set: 


With  the  change  of  variable 
1 

x  =  - 

y 

Eq.  (B-10)  can  be  written  in  the  form: 

y3  -  y  +  1  =  0  . 

With  a  further  change  of  variable. 

y  =  X  sin  l  . 

Eq.  (B-13)  becomes 

sin3  ^-4$in^+%so  . 

\Z  X 6 


(B-10) 


(B-ll) 


(B-12) 


(B-13) 


(B-14) 


(B-15) 


We  observe  that  this  equation  has  the  form  of  the  following  trigonometric  identity: 


3  3  i 

sin  s  i  n  ^  j  s  i  n  3  ^  =  o  , 


(B- 16) 


if  we  set 

_2 

X  =  ,  sin  3  ^  =  a  J 


(B-17) 


B-5 


From  Eq.  (B-17b).  it  follows  that 


*  =  \ 


sin”*  (aJ)  ♦  2  f 


(B— 18) 


Upon  combining  Eqs.  (B-llb).  (B-12).  (B-14).  (B-17)  and  (B-18).  we  can  write 

^  * ,  .  ir .  -i,  ~  2,i  ■ 

2  sin  1^  sin  (aJ)  * 


which  is  of  the  form  of  Eq.  (B-3).  To  apply  this  expression  to  a  specific  set  of 
data,  the  functions  a(a)  and  t[o)  must  be  specified.  A  procedure  for  doing  this  is 
described  below. 

Generally,  failure  data  are  obtained  from  triaxial  compression  and  triaxial 
extension  tests  and.  with  this  in  mind,  we  define  the  ratio  r  as 


r  = 


(f3c 
(F3  ' 


(B-20) 


where  the  subscripts  c  and  e  refer,  respectively,  to  compression  and  extension  and 
both  values  of  are  taken  at  the  same  hydrostatic  pressure.  Considering  Eq.  (B- 
19)  with  J  =  +1  for  compression  and  J  =  -1  for  extension,  it  can  be  shown  that 


r  = 


f  •  2r 

1 

.  -1  i 

[sin  -3  - 

3 

s  1  n 

“J 

f  •  2ir 

1 

•  -1  i 

(sin  -3  ♦ 

3 

sin 

“J 

(B— 21) 


Upon  solving  for  a,  we  find 
a  =  sin  |3  tan  * 


(B-22) 


B-6 


Therefore,  if  we  know  the  manner  in  which  (Jj^  and  (4~jT)e  depend  on  a  from  data, 
then  r(o)  is  known  from  Eq.  (B-20)  and  a(a)  can  be  found  from  Eq.  (B-22). 

To  determing  t(j).  we  solve  Eq.  (B-19)  for  r(a)  and  set  J  =  +1  to  obtain 


r(a) 


(B-23) 


Therefore,  since  a(a)  and  the  variation  of  (4-J^)c  with  a  are  known,  the  function  r(a) 
can  be  found. 

The  procedure  described  above  was  used  to  determine  the  forms  of  a(a)  and 
t(o)  for  plain  concrete,  using  the  failure  data  shown  in  Figure  B-2.  Here,  data  from 
both  triaxial  compression  and  triaxial  extension  tests  on  plain  concretes  conducted  by 
several  investigators  on  several  different  strengths  of  concrete  are  depicted.  In  this 
manner,  it  was  found  that  the  data  could  be  well  represented  if  we  set 

a  =  a 

(B— 24) 

T  ~  To  *  Pj~o  -  70  » 


where 


a  =  0.83 

o 

r  =  -1  91  f 

o  c 


(B-25) 


p  =  1.24  If 

o  c 

7  =  -2  69  f 

'o  c 


The  excellent  ability  of  the  failure  model,  with  only  four  material  parameters,  to 
correlate  the  failure  data  is  shown  in  Figure  B-2 


♦  ■Q# 


Correlation  of  concrete  fracture  data  by  present  failure  model. 


•49  VDS 


To  determine  the  corresponding  expression  for  the  shear  hardening  function  F  . 
we  first  note  that,  since  the  shear  kerne!  function  p( z  )  is  assumed  to  be  of  the  form 
of  Eq.  (6.7).  it  follows  from  Eq.  (B-5)  that 


M„= 


(B-26) 


Upon  substituting  the  values  for  and  a r  given  by  Valanis  and  Read  (1986).  we  find 


M_  =  17.5  ksi 


(B-27) 


Equations  (B-3).  (B-4)  and  (B-19)  may  be  combined  to  give  the  following  expression 
for  F  : 


F  _  - — 9 - 

s  .1.-1,..  2r 

sin  z  sin  (aJ)  —z 


(B-28) 


where  a(a)  and  r(a)  are  defined  by  Eqs.  (B-24)  and  (B-25)  and  cq  =  0.07  ksi’1. 

B.2  A  HYDROSTATIC  MODEL  FOR  HIGH  PRESSURES. 

Under  monotonically  increasing  pressure,  it  is  observed,  as  illustrated  in 
Figure  B-3.  that  the  hydrostatic  behavior  of  plain  concrete  is  characterized  by  an 
initially  convex  stress-strain  curve,  which  eventually  becomes  concave  and  finally 
(asymptotically)  linear.  The  linear  portion  is  reversible,  indicating  that  the  material 
has  reached  an  elastic  state.  As  the  elastic  state  is  approached,  the  stress  and  slope 
of  the  hydrostatic  pressure  -  plastic  volumetric  strain  curve  tend  to  infinity  as  the 
plastic  volumetric  strain  tends  toward  a  critical  value.  ep.  which  corresponds  to 
complete  compaction.  Further  compression  occurs  elastically  with  dep  =  o. 


The  role  of  the  hardening  function  FH  is  to  account  for  the  effect  of  compaction 
on  hydrostatic  behavior  and,  since  compaction  is  reflected  through  the  plastic 
volumetric  strain,  ep,  it  is  natural  to  take  F  =  FL1(ep).  It  can  be  shown  that,  for 

rl  n v  ' 


monotonic  loading  above  the  initial  convex  portion  of  the  stress-strain  curve,  the 
hydrostatic  stress  a  and  FH(ep)  are  related  by  the  expression 


a  = 


=  Ko  Fh(£P)  • 


(6-29) 


where  K  is  a  positive  material  constant.  Therefore,  except  for  the  constant  multiplier 


Kq.  the  function  FH(ep)  describes  the  monotonic  hydrostatic  loading  curve  for  a 


II 


pressures  above  the  initial  convex  part  of  this  curve. 


The  general  restrictions  that  the  function  F,.  must  satisfy,  in  addition  to  being  a 

D  ** 

monotonically  increasing  function  of  e  .  are: 


fh<°>  =  1 


(B-30) 


lim  F  =  « 

e'V,H 

C 


The  first  restriction  is  imposed  simply  for  convenience  and  allows  us  to  set  Kq  =  a q. 
where  aQ  is  the  intercept  on  the  a- axis  shown  in  Figure  B.13.  With  Kq  set  to  in 

Eq.  (B-29).  the  virgin  hydrostatic  compression  curve  from  A  to  C  is  then  described  by 


the  equation 

"  =  fh(£P)  ■ 

which  proves  convenient  in  applying  the  model  to  specific  materials. 


(B-31) 


The  second  restriction  arises  from  the  fact  that  it  takes  an  infinite  pressure  to 
close  all  pores,  that  is: 

(B-32) 

In  view  of  Eq.  (B-31).  condition  (B-32)  is  satisfied  by  taking  F^  such  that 
condition  (B-30b)  is  fulfilled. 

The  hydrostatic  data  to  which  the  endochronic  model  discussed  in  Section  4.1 
was  applied  covered  the  range  of  pressures  up  to  15  ksi.  Over  this  range,  it  was 
found  that  a  simple  linear  expression  for  FH,  which  satisfies  condition  (B-30a),  led  to 
an  excellent  description  of  the  hydrostatic  data.  For  higher  pressures,  however,  this 


lim  i 
ep+e 


linear  form  for  FH  is  inadequate  and  a  more  general  expression  which  also  satisfies 
condition  (B-30b)  is  required. 

To  extend  the  range  of  the  model  to  high  oressures.  including  elastic  behavior 
beyond  full  compaction,  a  suitable  expression  for  Fu  is  as  follows: 

ri 

FH  =  1  +  pS  +  7o  (<£  -  eT™  (B-33) 

where  /Jq.  and  m  are  positive  constants.  This  expression  satisfies  conditions  (B- 
30)  and  is  a  monotonically  increasing  function  of  ep.  There  are.  of  course,  other 
expressions  that  satisfy  these  general  requirements,  also,  but  that  given  by  Eq.  (B-33) 
appears  to  capture  the  behavior  of  plain  concrete  quite  well,  as  shown  in  Figure  B-4. 
Here,  the  response  of  the  endochronic  hydrostatic  model,  with  FH  defined  according  to 
Eq.  (B-33).  to  pure  hydrostatic  compression  of  plain  concrete  up  to  65  ksi  is 
compared  with  the  corresponding  data  reported  by  Green  and  Swanson  (1973).  which 
unfortunately  provided  very  limited  unloading  data.  The  agreement  between  the  model 
and  the  data  is  considered  quite  good. 

Finally,  if  pressures  beyond  point  B  in  Figure  B.3  are  not  of  interest,  the 
following  simple  expression  for  Fu  has  been  found  to  be  adequate  in  many  cases: 

Fh  =  </eP  (B-34) 

This  expression  satisfies  only  the  first  of  the  two  conditions  in  Eq.  (B-30)  and.  it 
reduces,  for  small  ep.  to  the  linear  expression  for  Fu  adopted  in  the  model  discussed 

n 

in  Section  4.3. 
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APPLICATION  OF  DIRECT  VARIATIONAL  METHOD  TO 
DENSE  STEEL  LAYOUT  -  3D  PROBLEM 

C.l  INTRODUCTION. 

This  Appendix  furnishes  a  detailed  derivation  of  an  advanced  mixture  model  for 
reinforced  concrete  ( R / C) .  The  development  is  based  on  the  use  of  the  Direct 
Variational  Method. 

The  theoretical  construction  process  to  be  discussed  involves  a  number  of  key 
steps;  these  include:  (1)  a  judicious  scaling  of  the  problem  variables;  (2)  introduction 
of  so-called  microcoordinates  and  a  two-space,  asymptotic  representation  of  all  field 
variables:  (3)  a  local  periodicity  condition;  and  (4)  introduction  of  a  variational 
principle.  Each  step  is  covered  in  some  detail. 

C.2  BASIC  RELATIONS. 

Consider  a  uniaxial  steel  layout  with  an  initially  local  periodic  array  occupying  a 
volume  Vq  with  boundary  9Vq  in  the  initial  configuration  as  illustrated  in  Figure  C.l. 
In  the  current  configuration  the  corresponding  volume  and  surface  will  be  denoted  V 
and  9V.  respectively. 

If  a  spatial  (Eulerian)  description  of  motion  is  adopted,  and  if  x  denotes  the 
global  position  vector  of  material  particles,  then  the  basic  conservation  equations  of 
each  constituent  are: 


D/>(a)/Dt  ♦  =  0  in  V  , 

(C-l) 

V  •  2^a)  ♦  =  p^D^Dt  in  V  , 

(C-2a) 

k-'-.i 


Figure  C.l  Densely  reinforced  concrete;  initial  configuration. 
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where  ^  is  the  gradient  operator  with  respect  to  x .  ^  is  the  velocity  vector. 
dH  =  (l/2)(v  +  v.  .)  is  the  rate-of-deformation  tensor,  g  denotes  the  Cauchy  stress 

tensor,  p  denotes  material  density,  f  is  a  body  force  per  unit  (current)  volume.  T  is 
the  traction  vector  acting  on  the  surface  9V  with  unit  outer  normal  f/.  e  denotes 
internal  energy  per  unit  mass.  F  is  the  heat  flux  vector  per  unit  time  and  current 
area,  q  is  a  heat  source  rate  per  unit  time  and  mass,  superscripts  a  =  1,2  denote 
material  a  with  a  =  1  representing  steel  and  a  =  2  concrete.  (  )T  denotes 
transposition  of  the  tensor  indices.  D(  )/Dt  is  the  material  time  derivative,  and 
is  the  current  volume  occupied  by  material  a  with  boundary  3V^.  In  what  follows, 
both  direct  and  index  notations  will  be  used.  Lower  case  indices  refer  to  a  fixed 
rectangular  Cartesian  coordinate  system  and  repeated  subscripts  are  summed. 
Superscripts  do  not  refer  to  tensor  components.  Repetition  of  unbracketed 
superscripts  implies  summation.  Bracketed  superscripts  are  not  summed.  Thus,  a  is 
not  summed  in  the  foregoing  equations.  Equations  (C-l)  -  ( C-4 )  represent  the 
conservation  of  mass,  linear  momentum,  angular  momentum,  and  energy,  respectively. 


In  (C-2).  the  components  of  V«g  and  j/«g  are  a...  and  v.a...  respectively.  In  (C-4)  the 
components  of  g:$i  are 


In  what  follows,  attention  will  be  focused  upon  physical  processes  that  occur 
over  time  intervals  that  are  small  compared  to  typical  thermal  diffusion  times. 
Consequently,  an  adiabatic  premise  will  be  adopted.  Accordingly,  in  (C-4)  one  sets 


q(fl)  ~  o  ,  £^2  0 


(C— 5) 


The  conservation  equations  must  be  supplemented  by  appropriate  constitutive 
relations.  A  reasonably  general  class  of  constitutive  behaviors  may  be  characterized 
by  assuming  that  an  objective  stress  rate  relates  linearly  to  the  rate-of-deformation 
tensor,  i.e.. 


{W.5W:iW 


(C-6) 


where  is  the  instantaneous  elastic-plastic  modulus  tensor  for  material  a  and  the 


components  of  are  D..k^dk^.  The  tensor  is  assumed  to  be  independent  of 
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the  stress  rate  and  the  rate  of  deformation,  but  may  depend  upon  the  current  stress 
state,  internal  energy,  and  certain  deformation  measures.  Isotropy  of  may  not  be 
assumed  in  general  due  to  plastic  flow  and/or  prior  cracking  of  the  concrete.  While 
is  symmetric  with  respect  to  the  exchange  of  i  and  j,  and  k  and  t,  it  is  not. 
in  general,  symmetric  with  respect  to  the  exchange  of  ij  and  k£_ 

The  stress  rate  in  (C-6)  deserves  comment.  The  Jaumann  rate  corotational  with 
the  material  is  most  frequently  adopted  in  constitutive  relations  of  the  type  (C-6).  In 
this  case 

=  Dg^/Dt  -  g^  ♦  g^»  w(a>  (C-7) 


where  components  of.  say.  w*g  are  w.^..  However,  a  number  of  other  objective 
stress  rates  have  been  introduced  by  various  authors.  While  each  rate  is  admissible 
from  the  standpoint  of  being  objective,  each  leads  to  a  different  constitutive  relation. 
An  appropriate  choice  depends  upon  the  problem  under  consideration. 

The  constitutive  form  (C-6)  covers  a  wide  class  of  elastic-plastic  and  elastic- 
plastic-brittle  fracture  behavior.  Yet.  it  is  not  intended  to  be  comprehensive.  Other 
constitutive  relations  can  be  incorporated  into  the  theoretical  development  to  follow. 
The  incremental  form  (C-6)  is  the  most  convenient  for  demonstration  purposes.  Some 
important  forms  such  as  the  endochronic  theory,  which  is  incrementally  nonlinear, 
necessitate  modification  of  (C-6). 

A  remark  concerning  the  foregoing  Eulerian  description  of  motion  is  appropriate 
at  this  point.  For  computational  purposes,  it  is  often  convenient  to  employ  an 
"updated  Lagrangian"  description  of  motion.  Accordingly,  the  Kirchhoff  stress  is 
adopted  as  the  stress  measure  and  the  current  configuration  is  used  as  the  reference. 
In  view  of  the  latter. 


r 


where  r  denotes  the  Kirchhoff  stress  tensor.  However.  Dr/Dt  i  Dg/Dt  since 


(C-8) 


t 


Dr/Dt  =  Dg/Dt  +  dg  ,  d  =  V  •  v 


i 


(C-9) 


Thus,  for  an  updated  Lagrangian  description,  one  need  only  replace  g  by  r  in  the 
conservation  equations  (C-l)  -  ( C-4 )  and  interpret  D(  )/Dt  as  a  material  time 
derivative.  However,  the  constitutive  relations  (C-6)  must  be  altered.  For  this 
purpose  (C-6)  is  replaced  by 


.CW:iW  ■ 


1 

(C-10)  ; 


From  (C-6).  (C-8).  and  (C-9).  one  finds 

B00  *  £«  .  jWj  ,  or  t(a)  =  B<a)  *  (C-ll) 

where  gj5  denotes  the  tensor  product,  i.e..  components  of  gg  are  0-6^  (this  product 
is  often  written  as  g  ®  £). 

In  addition  to  the  above  conservation  and  constitutive  equations,  relations 
between  tractions  and  velocities  for  each  constituent  across  the  steel-concrete  interface 
must  be  specified.  For  modeling  purposes  this  interface  is  idealized  as  smooth.  Let 
s/*  denote  that  part  of  9V^  in  contact  with  9V'  \  where  both  are  smooth.  Then, 
continuity  of  the  traction  vector  requires  that 


—  * 

on  . 


(C-12) 


With  respect  to  displacements,  relative  slip  across  ./*  must  be  allowed  to  simulate 
steel-concrete  bond  slip  and  dowel  action.  Accordingly,  a  relative  velocity  [v]  is 
defined  by 


on  ./ 


(C-13) 


An  interface  "constitutive  law"  is  jiext  postulated  in  the  form  of  a  linear  relation 
between  an  objective  traction  rate,  T*.  and  the  relative  velocity,  i.e.. 


•ft 


T  5  DT  /Dt  -  w  •  T  , 


(C-14b) 


where  ■■£  is  a  second  rank  tangent  modulus  tensor  which  may  depend  on  the  current 
states  of  stress  and  relative  slip  on  ■'/* . 

It  should  be  noted  that  the  normal  component  of  velocity  (to  the  surface  ■'/*) 
may  be  discontinuous  as  well  as  the  tangential  component  in  (C-14).  This  degree  of 
freedom  is  necessary  to  allow  proper  simulation  of  actual  steel-concrete  interaction 
under  certain  conditions  using  the  idealized  smooth  surface  •'/* . 

If  contact  between  9V^  and  9V^  is  lost,  then  one  sets  =  0.  =  0 

over  those  segments  of  9V^  and  9V^  corresponding  to  no  contact.  This  situation 
can  also  be  incorporated  into  the  form  (C-14a)  by  a  suitable  choice  of  ■£. 

The  foregoing  relations.  ( C- 1)  -  (C-14).  must  be  supplemented  by  initial 
conditions  at  t  =  0  and  boundary  conditions  on  9V. 

C.3  SCALING. 

It  will  be  convenient  to  nondimensionalize  the  basic  equations.  For  this  purpose 


=  typical  macrosignal  wavelength 


=  typical  steel  spacing 


c,  ,.p,  .  =  reference  wave  velocity  and  macrodensity 

(m)  r(m) 

E.  .  =  p.  .c?  .  =  reference  modulus 

M  M  (m) 


typical  macrosignal  travel  time  (=  A/c. 


=  A/A  =  ratio  of  micro-to-macrodimensions 


With  the  aid  of  the  above  notation,  nondimensional  variables  are  now  introduced 
according  to 


=  */*•  '  =  '/'h  • 

=  ®/W  x*  =  r/E(m). 

-  8/E(m).  (l/^=-S/(E|ml4). 

=  p(m)  ■  l  = 


(C-15) 


§s 
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C.4  MICROCOORDINATES. 

It  is  expected  that  stress  and  deformation  fields  will  vary  significantly  with 
respect  to  two  basic  length  scales:  (1)  a  "global"  or  "macro"  length  typical  of  the 
body  size  or  loading  condition:  and  (2)  a  "micro"  length  typical  of  a  "cell"  planar 
dimensions  as  depicted  in  Figure  C.2.  These  macro  and  micro  dimensions  will  be 
associated  with  the  variables  A.  A.  respectively.  Further,  it  is  expected  that  these 
scales  will  differ  by  at  least  one  order  of  magnitude  in  most  cases.  This  suggests 
the  use  of  multivariable  asymptotic  techniques,  Hegemier,  et  at.  (1979).  This 
approach  commences  by  introducing  new  independent  variables  according  to 


x*  =  f  (02  .  0(0  0  as  e  ■»  0 


2  =  0(Ox,  0(e)  -  1  as  e  ♦  0  . 


(C-16) 


For  the  present  analysis  it  will  suffice  to  set 
0(e)  =  e  ,  0(e)  =  1  ; 


(f-in 


thus,  all  field  variables  f(x.t)  are  now  functions  of  the  "microcoordinates  j*  iv  a. 
as  the  "macrocoordinates"  x  =  x.  i.e.. 


A  V  A  A 


(C-18a) 


Spatial  derivatives  of  a  function  f(j$.t)  then  take  the  form 

v  *  v  ♦  «'Yf  • 


(C-18b) 


For  notational  convenience  in  the  following,  the  functions  f  and  F  will  both  be  written 
as  f. 


For  some  problems  it  is  necessary  to  introduce  time-microvariables  as  well  as 
the  foregoing  space-microvariables.  These  additional  variables  will  not  be 
advantageous  for  the  present  formulation,  however. 


It  is  noted  that,  for  the  present  unidirectional  steel  layout,  it  is  appropriate  to 
constrain  the  vector  2*  according  to 


a*U  =  0 

where  t  denotes  the  tangent  vector  to  the  steel  centerline. 

C.5  SYNTHESIZED  FIELD  EQUATIONS. 

(C-18c) 

The  operations  (C-18),  when  applied  to  all  field  variables,  lead  to 
synthesized  governing  field  equations  in  nondimensionai  form: 

(a)  Conservation  Relations 

the  following 

D/j^/Dt  ♦  P^l.  •  ^  =  0  in  V^, 

(C— 19) 

l  •  2(a)  ♦  7  i-  e(a)  *  £(a)  =  0  in  VW  , 

(C-20a) 

Iw  -  •  e(a)  on  wW  , 

(C-20b) 

«W  •  2(a,T  VW  . 

(C-21) 

,«D.W/Dt  .  \  !*.<«)  =  ,«:(*«  .  i'W)  in  V<“> 

C-9 

(C-22a) 

sf(<0  =  \l  (j*2  *  !*Tj)  •  (C— 22b) 

(b)  Constitutive  Relations 

8(“)  ,  jjW;J(«)(  jj(«>  5  ,)(*>.  jj’to  inV(“)  (C— 23) 


(c)  Interface  Relations 


on  & 


(C-24a) 


w 


on  &  . 


(C-24b) 


(d)  Initial  Conditions  at  t  =  0.  Boundary  Data  on  8V. 

In  the  above.  £*  is  the  gradient  operator  with  respect  to  5*,  e.g..  components  of 
are  dv./8x*.  The  material  derivative  D()/Dt  refers  to  the  macrovariables  x.  only. 

C.6  CONSTRAINTS. 

In  the  following,  attention  will  be  restricted  to  processes  in  which: 

(1)  The  constitutive  and  interface  relations  do  not  depend  on  the 
internal  energy. 

(2)  The  microdeformation  gradient  is  sufficiently  small  that  micro-flux 
or  micro-convection  terms  may  be  neglected  in  the  material 
derivative. 

(3)  The  micro-rotation  is  sufficiently  small  that  its  contribution  to  the 
stress  rate  may  be  neglected. 

(4)  The  micro-deformation  is  sufficiently  small  that  microvariations  in 
the  mass  density  may  be  neglected. 


C-10 


Under  the  above  constraints,  which  are  considered  to  be  appropriate  for  R/C.  the 
last  term  in  (C-19)  can  be  neglected,  the  energy  equation  (C-22a)  need  not  be 
treated,  and  the  material  and  stress-rate  derivatives  are  to  be  interpreted  in  terms  of 
the  macrocoordinates  only. 

C.7  PERIODICITY  CONDITION. 

An  important  premise,  called  the  local  periodicity  condition,  is  now  introduced. 
This  condition  consists  of  the  assumption  that  local  periodicity  in  the  microvariable  x* 
may  be  invoked  for  ail  field  variables.  This  premise  allows  one  to  analyze  a  single 
cell  in  an  effort  to  determine  the  distribution  of  any  field  variable  with  respect  to  the 
microcoordinate  g*.  A  typical  such  cell  is  depicted  in  Figure  C.2. 

The  local  periodicity  premise  is  motivated  by  extrapolation  from  linear  analysis. 
It  is  expected  to  provide  a  good  representation  of  the  field  variables  for  dense  steel 
layouts. 

C.8  DIRECT  VARIATIONAL  PRINCIPLE  FOR  SYNTHESIZED  FIELDS. 

A  variational  (weighted  residual)  principle  is  now  constructed  for  the  synthesized 
field  variables.  For  this  purpose  consider  the  relation 


In  the  above.  d£*  =  dA*£,  where  £  is  the  unit  tangent  vector  to  the  steel  centerline 


in  a  given  cell.  Thus,  if  a  plane  with  normal  £  intersects  the  steel  centerline  within 
a  cell,  then  is  the  projected  area  of  material  a  on  this  plane.  The  integration  is 
carried  out  over  the  microcoordinates,  e.g..  if  x*  k  =  2.3  denote  local  Cartesian 
coordinates  in  the  above  plane,  then  dA*  =  dx*dx*.  The  quantity^  =  9A^.  i.e..'^ 
is  the  projection  of  the  steel-concrete  boundary  on  this  plane. 

Let  the  boundary  9V  =  9Vy  4-  9Vu-  Then,  the  following  is  assumed: 


T^  =  prescribed  on  9Vy 


=  0  on  9V 
~  u 


(C-26) 


In  addition,  let 


?(1)  -  _  t(2) 


X*  on  9A(1), 


*  fl} 

T  s  prescribed  on  9AV  J 


(C-27) 


Finally,  let 


eC*  on  A^  with  respect  to  ^<2  > 


(a.)  * 

<5vv  '  eg  -  periodic. 


(C-28) 


If  '  is  the  variation  of  the  exact  velocity  field,  then  the  first  term  of  (C-25) 
vanishes  due  to  (C-20a),  the  second  vanishes  due  to  (C-20b)  and  (C-26).  the  third 
vanishes  due  to  (C-20b),  and  the  last  term  vanishes  due  to  (C-24).  Thus,  the 
validity  of  (C-25)  is  established  if  is  the  exact  velocity  field. 


If  jj*  J  is  not  the  exact  field,  then  (C-25)  represents  a  weighted  residual 
procedure  for  generating  weak  solutions  of  the  governing  equations.  The  weighting 
functions  are  In  the  latter  case  the  solutions  are  of  course  approximate. 


With  the  use  of  Gauss'  Theorem,  and  (C-26)  -  (C-28).  (C-25)  can  be  placed  in 
the  following  form: 


♦ 


fe**S[v]ds 


lW-  S.WdA' 

dV  ♦ 

• 

Z 

U=i 
Jv  J 

*W 

• 

av. 

L«*i 

• 

r 

*w 

(C— 29) 


Equation  (C-29)  constitutes  the  principle  of  Virtual  Work  for  the  synthesized  fields. 
This  relation  represents  an  important  tool,  in  the  form  of  a  variational  principle,  for 
generating  mixture  field  equations  for  the  highly  nonlinear  class  of  problems  of  interest 
herein. 

C.9  MODEL  CONSTRUCTION:  DIRECT  METHOD. 

Kantorovich's  direct  variational  method  is  now  applied  to  generate  the  desired 
mixture  equations.  The  trial  velocity  field  (test  functions)  is  selected  in  the  form 

(***»«*;«]  =  S(fl)  (*>*)  ♦  9K^“V)  ♦  °(e2)  (C-3°) 

The  "weight”  functions  gK^(j*)  are  to  be  prescribed  for  each  material  (a  =  1.2). 
The  vector  0^  reflects  the  macro  contribution  to  the  velocity  field  whereas  the 
vectors  reflect  the  micro  contribution:  these  variables  play  the  role  of  new 

|  dependent-type  displacement  variables.  In  the  above,  and  in  what  follows,  the 

|  summation  convention  is  used  with  respect  to  the  index  K  with  range  1  to  n. 

1 

I  C.9.1  Equations  of  Motion.  Boundary  Conditions. 

! 

Substitution  of  the  test  function  (C-30)  into  the  variational  principle  (C-25) 
furnishes,  after  some  algebraic  manipulations  and  integrations  by  parts,  the  following 
relation: 

i 


C-13 


IS*"'-  e*-  • « •  &“*-  ■  •  >”w)l“ 


5(a)*  (2  •  M(a)*  E(a)  ♦ 


(C— 31) 


.  (1  •  /(«)-  gXW  .  £“<•>.  (-l)fl‘1gK)]dV  =  0  . 

4 

In  (C-30).  the  following  quantities  have  been  defined: 

(ew.?w]-[  (m9k(‘>)|(*V  . 

JA(«) 

[|M/W]  .  [  (x.€9K(a))2(4)«-  , 

J*« 


(C— 31) 


J  k  (<*) 


E  5  7  %*<***  >  8K  =  <S*  gK(a)ds*  (a  =  1  or  2) 


The  Euler-Lagrange  equations  associated  with  the  variational  equation  (C-30)  are 
X  •  ♦  £(a)  ♦  (-l)a+1£  =  0  in  V  ( 


X  •  ♦  EK(a^  ♦  (-1)0*1  8K  =  0  in  V 


(C— 32) 


(C-33) 


C-14 


The  boundary  conditions  are 


f(B)  =  *  •  »<“> 

or  =0  on  9V  , 

(C-34) 

f  <*>  ,  K  .  jjK(a) 

or  =0  on  3V  . 

(C-35) 

Equations  (C-32)  for  a  =  1.2  are  the  "smoothed"  equilibrium  relations.  They 
have  the  form  of  a  binary  mixture.  The  dependent  variables  represent  "partial 
stresses",  i.e..  stress  averages  over  the  cell  area  occupied  by  the  appropriate 
constituent.  The  quantity  £  represents  a  stress  interaction  vector  which  reflects 
stress  transfer  between  the  steel  and  the  concrete  across  the  interface. 

In  addition  to  direct  averages  of  the  stress  field  over  the  cell,  weighted  averages 
also  enter  the  formulation  via  the  variational  principle.  Equations  (C-33)  for  a  =  1.2 
govern  these  new  dependent  variables. 

5 

The  variational  principle  provides  appropriate  boundary  conditions  for  the 
displacement-type  variables  and  the  stress-type  variables 

These  conditions  are  expressed  by  (C-34)  or  (C-35). 


C.9.2  Mixture  Constitutive  Relations. 

The  test  functions  (C-30)  furnish  the  following  rate-of-deformation  tensor: 

d(“)  =!<°) .  egK  <“&(«> .  jj'w 


(C— 32) 


where 

i(“>  *  I  •  zTS(a))  .4K(0)  *  *  «TiK(a))  • 

i  -*  \  (sK(aYgK(0)  *  iK(a)j‘TgK(0))  ■ 


(C-33) 


C-15 


Substitution  of  (C-32)  into  the  constitutive  relation  (C-23).  and  the  latter  into 
(C-31b.c)  gives 


!<«>.  E(“>:  |C«)  .  jK(«)  ,J(a)  t  bK(.)  .  |K(«)  _ 

gK(«)  =  4K(“):|C«)  .  SKL(.):jjL(«)  .  £KL(„)  .  £(.(*)_ 

jj(*>  =,  /<«>  ;|(“)  .  feKL(a)^L(a)  ,  CKL(,»  .  |K(.) 


(C-34) 

(C-35) 

(C-36) 


where  i.  a.  §.  §.  g.  a.  fe.  £  denote  weighted  averages  of  the  material  tensor  Q  over 
the  cell.  The  latter  quantities  are  given  by 


6(«) _*!<(«)_  gK(«)i  sKL(.)j£KL(«)]  ,  |  ^  jB(«)  eD(«)gK(a) 


.VWg'We'”,  £9Kwew.2,gL(“>]d*- ; 

[/'">. |>kl<‘,>,£kl"»]  ,  f  [£xV(<bw,£V(“Y3k(“V“>, 

J*(°) 


(C-37a) 


£J*  gK(«).B(«).j*gt.(a)]dA. 


(C-37b) 


C.9.3  Interaction  Terms. 

Using  (C-31e),  (C-30).  and  (C-24b),  the  interaction  term  £  can  be  expressed  in 
the  form 


C-16 


5*2”  •  [0]  .  •  ■  (C-38a) 

f~-  >?  -[u]  ,  E  UK(1)«*’  •  (C-38b) 

C.9.4  Summary  of  Mixture  Equations. 

The  basic  equations  of  the  mixture  model  for  dense,  unidirectionally  reinforced 
concrete  can  now  be  summarized  as  follows: 

(a)  Equations  of  Motion 

S*H(<0  *  E(a)  *  (-i)*4lE  «  o  ,  (c-39) 

V-UK(a)-  gKC“)  .  EK(o)  .  (-l)a,19K  =  0  ;  (C-40) 

(b)  Constitutive  Equations 

g<»>  ,  .  bK(.).|K(«)  (  (C-41.) 

|K(«)  .  jKW  ;|.(«)  .  5KL(4):JL(o)  .  £KL(4).|L(a)  (C-41b) 

|K(«)  =  /(«)  .J <•)  .  feKL(«)  ^(«)  4  SKL(«).|L(<1)  (c.41c) 

(c)  Interaction  Terms 

l  ■  [S] .  SK  (c-42) 


c-ir 


(d)  Boundary  Conditions  on  dV 


I 


tl« 


$ 


j(«)  ,  r(j(“)  or  jy(°)  =0  on  9Vt  , 

iKW  =  rBKW  or  =  0  on  9VT 

(e)  Kinematic  Relations 

|W  ,  JfoM  .  ,T#M)  ,4K(.)  =  1(j5K(o)  .  ,Tf(o,] 

C.9.5  Weighting  Functions. 


(C-42b) 


(C-42b) 


(C-43a,b) 


In  order  to  utilize  the  proceeding  mixture  equations,  the  weighting  functions 
gK^(i*)  must  be  specified.  For  a  dense  steel  layout  with  a  cell  aspect  ratio  near 
unity,  an  adequate  description  of  the  microdeformation  can  be  obtained  in  conjunction 
with  a  "concentric  cylinders  approximation”.  Accordingly,  the  outer  cell  boundary  is 
approximated  by  a  circle  which  is  defined  based  upon  equal  area  of  the  original  and 
approximate  cell.  This  procedure  is  depicted  in  Figure  C.3  (for  a  hexagonal  steel 
array). 

Subsequent  to  the  introduction  of  the  approximate  cell,  the  weighting  functions 
are  defined  as  follows: 

41)(*2-xs)  *  *J/"(1)  • 


9(2)(v*;)  =  -  *2  - 


*2  *2 
x  2  *  x  3J 


g31)(x2'x3]  =  x3/n(1) 


(C-44) 


* 

f  *  *) 

*  *3 

lx2'x3J  =  ‘ 

x3  ‘  *2  *2 

x2  +  x3  . 

Vt  V  ;  V.  V  t  . 


Figure  C.3  Circular  cylinders  approximation  (shown  for 
periodic  hexagonal  steel  array). 


The  relations  (C-44)  satisfy  the  j*-periodicity  condition  on  the  (approximate)  outer  cell 
boundary.  The  form  (C-44)  represents  the  first  relevant  terms  of  a  Taylor  expansion 
in  the  singly  connected  region  and  a  Laurent  expansion  in  the  annular  region 
A^.  The  quantities  n^.n^  denote,  respectively,  the  volume  fraction  of  the  steel 
and  concrete.  If  r^.r^  represent  the  current  steel  and  cell  radii,  respectively,  then 
nl1)  =  (r^/r^)2.  n^  =  1  -  n^.  The  dimension  2r^  represents  essentially  the 
steel  spacing. 

If  more  accurate  measures  of  the  micro-deformation  process  are  desired,  they  can 
be  generated  by  a  more  sophisticated  selection  of  weighting  functions.  A  different  set 
of  weighting  functions  would  also  be  necessary  if  the  cell  aspect  ratio  differed 
considerably  from  unity. 

C.9.6  Constraints. 


For  some  applications,  it  may  be  appropriate  to  introduce  the  constraint: 


<jK(l)  =  gK(2) 


(C— 45) 


In  accordance  with  the  reduction  of  the  number  of  dependent  variables  represented  by 
(C-45).  the  relation  (C-40)  is  replaced  by 


V«MK  -  NK  +  FK  =  0 

w  w  A# 


(C-46) 


where 

(mV.eI  ■  l(!KW/('»!<Wj  (C— 47) 

K  K 

Constitutive  equations  for  the  new  dependent  variables  M  .  are  obtained  by 
adding  (C-35) ,  for  a  =  1.2:  and  by  adding  (C-36)  for  a  =  1.2.  The  result  is 

>  =  f;  [akW:|<”>  .  .  cKL(“>.5Ll  , 

0=1  1  J 


(C-48) 


(C— 49) 


Sk 

r 


K(a).|(«) 


♦  feKL(*).>00  +  £KL(a).jL] 


C.9.7  Coordinate  Systems. 

The  mixture  relations  (C-47)  -  (C-52)  are  based  on  the  premise  that  the  current 
configuration  is  also  the  reference  configuration.  Otherwise,  the  equations  are  general. 
In  particular,  since  direct  notation  has  been  used,  the  mixture  relations  are  valid  for 
arbitrary  global  (macro)  coordinate  systems.  Care,  however,  must  be  exercised  in 
performing  coordinate  transformations  since  the  local  or  microcoordinates  represent  a 
preferred  or  “ materiaT'-type  system.  For  example,  the  operator  .  and  the 
constitutive  tensors-:#.  3?  and ^{J  are  all  referred  to  a  local  coordinate  system  with  base 
vectors  as  illustrated  in  Figure  C.4. 

Let  (  )*  denote  tensor  components  referred  to  the  preferred  coordinate  system. 
For  simplicity,  let  this  system  be  rectangular  Cartesian.  Then  (see  Figure  C.4).  with 
respect  to  any  other  (rectangular  Cartesian)  system,  one  finds 


Pi  ■  ^ik^i»e4«.ftj  '  VVL  5  1 


4t ds  ■ 


"i  ■ '  "ki  B  ; 

y 

B*<?>  =  to  [  9K<*>dA*  ; 

•  j k  i  jk£. 

JA(«) 


(C-50a) 

(C-50b) 

(C-50c) 


KL(a) 

Cijk 


f  K(a)n  K  (a)..* 


(C-50d) 


Coordinate  systems:  g*.  g.  (i  =  1,2,3)  denote  microcoordinates 
and  macrocoordinates  respectively. 


SCA-1679 


w  -  * 


(C-50f ) 


cKL(a)  K 

ij  r£mrkn  8  ,m  ,m  £.ijk 

JAw 


In  the  above.  denote  angle  cosines  defined  by 


d. .  =  e- •  e. 

r,  j  ~ ,  ~j 


(C-50g) 


(C-50h) 


where  £..  £.*  are  unit  orthogonal  base  vectors  associated  with  the  general 
macrocoordinate  system  and  the  microcoordinate  system,  respectively.  The 
components  of  all  other  tensors  in  (C-39)  -  (C-43)  do  not  involve  the 
microcoordinates  and  hence  can  be  written  directly  without  involvement  of 


C-23 


I, 

f 


APPENDIX  D 

MIXED  VARIATIONAL  PRINCIPLE  FOR  SYNTHESIZED  FIELDS 

The  direct  variational  method,  as  outlined  in  Appendix  C.  is  convenient  from  a 
theoretical  construction  standpoint.  However,  there  are  certain  disadvantages 
associated  with  such  an  approach.  For  example,  the  constitutive  relations  for  the 
interaction  terms  do  not  evolve  naturally.  Also,  the  microstresses  resulting  from  the 
direct  method  will  be  rather  crude  unless  complex  weighting  functions  are  utilized. 
Further,  the  system  of  equations  generated  by  the  direct  method  may  be  too  "stiff'. 

In  an  effort  to  overcome  problems  such  as  those  noted  above,  the  use  of  other 
variational  methods  was  explored.  One  such  method  —  a  Reissner-type  mixed 
variational  formulation  —  was  found  to  be  particularly  useful.  This  technique  is 
discussed  next. 


According  to  (C-25)  of  Appendix  C.  the  conservation  of  linear  momentum  is 
satisfied  in  an  average  sense  via  a  weighted  residual  procedure  wherein  the  variations 
in  velocity,  constitute  the  weighting  functions.  In  the  mixed  variational  method, 

both  the  conservation  of  linear  momentum  and  the  constitutive  relations  are  to  be 
satisfied  in  an  average  sense  via  the  following  weighted  residual  procedure: 
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According  to  (D-l).  both  velocity  and  stress  (increment)  are  to  be  varied  in  this 
principle.  Equation  (D-l)  is  (25)  of  Appendix  C  with  stress  in  incremental  form 
together  with  constraint  conditions  which  are  represented  by  the  second  and  third 
terms  of  (D-l).  In  the  second,  is  to  be  computed  from  the  constitutive  relation 
(C-30)  given  the  stress  increment,  i.e.,  formally  by 


-  M 
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where  Q  denotes  the  inverse  of  Q.  Similarly,  in  the  third  term,  [j^j  is  to  be 
computed  from  the  interface  constitutive  relation  (C-31).  i.e..  formally  by 
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The  variations  of  5%  are  to  satisfy  the  same  conditions  as  described  for 
(C-25).  Appendix  C.  The  variations  6g  are  arbitrary  over  V.  but  are  to  satisfy  the  j* 
-  periodicity  condition.  The  variation  £J*  is  arbitrary  over  the  interface 

With  application  of  Gauss’  Theorem,  use  of  (C-26)  -  (C-28)  together  with  the 
j5*-periodicity  of  Equation  (D-l)  can  be  written 
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Equation  (D-4)  is  an  incremental  form  of  the  Principle  of  Virtual  Work  for  the 
synthesized  fields.  The  terms  of  (D-4)  involving  and  6J*  are  constraint 

conditions  that  reflect  satisfaction  of  the  material  and  interface  constitutive  relations 
respectively:  here  and  5J*  play  the  roles  of  Lagrangian  multipliers. 

The  mixed  variational  method  has  been  used  to  generate  mixture  models  for 
several  important  cases.  This  work  is  presented  in  Appendices  E  -  H.  In  Appendix 
E.  a  model  is  derived  for  a  dense  ID  rebar  layout,  and  for  3D  deformation  of  the 
resulting  composite.  In  Appendix  F.  a  model  is  derived  for  a  R/C  plate:  this  is  a  2D 
problem.  In  Appendix  H.  the  basic  equations  for  a  R/C  rod  (ID  problem)  are 
obtained. 
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APPENDIX  E 


APPLICATION  OF  MIXED  VARIATIONAL  METHOD  TO 
DENSE  STEEL  LAYOUT  -  3-D  PROBLEM 

The  purpose  of  this  appendix  is  to  demonstrate  the  application  of  the  mixed 
variational  method  by  a  specific  example  calculation.  For  this  purpose,  a  dense 
unidirectional  steel  layout.  Figure  E.l.  will  be  considered.  For  simplicity,  the 
development  will  be  restricted  to  small  deformations.  Further,  the  concrete  will  be 
modeled  as  elastic-brittle  fracture,  and  the  steel  as  elastic.  Only  tangential  steel- 
concrete  slip  will  be  allowed;  otherwise,  the  interface  constitutive  relation  will  remain 
reasonably  general. 

If  the  rectangular  Cartesian  reference  system  illustrated  in  Figure  E.l  is  selected, 
a  quasi-static  process  is  considered  without  body  forces,  and  small  deformations  are 
assumed,  then  the  basic  equations.  (C-19)  -  (C-24)  of  Appendix  C.  of  the  synthesized 
fields  reduce  to: 


(a)  Equilibrium 
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In  the  above,  X^  and  are  constants,  whereas.  Kj.  is  a  (variable)  tangent 
modulus. 

Under  the  conditions  noted  above,  the  variational  principle.  (4)  of  Appendix  C. 
can  be  written  in  component  form  as  follows: 


In  the  above,  u.  is  the  displacement  vector.  (  )  =  9(  )/9t  and  (  ).*  =  9(  )/9x. 
with  9(  J/dXj  =  O.  In  what  follows,  all  stress  components  will  be  vaVied  with  the 
exception  of  a  ;  thus  =  0.  The  notation  e.!a'(»»«)  implies  inversion  of  the 
constitutive  relations  for  (i.j)  t  (1.1).  Similarly,  the  notation  [u. ](•••)  implies 
inversion  of  the  interface  constitutive  relation.  Herein,  only  t.will  be  varied:  thus  dt_ 
O.  '  2 

Appropriate  trial  displacement  and  transverse  stress  fields  for  use  in  (E-5)  can  be 
obtained  by  an  asymptotic  procedure  (Murakami.  Maewal,  and  Hegemier  [1981]). 
Accordingly,  one  invokes  the  premise  that  the  macrodimension  A  is  much  smaller  than 
the  microdimension  A.  i.e..  e  =  A/A<<  1.  This  premise  and  the  form  of  the  scaled 
equations  (E-l)  and  (E-2).  suggest  an  expansion  of  the  dependent  variables  in  the 
asymptotic  series: 
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If  (E-6)  is  substituted  into  (E-l)  -  (E-4).  and  the  coefficients  of  different  powers 
of  e  are  equated  to  zero,  a  sequence  of  problems  defined  on  the  cell  is  obtained. 
Consideration  of  this  sequence  to  0(c)  provides  the  following  form  for  the  trial  fields 
(see  Murakami  and  Hegemier  [1986]): 


(a)  Displacement  trial  functions 
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(b)  Stress  Trial  Functions 
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Following  substitution  of  the  trial  fields  (E-7)  and  (E-8)  into  the  variational 
principle  (E-5),  use  of  Cauchy's  formula  to  evaluate  T*according  to 
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and  appropriate  integrations  by  parts,  one  obtains  the  Euler  equations  of  the 
variational  principle  in  the  form 
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(a)  Mixture  Equilibrium  -  stress  averages 
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(b)  Mixture  Equilibrium  -  stress  moments 
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(c)  Constitutive  Relations  -  transverse  stress  averages 
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(d)  Constitutive  Relations  -  interaction  terms 
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(e)  Constitutive  Relations  -  stress  moments 
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where  it  is  understood  that  2M. .  =2M..  3M..  =3M.  . 
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Additional  definitions  used  in  the  foregoing  equations  include 
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(f)  Constitutive  Relations  -  axial  stresses 

The  remaining  constitutive  relations  associated  with  a  Jfl ^  are  obtained  from 
(E-2a)  with  i=j=l:  the  result  is 
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(g)  Boundary  Conditions 

The  boundary  conditions  on  9V  associated  with  Euler  equations  are 
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Equations  (E-10)  -  (E-23)  define  a  well-posed  boundary  value  problem  with 

respect  to  the  macrocoordinates  x, . 
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APPENDIX  F 


APPLICATION  OF  MIXED  VARIATIONAL  METHOD  TO 
R/C  PLATE  -  2D  PROBLEM 

The  purpose  of  this  appendix  is  to  furnish  a  second  example  application  of  the 
mixed  variational  method.  In  this  example,  mixture  equations  are  derived  for  the 
planar  deformations  of  a  R/C  plate.  Figure  F.l.  The  rebar  layout  is  unidirectional, 
small  deformations  will  be  considered.  The  concrete  is  modeled  as  elastic-brittle 
fracture,  and  the  steel  remains  elastic.  The  development  includes  tangential  steel- 
concrete  slip. 


As  in  the  first  example,  the  basic  equations  of  the  synthesized  field  (C- 19) - 
(C-24)  of  Appendix  C  are: 

(a)  Equilibrium 
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(b)  Constitutive  Relations 
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In  addition,  the  variational  principle  (4)  of  Appendix  D  can  be  written 
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Consider  now  a  preferred  rectangular  Cartesian  Coordinate  system  as  depicted  in 
Figure  F.l.  The  x.-axis  in  this  system  is  aligned  with  the  rebar  axis,  and  x3  is 
normal  to  the  plate. 

The  analysis  commences  by  invoking  the  following  premise  concerning  the  global 
stress  Fields: 
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where,  as  in  Appendix  C 


m!«)  h 
•j 


'J 

a(“) 


dA* 


(F-7) 


Equation  (F-6)  represents  a  type  of  global  "plane  stress". 


Next,  the  following  trial  displacement  field  is  selected: 
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where  r.!a^.  t.!a^.  and  P.  are  functions  of  x„.  x,  and  t. 
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Equations  (F-10)  are  a  special  case  of  (E-8),  Appendix  E.  for  which 
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Equations  (F-lla.b)  result  from  the  premise  (F-6)  whereas  (F-llc)  arises  from  the 
constraint 
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i.e.,  the  top  and  bottom  concrete  surfaces  are  assumed  to  be  traction  free. 


Substitution  of  the  trial  displacement  and  stress  fields  into  the  variational 
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principle  (F-5) .  noting  that  '=  0.  evaluating  T.  by 
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and  performing  necessary  integration  by  parts  furnishes  the  Euler  equations  of  the 
variational  problem  is  the  form 


(a)  Mixture  Equilibrium  -  Stress  Averages 
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(c)  Constitutive  Relations  -  Stress  Averages 
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(d)  Constitutive  Relations  -  Interaction  Terms 
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(e)  Constitutive  Relations  -  Stress  Moments 


APPENDIX  G 


MIXTURE  EQUATIONS  FOR  ROTATED  COORDINATE  SYSTEM 

The  coordinate  system  selected  in  Appendix  F  is  a  preferred  system.  For  many 
problems,  another  reference  system  may  be  desirable. 

Consider  Figure  G.l.  where  the  arbitrary  coordinate  system  is  related  to  the 
preferred  system  by 


x  =  Xj  cos0  -  Xj  sin0 
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Under  the  transformation  (G-l).  the  basic  mixture  relations  for  the  two-dimensional 
R/C  plate  take  the  form 
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(b)  Mixture  Equilibrium  -  stress  moments 
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(c)  Constitutive  Relations  -  stress  average 
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(d)  Constitutive  Relations  -  stress  moments 
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(d)  Constitutive  Relations  -  interaction  terms 
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APPENDIX  H 


APPLICATION  OF  MIXED  VARIATIONAL  METHOD  TO 
R/C  ROD  -  ID  PROBLEM 

In  this  appendix,  the  governing  mixture  relations  for  a  R/C  Rod.  Figure  1.  are 
presented.  This  geometry  is  relevant  since  it  represents  specimens  that  are  most 
often  used  to  investigate  steel-concrete  bond  properties  in  the  laboratory.  Tension 
tests  are  also  frequently  performed  using  this  geometry. 


The  mixture  equations  for  the  R/C  rod  can  be  obtained  from  the  analysis  of 
Appendix  E  by  imposing  a  number  of  constraints.  With  reference  to  Equation  (8). 
Appendix  E.  the  appropriate  displacement  constraints  are 


while  the  necessary  stress  constraints  (see  Equation  10  of  Appendix  E)  are 
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Therefore,  the  trial  displacement  field  becomes 
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whereas  the  trial  stress  field  becomes 
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In  view  of  (H-4a).  the  appropriate  constitutive  relations  are 
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where  E^. 


denote  Young's  modulus  and  the  shear  modulus,  respectively. 


The  desired  mixture  relations  can  be  obtained  by  substituting  (H-3)-(H-5)  into 
the  mixed  variational  principle  of  Appendix  0.  This  same  result  can  be  obtained  by 
invoking  (H-l)  and  (H-2)  in  the  R/C  plate  equations  of  Appendix  F,  and  by  noting 
that  (  )2=0.  The  result  is 

(a)  Mixture  Equilibrium 
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In  the  above.  =  n^  Nj|ap^  denotes  a  "partial"  stress. 

The  foregoing  relations  apply  prior  to  local  bond  slip  and  prior  to  plastic 
deformation  of  the  rebar.  Where  the  latter  nonlinearities  occur,  the  quantities  p^  and 
must  be  interpreted  as  tangent  moduli  as  in  Appendix  C. 
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APPENDIX  I 


ON  DOWEL  ACTION  IN  REINFORCED  CONCRETE 
1.1  INTRODUCTION. 

The  nonlinear  response  of  reinforced  concrete  subsequent  to  initial  cracking  is 
often  dominated  by  complex  mechanisms  associated  with  steel-concrete  interactions. 
These  mechanisms  include  aggregate  interlock  or  interface  shear  transfer  (1ST),  dowel 
action  (DA),  nonlinear  dowel  effects  (NDE).  and  steel-concrete  bond  action  (BA). 
They  are.  in  general,  interdependent.  This  interdependence  places  a  severe  limitation 
on  experimental  attempts  to  quantify  the  contribution  of  each  mechanism. 

The  foregoing  experimental  difficulties  are  amplified  by  a  dependence  of  test  data 
on  the  test  set-up  For  example.  Fenwick  and  Paulay  [1]  have  shown  that  the 
measured  global  force  versus  displacement  relations  from  DA  and  DA  plus  1ST  tests 
depend  strongly  on  the  geometry  of  the  test  specimen  and  the  test  set-up.  i.e..  on 
structural  constraints.  The  inability  to  separate  structural  effects  and  basic  material 
properties  within  a  framework  of  continuum  mechanics  renders  the  development  of 
empirical  relations  for  the  response  of  reinforced  concrete  members  of  arbitrary 
geometry  and  boundary  conditions  intractable. 

Problems  such  as  those  noted  above  demonstrate  the  need  for  analytical 
descriptions  of  the  basic  interaction  mechanisms.  Even  one-dimensional  analytical 
models  (Friberg  [2];  Rehm  [3];  Hegemier.  et  at.  [4])  have  proved  to  be  very  useful  in 
studies  of  certain  coupled  interaction  mechanisms  in  reinforced  concrete  members. 

To  achieve  optimum  utility  and  simulation  capability,  the  analytical  descriptions 
noted  above  should  be  nonphenomenological.  In  particular,  the  global  behavior  of  a 
reinforced  concrete  member  should  be  synthesized  from  the  properties  and  geometries 
of  the  composite  (steel,  concrete,  interfaces).  Further,  the  description  should  provide 
measures  of  local  fields  (e  g.,  average  stress  states  in  the  steel  and  concrete,  stress 
states  at  interfaces,  interface  slip,  etc.)  as  well  as  global  fields 


The  construction  of  an  analytical  model  of  reinforced  concrete  that  can 
synthesize  global  behavior  from  component  responses  is  a  difficult  task  which  requires 
some  level  of  micromechanical  analysis.  This  is  especially  true  when  a  DA-type 
interaction  mechanism  is  active.  In  this  paper  one  procedure  for  accomplishing  this 
task  is  discussed.  For  this  purpose  a  specific  problem  is  selected.  The  latter,  and 
the  developed  model,  involve  combined  DA,  NDE.  1ST  and  BA.  This  work  represents 
a  generalization  of  a  linear  dowel  action  model  (Murakami  and  Hegemier  [5])  to 
include  large  deflections,  moderately  large  rotations,  and  small-strain  plastic 
deformation  of  the  rebars.  As  a  result,  and  in  contrast  to  the  previous  work,  the 

developed  analytical  description  is  capable  of  simulating  large  slip  dowel  tests 

(Karagozian  and  Case  [6];  Paulay.  Park  and  Phillips  [7]). 

Following  model  development,  the  simulation  capability  of  the  theory  is 
demonstrated  for  both  monotonic  and  hysteretic  dowel-deformation  time  histories.  The 
monotonic  cases  are  conducted  with  and  without  1ST.  The  hysteretic  example 
excludes  1ST.  BA  and  NDE  effects  are  included  in  all  examples. 

The  discussion  herein  focuses  on  a  special  case  of  a  more  general  problem.  In 
particular,  the  crack  plane  is  assumed  to  be  normal  to  the  principal  rebar  direction. 
The  general  case  in  which  concrete  crack  planes  occur  at  arbitrary  angles  with  respect 
to  the  principal  steel  directions  will  be  treated  in  a  subsequent  work. 

1.2  FORMULATION. 

Consider  a  reinforced  concrete  element  of  length  2t  with  a  pre-existing  crack  (in 
the  concrete)  normal  to  the  rebar  at  the  center  of  the  undeformed  specimen.  (This 
element  may  represent  an  actual  structure  or  a  "cell"  in  a  periodic  undirectional  steel 
layout,  see  Murakami  and  Hegemier  [5]).  Let  rectangular  Cartesian  coordinates  Xj. 
Xj,  Xj  be  selected  with  X}  in  the  axial  (rebar)  direction  as  shown  in  Figure  1.1.  The 

crack  surface  in  the  undeformed  configuration  is  thus  located  at  Xj  =  0.  For 

notational  convenience.  (  )^.  a  =  1.2.  will  denote  quantities  associated  with  material 
a  with  a  =  1  representing  the  rebar  and  a  =  2  the  concrete.  In  the  cross-sectional 
Xj.  Xj-plane  the  rebar  center  is  located  at  the  origin  of  the  coordinates,  and  material 


a  occupies  the  domain  A^;  the  exterior  boundary  of  U  A^  and  the  interface 
between  the  steel  and  concrete  are  denoted  by  9A  and  /.  respectively. 

In  what  follows,  a  dowel-type  problem  will  be  considered  for  which  the 
displacements  on  9A  are  prescribed  as  depicted  in  Figure  1.1.  For  a  given  history,  the 
corresponding  global  forces  in  the  coordinate  directions  are  sought,  as  well  as  details 
concerning  the  stress  states  in  material  a  which  relate  to  failure. 

Model  Construction  for  DA 

A  model  for  deformation  in  the  X^.X^plane  is  developed  by  application  of  the 
direct  variational  method  in  conjunction  with  an  appropriate  trial  displacement  field, 
and  by  invoking  an  approximation  of  plane  stress  in  the  X^X^plane  at  a  certain 
point  of  the  analysis.  Both  the  trial  functions  and  the  plane  stress  assumption  have 
been  previously  adopted  in  the  derivations  of  classical  beam-type  theories. 

Analysis  of  the  dowel  problem  commences  by  application  of  the  principal  of 
virtual  work  to  each  material  in  the  undeformed  configuration.  Addition  of  the  results 
furnishes 


where  s!^  is  the  second  Piola-Kirchhoff  stress.  E is  the  Green-Lagrange  strain,  u 
is  the  displacement,  [u.]  =  uj2^  -  uj1^  is  the  displacement  discontinuity  (slip)  at  the 
surface  /.  T.  is  the  stress  vector  defined  on  the  surface  with  exterior  normal  u ,  and 

•  I  •  | 

T.  denotes  the  stress  vector  on  I;  both  T  and  v (  are  defined  in  the  undeformed 
configuration  In  (1-1)  the  virtual  displacement  and  strain  are  expressed  by  5u  and 


5 E..,  respectively,  and  ds-  denotes  an  infinitesimal  line  element  on  9AT  and  /;  the 
ij  0  T 

quantity  9AT  denotes  the  subdomain  of  9A  on  which  the  traction  is  specified. 

In  the  development,  it  is  assumed  that  the  rebar  experiences  large  deflections 
and  moderate  rotations  with  small  strain  plastic  deformation.  At  the  interface  /.  the 
rebar  is  allowed  to  suffer  axial  slip  [uj;  the  concrete  in  the  neighborhood  of  /  may 
experience  cracking  or  crushing  due  to  the  relative  transverse  displacement  [ u 2 ] 
between  the  rebar  and  concrete.  The  interface  /  is  an  "effective"  surface  associated 
with  a  bar  that  is  actually  deformed. 

The  assumed  trial  displacement  field  for  use  in  (1-1)  has  the  form 

“ia)  [Ws)  =u'a)(xi)  (XJXj  .  (1-2) 

where  6..  is  the  Kronecker  delta. 

1 1 

Under  (1-2),  the  Green-Lagrange  strains  become 
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For  simplicity,  and  consistent  with  the  premise  of  moderate  rotations,  the  von  Karmon 
approximation  [8]  has  been  adopted  in  which  only  the  nonlinear  terms  associated  with 
U2  j  have  been  retained  in  (1-3). 


Substitution  of  (1-2)  and  (1-3)  into  (1-1)  and  integration  by  parts  furnishes 
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and  where  (  )*  implies  that  (  )  is  evaluated  on  the  interface  I  of  the  undeformed 
configuration. 


In  the  derivation  of  (1-4).  the  following  expressions  were  employed: 
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where  T..  is  the  first  Piola-Kirchhoff  stress  and  =  S..u.  is  the  stress  vector 
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related  to  the  second  Piola-Kirchhoff  stress.  In  the  spirit  of  the  von  Karmon 


approximation  for  moderate  rotations,  the  undeformed  coordinates  X  were  related  to 


the  deformed  coordinates  x.  by 


dXj/dXj  =  1  .  9V3X2  =  0 
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Equation  (1-4)  indicates  that  the  Euler-Lagrange  relations  of  the  variational 
problem  on  0  <  X}  <  £.  are: 
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where  (l-7a)  was  used  to  simplify  ( I- 7b) .  In  addition,  the  appropriate  boundary 
conditions  to  (1-7)  at  X}  =  0  and  X^  =  ±  t  ate 
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either  5uja^  =  0  or  =  T^aa^  , 

either  <5U<a)  =  0  or  [n^  ♦  N<J>]  =  T^a)  (I-8b) 

either  5j^a^  =  0  or  =  j^(aa)  (I-8c) 

Equations  (1-7)  are  of  the  mixture  form.  The  quantities  P  ,  Pj.  P^  are 
interaction  terms,  which  reflect  stress  transfer  between  the  rebar  and  concrete  across 
the  interface  /.  The  quantities  Q  .  Q2.  denote  forces  or  moments  resulting  from 
tractions  on  the  exterior  cell  boundary  3A^.. 

In  the  remainder  of  the  analysis,  attention  is  focused  on  a  special  dowel  problem 
for  which 

U[2)  =  0.  U<2)  =  sgn  (Xj)  A/2.  *[2)  =  0  .  (1-9) 

where  sgnfXj)  =1  for  Xj  >  0  and  =  -1  for  Xj  <  0. 


This  case  corresponds  to  a  uniform  displacement  of  the  exterior  boundary  9A  of  the 
concrete  cover  (see  Figure  1.1).  Under  (1-9)  6u!2^  =  0  on  0  <  X}  <  £.,  and  the 
relevant  mixture  relations  reduce  to 
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For  the  case  of  DA  only,  a  lubricated  interface  is  assumed  at  X  =  0.  Further, 
a  test  set-up  is  considered  for  which  the  concrete  is  traction  free  at  trie  termini  Xj  = 
+  £.  This  condition  corresponds  to 
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at  Xj  = 
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Tj2a)  =  0  at  Xj  =  ♦  l  .  (I-12b) 

In  addition  to  (1-10),  constitutive  relations  are  required  for  N|j^.  n|^.  m|^.  Pj. 
P  pj*  to  close  the  model.  For  this  purpose  it  is  necessary  to  relate  tne  time  rate 
of  change  of  the  second  Piola-Kirchhoff  stress  §  and  the  time  rate  of  change  of  the 
Green-Lagrange  strain  £  =  9||/9t  at  time  t  to  the  corresponding  quantities  in  the 
deformed  configuration.  For  small  strain  problems  it  is  possible  to  interpret  the 
components  of  §  as  those  of  the  Cauchy  stress  £  with  respect  to  the  covariant 
bases,  and  the  components  of  |  as  those  of  the  rate  of  deformation  tensor  ().  with 
respect  to  the  convected  contravariant  bases  in  the  deformed  configuration.  The 
premise  of  small  strain  implies  that  the  convected  covariant  bases  (gj.g2)  can  be 
approximated  by  the  rotated  rectangular  Cartesian  bases  (g*jg*2)  at  ^  =  0.  one  of 
which  is  gj  and  the  other  orthogonal  to  g1  to  form  a  right-handed  system.  Figure  1.2. 
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As  a  result  of  these  approximations,  one  may  interpret  N(!^.  Mjj^.  Pf  and  P^1  as  the 
corresponding  quantities  defined  in  terms  of  the  Cauchy  stress  with  respect  to  rotated 
rectangular  Cartesian  bases. 

For  elastic  response,  the  appropriate  stress  constitutive  relations  are.  under  the 
approximation  of  plane  stress  in  the  X^.  X  plane. 


§(!) 

511 

_  E(l)  t(l) 

'  L11  * 

(I-13a) 

§(!) 

b12 

=  *<l>  6g>  - 

(I-13b) 

Substitution  of  (1-13)  into  (l-5a)  with  use  of  (1-3)  for  a  =  1  furnishes  the  elastic 
segment  of  the  constitutive  relations  in  the  form: 


■S’ ,  -  w  ««  .  «<«  oo>)  , 


(I-14a) 


-  n(1>A(1)A(0^  ♦  »«)  • 


(I-14b) 


=  E(1)  I(1)  (I— 14c) 

where  n^A  is  the  area  of  the  subdomain  A^  and  n^  denotes  the  area  fraction  of 
the  rebar  in  the  X2,  Xj-plane,  A  is  the  total  area  of  the  slip  plane,  and 


i(1)  E  Wad  py2  dx2  dx3  • 


(1-15) 


For  plastic  response  of  the  rebar,  the  von  Mises  yield  criterion  with  the 
associated  flow  rule  is  employed  with  respect  to  the  rotated  coordinate  bases  as 
follows: 


f(2)  =  J3^  -  *y00  =  0 


(1-16) 


where,  for  plane  stress, 


(1-17) 


and  o  is  the  yield  stress  in  simple  tension:  the  latter  is  a  function  of  the  plastic 
work  k  where 


‘  =  (1-18)  | 

The  quantity  Dp  is  the  plastic  part  of  the  rate  of  deformation  tensor.  Application  of  1 

the  associative  How  law  gives 


,p  _  _3X 

12  "  2a  °12 
y 


(1-19) 


where  X  is  the  plastic  multiplier.  By  invoking  the  consistency  condition.  X  can  be 
evaluated.  By  further  invoking  the  premise  of  small  strains.  £  and  p  can  be 
interchanged  with  §  and  £.  respectively,  as  was  noted  previously.  Following  use  of 
(1-13)  for  the  elastic  part  of  Q.  i.e..  p*  =  Q  -  Qp,  one  obtains  the  following  local 
constitutive  relations: 


where 


> 

> 


(1-21) 


Equation  (1-20)  can  now  be  employed  to  furnish  constitutive  relations  for  the 
global  quantities  N^Jj.  and  ^  by  substitution  into  (l-5a).  To  facilitate  the 

resulting  integrations,  and  in  order  to  arrive  at  a  result  of  practical  utility,  three 
approximations  are  made  prior  to  integration.  The  first  consists  of  dropping  terms  of 
0(Sj2)  compared  to  unity  in  the  diagonal  elements  and  denominator  of  (1-20):  the 
second  consists  of  approximating  a2  by  in  (I-2U).  i.e..  assuming  S*2  <<  sir  lhe 
third  involves  the  assumption  tha^  the  cross-section  is  fully  plastic.  The  result  of 
these  approximations,  which  are  judged  to  be  reasonable  for  the  problem  class  under 
consideration,  is 
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where 


c  =  E(1)H  .(1) 
'•1  1 


n  -  E(i).  H"  A  •  ci2  *  ,  (Ed).  H) 
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(1-23) 


The  quantity  E^H/E^  +  H  in  (1-23)  is  the  elasto-plastic  tangent  modulus.  The 
relation  (1-22)  was  obtained  with  the  aid  of  a  quadratic  interpolation  of  SJ2  over  the 
cross-section  A^.  i.e. 


(1-24) 


cd) 

i2 


=  ^(a2-X^]  ,  1.  [a2  -  Xl) 


X2JdX2dX3 


(1) 


where  a  is  the  radius  of  the  domain  A^. 


In  addition  to  the  foregoing  global  stress  constitutive  relations,  expressions  must 

be  considered  that  define  the  constitutive  behavior  of  the  interaction  terms  P.(i  = 

- 

1.2).  Here  P..  eg.,  in  (l-10a)  and  (l-10b).  represent  the  forces  exerted  on  the  rebar 
by  the  concrete.  In  general,  one  may  write 

P;  =  Fj ([u1]]  (1-25) 

where  F.  denote  functionals  and 

i 

H  *  h]  •  1“2]  *  u2,iK)  *  lu2  ]  •  t1-26') 

For  the  present  problem 

[U.]  =  <5i2  sgn(X2]A/2  -  uf!)  (I-26b  ) 


where  6..  is  the  Kroncher  delta.  In  (1-26).  [U.]  is  the  interface  displacement 
discontinuity  defined  in  the  undeformed  configuration;  the  transformation  (l-26a)  has 
been  introduced  such  that  [u1]  are  components  with  respect  to  the  rotated  coordinate 
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ii 


yiViViy 


bases  in  the  deformed  configuration. 


The  interaction  term  P  in  (l-10a)  reflects  "bond"  action  between  the  steel  and 
concrete,  and  represents  a  tangential  force  on  the  steel  due  to  the  concrete.  In 
previous  studies  by  Hegemier  et  at.  [4]  and  Hageman  et  at.  [9]  it  was  found  that  this 
interaction  term  can  be  adequately  expressed  in  the  form 


P^K^u1]  (1-27) 

where  Kj  denotes  a  (non-constant)  tangent  modulus.  For  monotonic  deformation,  the 
following  bilinear  form  of  (1-27)  has  been  found  to  be  adequate: 

*1  ■  for  IPjl  l  (P1}y  . 

(I-28a) 

\  for  IPjl  >  (P1)y 


where  (Pj)  denotes  a  critical  value  at  which  the  modulus  changes  from  p^  to  p^ . 
If  one  introduces  the  hardening  parameter  H^,  then  p  **  can  be  written  as 


The  constant  (3  in  (1-28).  which  corresponds  to  no  bond  slip  (perfect  bond),  has  been 
evaluated  analytically  in  a  previous  work  by  Murakami  and  Hegemier  (5).  For  dowel 
problems  involving  hysteretic  bond  slip,  a  more  complex  representation  of  is 
necessary.  The  form  used  for  computational  purposes  herein  is  depicted  in  Figure  1.3. 


The  interaction  term  P 2  reflects  a  normal  force  on  the  steel  due  to  the  concrete. 
Again,  for  monotonic  deformation  the  bilinear  expression  below  has  been  adopted 


P2  =  P2W]  for  IPjl  i  (P2)y 

P2  =  />'P  It2]  for  IP2I  >  (PJ, 


(I-28c) 


where 


=  h 


Pc 


Pi  *  «, 


(I-28d) 


and  where  /L  has  been  evaluated  analytically  in  Murakami  and  Hegemier  [5],  For 
hysteretic  deformation,  the  form  of  K2  depicted  in  Figure  1.4  has  been  found  to  be 
adequate:  in  general,  the  incremental  form 


P2  =  K2Iu2] 


(I-28e) 


where  K2  is  a  tangent  modulus,  should  suffice  for  most  applications. 

It  is  noted,  that  in  contrast  to  P  where  the  nonlinear  part  of  the  bond  behavior 
must  be  postulated,  the  nonlinear  part  of  P2  should  be  derivable  from  the  component 
constitutive  behavior  and  geometry.  This,  however,  requires  the  introduction  of  a 
nonlinear  model  of  the  concrete  and  the  conduct  of  a  detailed  3-D  microanalysis. 


I 

) 


This  is  a  formidable  task  which  is  outside  the  scope  of  the  present  paper.  For 
purposes  of  examining  the  general  features  of  the  current  modeling  process,  the 
constitutive  relation  for  has  been  postulated  beyond  the  linear  range. 

Finally,  the  interaction  term  must  be  considered.  From  (l-5b)  it  is  evident 
that  this  quantity  represents  a  moment  generated  by  the  steel-concrete  interface  shear 
stress  S2j.  In  general,  this  moment  can  be  shown  to  be  negligible  for  problems  of 
practical  interest  which  necessitate  bond  breakage  prior  to  independent  bending  action 
of  the  rebar.  As  a  consequence  of  this  point,  it  is  assumed  that 

pJ1  =  0  .  (1-29) 


The  formulation  of  the  model  is  now  complete  for  DA  and  IMDE  without  1ST. 
The  basic  equations  are  (1-10),  (1-14).  or  (1-22).  (1-28)  and  (1-29).  and  the  boundary 
conditions  (1-8)  with  a  =  1.  For  convenience,  these  relations  are  summarized  below; 
the  interaction  terms  listed  refer  to  the  case  of  monotonic  deformation: 

Summary  of  Basic  Equations  for  DA 


(1)  Equilibrium: 


(2)  Stress  Constitutive  Relations: 
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where 

C  -  E^n^A  C  -  0  c  tfd)n(l)A  C  =  gd)jd) 
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for  f  <  0  or  f  =  0  and  ?  <  0 
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f  .  TT£  -  »  ,  3J2  =  N<J> 
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for  f  =  0  and  ?  =  0 


n(1)A  ,  C12  =  -3 
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(3)  Interaction  Constitutive  Relations: 

Pj  =  p.  [u 1 ]  for  IP. I  <  (P.) 

=  /?®P[u’]  for  IP.  I  >  ^P.j  ,  (i  =  1,2,  no  sum  on  i) 


where 
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(I-31a) 


(I-31b) 


(I-31c) 


(I-32a) 


C  •/;  (‘  -  . 


(I-32b) 


[u1]  E  [uj,  [u2]  =  [u2]  .  U2il[Uj] 


(I-32c) 


(4)  Boundary  Conditions: 


-  o 
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0  ,  =0  at  Xj 


(1-33) 


(5)  "Driving"  terms: 


[Ui]  =  *i2  *9"  (X2]  A'2  -  Ui 


(1-34) 


Addition  of  1ST 

An  elementary  procedure  for  inclusion  of  1ST  in  the  foregoing  dowel  problem  is 
now  considered  as  follows:  First.  (1-30)  -  (1-34)  are  solved  incrementally  for  LlJ1^  and 
T.(la)  at  X.  =  0  where  l!la)  =  N1J1).  +jla)  =  N  +  U,  .(1)N.  .(1).  Second,  an 
effective  normal  stress  p  =  pQ  +  (T |la v  A)  is  defined  on  the  crack  face  at  Xj  =  0. 
where  pQ  is  a  confining  pressure  (see  Figure  1.1)  which  is  part  of  the  test  set-up  for 
combined  DA  -I-  1ST  experiments.  Third,  the  effective  normal  stress  p  is  applied  to 
the  concrete  interface  at  X.  =  0  and  an  1ST  law  is  postulated  in  the  form 

i  (2a\  *  i  (2al 

f(T'  \p.A)  =  0  and  this  relation  is  solved  for  T;  1  for  a  specified  p  and  A.  This 


pecified  p  and  A.  This 


latter  quantity  reflects  1ST.  According  to  this  procedure,  1ST  influences  the  total 
effective  shear  stress,  r.  across  X.  =  0  but  does  not  alter  the  DA  contribution. 
T<la>  Here 
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APPENDIX  J 


DESCRIPTION  OF  CONCRETE  UNDER  COMPRESSIVE 

STRESS  STATES 

J.l  INTRODUCTION. 

It  is  well  known  that  the  strength  of  concrete  increases  with  increasing  mean 
confinement  stress.  This  effect  can  be  observed  in  the  triaxial  failure  surface  of 
concrete  in  stress  space.  Figure  J.l. 

In  addition  to  strength,  the  ductility  of  concrete  is  observed  to  increase  with 
increasing  mean  compressive  stress.  Typical  increases  in  ductility  during  standard 
triaxial  tests  are  shown  in  Figure  J.2.  The  effect  of  confinement  on  ductility  is  seen 
to  be  much  greater  than  that  on  strength.  This  can  also  be  deduced  from  the 
behavior  of  reinforced  concrete  in  the  presence  of  confinement  steel.  Figure  J.3. 

Analytical  and  experimental  studies  (see  Nemat-Nasser.  et  at..  [1])  on  brittle 
materials  such  as  concrete  suggest  that,  at  low  confining  stress  levels  failure  occurs 
by  the  growth  and  connection  of  a  narrow  zone  of  microcracks  to  form  one  or  at 
most  several  macrocracks  or  macrofaults.  Figure  J.4a.b.  The  onset  of  this  event 
(formation  of  a  macrocrack)  corresponds  to  the  onset  of  "strain  softening"  [2].  In 
contrast,  at  higher  levels  of  confining  stress  microcracks  grow  in  a  more  uniform 
manner  and  a  distribution  of  microcracks  evidently  develops.  Figure  J.4c.  Although 
failure  may  again  occur  via  localization,  the  formation  of  this  distribution  appears  to 
be  responsible  for  the  observed  increase  in  ductility. 

Below  the  brittle-ductile  transition,  slip  along  the  planes  that  constitute  the 
above  distribution  is  a  more  plausible  explanation  of  increased  ductility  than  is  true 
plastic  flow.  Consequently,  conventional  phenomenological  elastoplastic  descriptions  of 
concrete  under  moderate  compressive  stress  states  appear  to  be  inappropriate. 

In  view  of  the  above  discussion,  an  effort  was  made  to  construct  a  theory  of 
concrete  based  on  a  "slip  system"  concept.  A  description  of  progress  made  to-date 
on  this  subject  follows. 


J.2  THE  SUP  SYSTEM. 


Slip  system  models  have  been  used  previously  to  describe  both  metals  and 
frictional  materials.  A  representative  cross-section  of  such  models  includes  the  works 
of  Mandel  [3],  Spencer  [4].  Asaro  [5].  and  Nemat-Nasser,  et  at..  (6],  An  excellent 
discussion  and  in-depth  treatment  of  the  subject  is  given  by  Nemat-IMasser  [7], 

Figure  J.5  describes  the  slip  concept  in  the  case  of  a  single  slip  system.  Locally 
n  represents  a  unit  normal  to  the  slip  surface  and  s  is  a  unit  tangent  vector  which 
describes  the  direction  of  slip.  The  quantity  A  denotes  an  effective  slip  surface 
spacing.  The  slip  system  and  its  geometry  is  presumed  to  form  at  a  certain  critical 

stress  state.  The  collection  of  these  states  forms  a  surface  in  stress  space. 

In  the  sequel,  a  theory  is  described  for  an  arbitrary  number  of  slip  systems. 

Associated  with  the  ath  system  is  a  normal  vector  n^.  tangent  vector  and 

spacing  A^.  Figure  J.5.  where  a  =  1.2  ....  N. 

According  to  the  slip  concept.  aH  irreversible  deformation  (i.e..  "plastic  slip")  is 
assumed  to  take  place  along  the  slip  surfaces.  The  material  between  these  surfaces 
will  be  modeled  as  elastic. 

J.3  DEFORMATION. 


Let  the  relative  rate  of  slip  across  a  slip  surface.  Figure  J.5.  be  denoted  as 
[vjsj].  Further,  let  g.  (i  =  1  to  3)  denote  the  unit  orthogonal  base  vectors  of  a 
rectangular  Cartesian  reference  system.  Then,  if  the  rate  of  deformation  tensor 


(J-l) 


where  v  is  the  velocity  vector,  is  decomposed  into  elastic  and  inelastic  parts  according 
to 


$ 


d. .  =  de.  ♦  dP. 

'J  ij  U 


one  obtains 


(J-2) 


hMl 


dP.  .1  fs.n.  ♦  s.n.U-* 

ij  2  l  i  j  j  iJ  A 

where 

s-  =  s  •  e.  .  n.  =  n  •  e. 

|  ~  N  |  9  |  ~  *•  | 

In  the  case  of  multiple  slip  systems,  the  effects  are  superposed  to  give 
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Since  s»n  =  s.n.  =  0,  it  is  evident  that  the  inelastic  deformation  characterized 
by  (J-5)  does  not  involve  volume  change.  There  are  three  plausible  ways  that 
inelastic  volume  change  may  accompany  plastic  flow  within  the  context  of  a  slip 
system  (Nemat-Nasser[7]):  (a)  plastic  volumetric  expansion  may  occur  normal  to  the 
slip  plane  as  slip  takes  place;  (b)  uniform  isotropic  expansion  or  contraction  may 
accompany  slip;  and  (c)  uniform  inelastic  volumetric  deformation  may  take  place 
independently  of  slip. 

The  physical  cause  of  (a)  is  asperities  in  the  form  of  aggregate  interlock.  Let 
[v.  .]  denote  the  rate  of  separation  of  the  +  slip  surfaces  in  the  normal  direction. 
Then 
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represents  the  dilatancy  due  to  (v(nj]  ^  l^e  latter  is  related  to  the  slip  [v(sjl 
according  to 


[V  (n)  ]  =  Mv(s)] 


(J-6) 


J-3 


where  6  is  a  scalar  function,  then  (J-5)  can  be  generalized  to  give 
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The  physical  basis  for  (b)  is  the  formation  of  microvoids  and  microcracks  due  to 
slip,  while  that  for  (c)  is  the  collapse  of  voids  under  hydrostatic  compression.  These 
effects  can  be  modeled  by  adding  to  (J-7)  terms  of  the  form 
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Consequently,  the  inelastic  deformation  rate  tensor  becomes 
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In  addition  to  the  rate  of  deformation  tensor  dp.,  a  spin  tensor  «P.  is 

IJ  IJ 

defined  according  to  (Asaro  |5]) 
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where 
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J.4  INTERFACE  BEHAVIOR. 


The  ath  slip  system  is  assumed  to  be  activated  if  the  resultant  shear  stress  ra 
has  attained  a  critical  value  defined  by 


I 


f a[ra ,oa ,p,eaP)  =  0  ;  a,p  =  1  to  N  (J-10) 

where  aa  denotes  the  normal  stress  on  the  slip  surface,  p  is  the  mean  pressure,  and 
eaP  are  measures  of  the  history  of  slip.  Slip  occurs  if.  in  addition  to  (J-10). 

=  0  .  (j-11) 


Equation  (J-11)  can  be  expressed  in  the  form 


•  a  *a  • 

r  ♦  a  tan r/^  ♦  p  tan ij^ 


N 

-  z 

p= 1 


aptisl 0 


(J-12) 


where  tan  r)^  and  tan  rj 2  are  material  parameters  that  characterize,  respectively,  the 
eiiects  of  normal  stress  and  hydrostatic  stress  on  slip.  Here  rj.  corresponds  to  the 
usual  friction.  The  quantity  r;2  reflects  the  influence  of  pressure  on  plastic  flow:  e  g., 
in  crystal  plasticity  this  term  would  represent  the  influence  of  pressure  on  the  motion 
of  dislocations  2.  (In  two-dimensional  problems  (J-12)  can  be  re-written  in  terms  of 
aa  only.)  The  parameters  haP  characterize  the  manner  by  which  the  slip  system's 
resistance  changes  due  to  slip.  In  what  follows  ha^  is  assumed  to  be  symmetric  (but 
not  necessarily  positive  definite).  Finally,  the  effect  of  cohesion  has  been  excluded  in 
the  above  formulation. 

Consider  now  the  problem  of  calculating  the  time  rate  of  change  of  ra  and  oa 
in  (J-10)  -  (J-12).  Since 


a  a  a 

s  =  a.  -n. s. 
*  ij  i  j 


„  a 

a .  . n. n . 

'J  '  J 


(J-13) 


it  is  clear  that  one  must  determine  and  in  order  to  calculate  ra  and  aa .  The 
former  rates  must  be  objective.  Based  on  the  physical  processes  under  consideration, 
the  following  rates  of  change  of  the  unit  vectors  q0  and  are  defined 


J-  5 


(J-14) 
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n.  =  w .  .  n  .  .  s.  —  w .  ■  s  • 
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According  to  (J-14).  the  unit  vectors  are  corotational  with  the  elastic  deformation  of 
the  slip  system.  With  use  of  (J-14).  equations  (J-13)  furnish 


•a  V  a  a  »a  V  a  a 
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where 
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If  one  defines  the  quantity  q..  according  to 


i  .  1  [  a  i  a  a)  a  a  .  e 

qij  =  5  IVj  *  *j"ij  ♦  Vj  *•»  ii  *  5ij  *■" 
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then  (J-12)  furnishes  the  following  relation  which  can  be  used  to  calculate  $ij: 
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J.5  GLOBAL  CONSTITUTIVE  RELATIONS. 

In  this  section  the  global  constitutive  relation  for  a.,  is  assembled.  For  this 

|j 

purpose  the  elastic  part  of  the  slip  system  deformation  is  postulated  in  the  form 


l..  =  L.  ..0d® 

i  j  i  j  k£.  kfc. 


(J-19) 


where  L..^  is  the  instantaneous  elastic  tensor  modulus.  If  the  material  between  slip 
planes  is  isotropic,  then 


Li j kt  =  G(5ik5j£.  4  5i£.5jk)  4  X<5i  j5k£ 


(J-20) 
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It  will  be  assumed  that  L..^  is  independent  of  the  rate  of  elastic  deformation. 

With  use  of  (J-2).  (J-8a).  and  ( J-18) .  the  relation  ( J- 1 9)  can  be  cast  in  the 
following  form: 


■■  ~  Cijk£  (dk£  ‘  6kl  6 ) 
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where  C*.k£  is  an  "elastic-plastic"  modulus  defined  by 


(J-21) 
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where 

M*''5  •  (J'**) 

J.6  REMARKS. 

The  constitutive  relations  (J-21)  are  fully  nonlinear  and  incorporate  an  arbitrary 
number  of  slip  systems.  Their  form  is  identical,  for  all  practical  purposes,  to  the 
relations  obtained  by  Nemat-Nasser  [7] .  The  only  relevant  difference  concerns  the 
interpretation  of  the  tensor  L..^.  The  definition  of  this  quantity,  which  may  depend 
on  the  current  stress  and  the  history  of  deformation,  is  somewhat  vague  for  brittle 
materials  such  as  concrete. 


The  constitutive  relations  (J-21)  include  classical  plasticity  as  a  special  case. 
However,  there  is  an  important  difference.  It  can  be  easily  shown  that  the  current 
model  is  such  that  the  stress  tensor  a.,  is  noncoaxial  wilt,  the  inelastic  deformation 
rate  tensor  dp.  This  feature  furnishes  an  ability  to  model  localization  phenomena. 

A  number  of  special  case  studies  using  the  foregoing  results  are  under  current 
study.  These  investigations  are  intended  to  determine  the  modeling  capabilities  and 
limitations  of  the  slip-system  description  of  concrete  behavior. 
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Figure  J.2.  Influence  of  confinement  on  strength  and  ductility. 
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Figure  J.3.  Behavior  of  reinforced  concrete  in  compression. 
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